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Abstract 


The  realization  of  all-optical  devices  requires  identification  and  characterization  of 
materials  possessing  large  nonlinearities  and  low  loss.  Although  gallium  arsenide  — 
aluminum  gallium  arsenide  (GaAs-AlGaAs)  quantum  well  structures  have  very  large 
third  order  resonant  nonlinearities  just  below  the  bandgap,  the  absorption  in  this  region 
makes  them  unusable  for  all-optical  processing  devices.  Using  these  materials  at  below 
Vi  of  the  bandgap  reduces  the  losses,  but  also  reduces  the  nonlinearity. 

In  this  research,  we  investigate  a  different  method  of  accessing  the  nonlinearity.  By 
fashioning  GaAs-AlGaAs  double  heterostructure  devices  into  semiconductor  optical 
amplifiers,  a  strong  positive  nonlinearity  correlated  to  the  gain  of  the  device  was 
identified  and  characterized.  This  nonlinearity  is  all  the  more  valuable  as  these  devices 
can  provide  gain  to  the  signal  as  opposed  to  loss  thereby  boosting  the  effective 
nonlinearity  of  the  devices.  The  detriment  in  using  this  type  of  carrier  induced  effect  is 
that  the  nonlinearity  is  not  instantaneous  and  requires  over  1  ps  to  stabilize. 

Standard  interferometric  techniques  were  used  to  characterize  the  nonlinearity  in 
these  forward  biased  devices.  The  resulting  measurements  when  compared  to  Ken- 
nonlinearity  of  generic  form  +  ^  returned  an  effective  nonlinear  coefficient  of 

refractive  index  on  the  order  of  5xlO'*°cm^AV. 


Further  characterization  of  these  spatial  solitons  proved  that  they  were  capable  of 
guiding  a  second,  weaker  non-coherent  and  cross-polarized  signal  beam  within  the 
refractive  index  trail  created  by  the  soliton  writing  beam. 

Following  that,  a  soliton  beam-steerer  was  designed  and  fabricated  from  the  same 
semiconductor  double  heterostructure  material.  The  device  employed  a  single  mode 
waveguide  to  define  the  soliton  launch  path.  The  beam  then  passed  through  an  electro¬ 
optic  steering  prism  before  landing  within  a  40-pm  wide,  active  slab  waveguide  section. 
The  solitons  that  were  formed  within  this  last  section  of  the  device  were  readily  steered 
across  a  20-nm  wide  region  at  the  output  facet  of  the  slab  using  current  injection  into  the 
steering  prism.  The  mechanism  for  the  beam  steering  was  due  to  the  carrier-induced 
nonlinearity  and  our  measurements  confirmed  that  the  electrically  induced  nonlinearity 
was  of  the  same  magnitude  as  the  optically  induced  nonlinearity  (6.05x10^^  cm^/carrier). 
This  device  may  have  potential  applications  as  beam-routers  in  telecommunications 
networks. 
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CHAPTER  1:  INTRODUCTION 


Until  the  development  of  the  laser  by  Maiman  in  I960,*  the  field  of  optics  was  limited 
to  linear  systems.  At  this  time,  copper  cables  held  a  monopoly  on  long  haul 
communications,  and  the  cladded  core  fiber  optic  cable  had  just  been  introduced^  mainly 
for  the  purpose  of  image  transmission  through  a  bundle  of  glass  fibers.^ 

The  laser  provided  the  first  stable  source  of  intense,  coherent  light,  and  the  advent  of 
the  laser  launched  a  new  revolution  in  optics.  Soon  after  the  laser’s  development, 
Franken  demonstrated  second  harmonic  generation  for  the  first  time  using  crystalline 
quartz.'*  Over  the  next  two  years  optical  rectification,^  sum  frequency  generation,^  optical 
parametric  amplification,^  stimulated  Raman  scattering*  and  third  order  intensity 
dependent  refractive  index  changes^  were  all  reported. 

The  first  semiconductor  laser  was  reported  in  1962,  however,  it  wasn’t  until  the  later 
years  of  the  1970’s  that  improved  material  growth  combined  with  sophisticated  micro¬ 
fabrication  techniques  paved  the  way  for  the  mass  production  of  reliable  cw  GaAlAs 
laser  diodes.*®  Techniques  such  as  gain  switching,**’*^  proton  implanted  quantum  well 
passive  mode  locking,*^  active  mode-locking,*'*  and  hybrid  mode  locking*^  have  been 
utilized  to  shorten  output  pulses  and  raise  peak  powers  of  semiconductor  diode  lasers. 
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Pulses  as  short  as  200fs  with  over  160W  of  peak  power  have  been  demonstrated  using 
external  cavity,  semiconductor  laser  systems  with  intra-cavity  elements.*^ 

At  the  same  time,  improved  purity  and  fabrication  processes  reduced  transmission 
losses  in  optical  fibers  from  an  initially  high  value  of  >1000  dB/km  to  20dB/km  by 
1970,^^  and  to  0.2  dB/km  by  1979.**  Nonlinearities  in  fibers  such  as  stimulated  Raman 
and  Brillouin  scattering,*^  optically  induced  birefringence,^**  four-wave  mixing,^*  self¬ 
phase  modulation  and  temporal  solitons^^  were  demonstrated  in  fibers  throughout  the 
1970’s.  These  techniques  have  been  used  for  a  variety  of  systems  such  as  the  erbium- 
doped  fiber  amplifier  developed  in  1990  and  made  commercially  available  in  1992,^^  and 
in  fiber  pulse  compression  techniques  allowing  the  production  of  a  6  fs  pulse.^  Most 
recently,  low  loss  plastic  fibers  are  under  development  for  use  in  short  haul,  high  traffic 
applications.  These  links  permit  the  use  of  650  nm  laser  diodes  and  have  been 

25 

demonstrated  to  operate  in  excess  of  2.5GBps. 

1 . 1  DEVICES  FOR  INFORMATION  TECHNOLOGY 

It  should  be  no  surprise  that  research  pushes  the  envelope  of  optical  processing  speed. 
With  the  explosion  of  data  traffic  passing  through  land,  ocean,  and  satellite  links,  a  real 
need  exists  for  systems  that  pass  as  much  data  through  every  communications  link  as  fast 
as  possible.  Fiber  optic  lines  with  a  theoretical  effective  bandwidth  of  25THz^^  can 
replace  thousands  of  coaxial  links,  lOOMHz  bandwidth,  and  even  more  twisted  pair  lines, 
IMHz  bandwidth.*"  Unfortunately,  the  ability  of  electrical  components  to  load  25THz  of 
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information  onto  a  fiber  optic  link  has  not  caught  up  with  the  link’s  capability,  and  most 
likely  never  will.  Optical  systems,  on  the  other  hand,  have  the  potential  of  routing 
information  onto,  through,  and  off  of  these  high-speed  links.  Recently,  researchers  at  the 
NEC  Corporation  demonstrated  a  2.6  Terabit/sec  (20-Gbit/sec  x  132  channels) 
communication  link.^^  Devices  such  as  an  all-optical  clock  recovery  system**  and  a  20- 
channel  optical  transmitter  fabricated  firom  a  single  laser  diode  assist  in  loading  and 
unloading  information  fi'om  these  links.  High  speed  demultiplexers  such  as  the 
Terahertz  Optical  Asymmetric  Demultiplexer  (TOAD),^"  nonlinear  optical  loop  mirror 
(NOLM),**  and  soliton  switch**  can  optically  route  data  up  to  several  hundred  gigabits 
per  second.  A  prototype  Y-junction  device  has  also  been  demonstrated  to  provide  a 
switching  response  capable  of  demultiplexing  data  at  a  rate  of  588  GHz.**  Many  other 
switching  devices  using  the  nonlinear  effects  listed  above  have  been  fabricated  and 
demonstrated.*'* 

1.2  MATERIAL  NONLINEARITIES 

The  performance  of  optical  switches  depends  strongly  on  the  materials  they  are 
fabricated  fi'om.  Several  figures  of  merit  (FOM)  exist  as  ratings  of  a  material’s  suitability 
for  use  as  a  nonlinear  optical  device.  For  example,  the  FOM  for  a  linear  loss  system  is 


where  the  FOM,  W,  is  to  be  maximized,  nj  is  the  nonlinear  coefficient,  /(OJ  is  the  launch 
intensity,  oo  is  the  absorption,  and  A  is  the  wavelength.  As  all  nonlinear  optical  devices 
require  at  least  a  2k  nonlinear  phase  shift,^^  this  FOM  must  exceed  1.0.  The  form  of  the 
FOM  varies  based  upon  the  mechanism  responsible  for  the  nonlinearity,  however,  all 
maximized  nonlinear  FOM’s  have  two  things  in  common:  a  high  nonlinear  coefficient 
and  low  absorption. 

Vast  improvements  have  been  made  in  creating  materials  with  large  nonlinearities 
and  low  loss,  as  well  as  harnessing  the  nonlinearities  to  their  fullest  extent  possible. 

Franken  used  crystalline  quartz  («?„  =  4*10’”  for  his  second-harmonic 

generation  experiment.  Today,  lithium  niobate  (LiNbOa),  barium  titanate  (BaTiOa),  and 
lithium  iodate  (LilOa)  possess  second  order  nonlinear  coefficients  one  to  two  orders  of 
magnitude  larger  than  quartz.  The  nonlinear  coefficient  of  Cd(jeAs2  is  three  orders  of 
magnitude  larger  than  that  of  quartz.  Nonlinear  materials  have  also  been  developed  with 
outstanding  third  order  nonlinearities.  The  high-speed  conjugated  polymer  PTS  has  a 

third  order  nonlinearity  of  2.2x10’’^  and  an  unusually  high  fifth-order 

nonlinearity  of  -0.8x10'”  2  •  This  combination  permitted  researchers  to 

demonstrate  stable  2-dimensional  first  and  second  order  solitary  waves,  typically  very 
unstable,  within  this  material^^.  In  comparison,  silica  glass  fibers  have  a  third  order 

nonlinearity  of  only  2x10'”  •  The  extremely  low  losses  in  these  fibers,  however. 
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allows  very  efficient  length  scaling  and  help  maintain  Si02  fibers  as  a  viable  nonlinear 
material^^  for  all-optical  devices. 

Device  geometries  also  play  an  important  role  in  developing  nonlinear  devices  and 
techniques.  Second  harmonic  fields,  for  instance,  will  back  convert  to  the  primary 
wavelength  then  back  to  the  second  harmonic  in  a  periodic  manner.  Quasi-Phase 
matching  eliminates  this  back  conversion  by  stopping  or  reversing  the  down  conversion 
during  that  part  of  the  propagation.^*  Cascading  nonlinearities  provides  another  avenue 
to  advance  the  state  of  the  art.  Cascading  makes  use  of  the  second  harmonic  down 
conversion  to  create  an  effective  third  order  nonlinearity.  This  process  takes  advantage 
of  the  large  second  order  coefficient  to  generate  large  third  order  nonlinearities.  Using 
this  approach,  researchers  have  created  spatial  solitons  in  KTP  crystals.'*®’^*  This  method 
provides  stable  beam  self-focusing  and  second  harmonic  generation  together  in  one 
automatically  aligned  path.  Using  this  technique,  the  index  tuning  range,  usually  small 
for  typical  second  harmonic  crystals,  was  enhanced  by  a  factor  of  ten.'*^ 

1 .3  NONLINEARITY  IN  SEMICONDUCTOR  OPTICAL  AMPLIFIERS 

One  of  the  drawbacks  of  third  order  nonlinear  coefficients  is  that  they  are  typically 
small,  usually  1x1  O'**  tolxlO‘'^cyw^  /W .  Creating  an  index  change  large  enough  to 
form  spatial  solitons  in  these  materials  requires  intense  optical  intensities,  typically 
generated  by  mode-locked  or  Q-switched  short  pulse  lasers.  Materials  with  larger 
nonlinearities  have  been  found,  however,  they  typically  have  large  absorption.^^  This 


5 


absorption  degrades  device  performance  making  the  nonlinearity  an  interesting,  but 
useless,  phenomenon. 

Semiconductors  also  possess  large  nonlinearities  near  their  bandgap.^^’'*^  Third  order 

%45  * 

r  for  exciton  saturation  have  been  reported.  Like 

many  other  materials,  semiconductors  also  possess  a  large  absorption  near  the  bandgap, 
when  they  are  inactive.  One  potential  is  to  use  a  passive  AlGaAs  waveguide  at  below  Vx 
of  the  bandgap  to  access  the  Kerr  like  nonlinearity  while  avoiding  two-photon  absorption. 
At  this  wavelength,  the  positive  third  order  nonlinearity  was  measured  to  be 

l.SxlO'*^  and  1.43x10''^  TE  and  TM  polarization  states  relative  to 

the  slab  waveguide.'*^  As  the  ratio  of  cross-phase  modulation  to  self-phase  modulation  is 
almost  one,  several  varieties  of  spatial  solitons  such  as  primary,  vector  and  Manakov 
have  been  observed.'*'^  Unfortunately,  the  small  magnitude  of  this  nonlinearity  requires 
large  optical  powers  to  access  enough  self-phase  modulation  to  observe  self-focusing. 

There  exists,  however,  a  strong  resonant  nonlinearity  at  the  band  edge.  The 
nonlinearity  lies  well  within  the  absorption  edge  of  the  material,  but  by  fashioning  the 
material  into  a  semiconductor  optical  amplifier  the  nonlinearity  can  be  used  in  the  gain 
region  of  the  material.'**  ■ 

There  are  many  mechanisms  involved  in  this  gain  nonlinear  process.  Figure  1.1 
shows  the  time  evolution  probe  response  of  the  nonlinearity  to  a  440  femto-second  pump 
pulse  for  GaAs-AlGaAs  amplifier  in  the  gain  (a),  transparency  (b),  and  absorption  (c) 
regions.  In  each  of  these  regions,  a  strong  instantaneous  negative  nonlinearity,  lasting  on 
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the  order  of  150  fs,  was  attributed  to  rapid  electronic  or  virtual  processes."*^  A  positive 
going  nonlinear  spike  caused  by  free  carrier  heating^”’*'  and  by  spectral  hole  burning 
within  the  homogeneously  broadened  amplifier^^’’^  quickly  replaced  this  nonlinearity. 
The  long-term  behavior  of  the  probe  phase  shift  depends  upon  the  slowly  recovering 
carrier  density  change  caused  by  the  pump  beam.^°  Both  the  Burstein-Moss  shift  and 
Free  Carrier  Plasma  models  describe  this  shift.  The  sign  of  this  long-term  nonlinearity 
depends  on  whether  the  device  is  in  the  absorption,  transparency,  or  gain  region. 


-50  5  10 

Pump-Probe  Delay  (psec) 


Figure  1.1  Time  Resolved  Probe  Phase  Change  in  AlGaAs  Optical  Amplifiers  for  the 
(a)  gain,  (b)  transparency,  and  (c)  absorption  regimes.  Pump  duration  was  440  fs^°. 
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These  measurements  indicate  that  the  short  pulse  bright  soliton  is  doomed  to  fail,  as 
the  nonlinearity  is  always  negative  for  short  pulses.  Delfyett,  however,  performed  the 
same  measurements  using  longer  pulses  and  found  that  the  negative  going  initial  spike 
was  not  detectable  for  pulses  of  one  picosecond  in  duration.*®  A  long,  one-microsecond 
pulse  therefore,  was  selected  to  create  the  spatial  solitons  characterized  in  this 
dissertation.  The  two  mechanisms  leading  to  a  positive  nonlinearity  are  the  Burstein- 
Moss  Shift,  and  the  index  shift  caused  by  the  change  in  the  free  carrier  plasma.  Both 
provide  a  positive  nonlinearity  in  the  spectral  gain  region. 

1.3.1  DYNAMIC  BURSTEIN-MOSS  SHIFT 

Absorption  and  gain  in  a  semiconductor  is  the  process  of  either  raising  an  electron 
from  the  valence  band  to  the  conduction  band  or  vice  versa  through  interaction  with  a 
photon.  In  optical  absorption,  the  energy  in  a  photon  is  transferred  to  an  electron  in  the 
valence  band  thus  raising  the  energy  level  of  the  electron.  As  there  are  no  levels  to  be 
occupied  between  the  bands,  photons  with  energies  below  the  bandgap  are  not  absorbed 
and  pass  through  the  material.  Photons  with  energies  exceeding  the  bandgap  are  readily 
absorbed.  Upon  entering  the  conduction  band,  the  electron  will  rapidly  thermalize  in 
times  much  less  than  the  carrier  lifetime  to  fill  states  near  the  band  extrema.  As  only  one 
carrier  can  occupy  a  state  at  any  time,  the  next  absorption  event  must  raise  an  electron 
from  a  slightly  lower  level  in  the  valence  band  to  a  slightly  higher  level  in  the  conduction 
band.  Although  the  effect  is  negligible  for  a  single  carrier,  a  significant  amount  of 
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absorption  will  cause  the  band  edge  to  shift  towards  the  higher  energy  side  (blue  shift). 
This  is  the  Burstein-Moss  shift.*'*’”  The  reverse  shift  is  also  true  in  the  emission  region. 
When  a  distribution  of  excited  carriers  in  the  conduction  band  is  significantly  depleted  by 
stimulated  gain,  the  band  edge  shifts  towards  lower  energies  or  a  red  shift. 

A  periodic  crystal  lattice,  such  as  a  semiconductor,  has  distinct  properties  that  effect 
the  number,  distribution,  and  energy  of  the  electrons  and  lack  of  electrons  (holes)  in  the 
ground  state  (valence  band)  of  the  material.  Characteristic  of  semiconductors  is  the  band 
diagram.  The  diagram  depicts  a  series  of  closely  spaced  states  at  and  below  the  valence 
band  for  holes,  and  at  and  above  the  conduction  band  for  electrons.  Separating  these 
bands  is  the  forbidden  bandgap  where  there  are  no  states.  The  number  of  carriers,  n  for 
electrons  in  the  conduction  band  and  p  for  holes  in  the  valence  band,  is  characterized  by 
the  expression*^ 


«,=\Na)fa)‘k 
p>  =  \mA^-fa)\k 


(1.2) 


N(^  is  the  density  of  states  function  that  defines  the  number  of  states  available  for 
potential  occupation,  and/(^  is  the  Fermi  function  which  defines  the  probability  that  a 
given  state  will  be  occupied.  Under  normal  equilibrium  conditions,  the  Fermi  function  is 
dependent  on  the  Fermi  energy  ^  and  temperature  T  and  is  identical  for  both  n  and  p 
concentrations. 
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Nc  and  Nv  are  the  effective  density  of  states  constants  for  the  semiconductor.  For 
GaAs  =  4.4xl0"c/w“^  and  =  8.2x10"cot”^  .  The  rate  at  which  absorption  or  gain 

can  occur  is  dependent  upon  the  probability  that  states  exist  and  are  filled  that  satisfy 
conservation  of  momentum.  These  rates  are  defined  by  the  expressions 
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)VA<  )NA> 


(1.6) 


In  these  expressions,  P  is  the  transition  probability  and  Np  is  the  density  of  photons  of 
energy  4,  and  Nj  is  the  joint  density  of  states  function.  Weighted  energies  and  %2  are 
selected  to  maintain  conservation  of  energy,  conservation  of 

momentum,  The  photon  momentum  is  negligible  compared  to  the 

mass  of  the  electron  and  hole,  and  we  approximate  this  conservation  as  k  electron  ^  ^hole  ’ 

Assuming  the  energy  bands  are  parabolic  with  characteristic  masses  me*  and  n^  . 

* 

4.  =5. 

(17) 


m*  is  the  system  reduced  effective  mass  given  by  — r  =  — r  +  — r 

ntr  m, 

The  net  stimulated  emission  rate  per  energy  is  the  difference  between  the  absorption 
rate  and  the  emission  rate.  Dividing  this  difference  by  the  net  number  of  photons 

c  • 

crossing  a  given  area  of  speed  —  and  assuming  Pem  and  Pst  are  identical,  the  spectral 


form  of  the  gain  g(^)  becomes 


g(V  = 


•J2(  ml)‘^q^p^cr 
3m,  Go 


(1.8) 
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In  an  active  semiconductor,  the  carrier  distribution  is  created  through  electrical 
injection  of  carriers  into  the  gain  region.  Knowing  the  number  of  carriers  injected  into 
the  active  gain  region,  and  the  average  lifetime  the  carriers  exist  prior  to  recombination, 
we  can  calculate  the  quasi  Fermi  levels  and  finally  the  gain  spectrum.  If  stimulated 
emission  occurs,  the  average  carrier  density  will  decrease  locally  causing  a  change  in  the 
gain  spectrum. 

This  change  can  be  directly  related  to  the  change  in  the  index  of  refraction  through 
the  Kramers-Kronig  relationship,  appendix  B.  Figure  1.2  shows  the  effect  that  a  change 
in  the  gain  spectrum  has  on  the  refi'active  index  of  the  material,  assuming  an  incident 
optical  excitation  at  the  peak  of  the  gain  curve.  It  is  important  to  note  that  the  change  in 
index  is  not  the  same  for  all  wavelengths,  but  that  the  change  in  refractive  index  is 
positive  for  all  energies  below  the  peak  of  the  gain.  Furthermore,  the  change  in  the  index 
of  refraction  calculated  by  this  method  is  a  change  in  index  per  depleted  carrier  (An/AN). 
Thus,  the  actual  change  in  the  index  of  refraction  is  proportional  to  the  change  in  carrier 
population  within  the  amplifier. 
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Figure  1 .2  Index  Change  Due  to  Absorption  Change,  Kramers-Kronig  Relationship. 
Inset  shows  gain  change  and  delta  gain  in  the  region  about  the  peak  of  the  gain. 


1.3.2  FREE  CARRIER  PLASMA 

The  effect  of  free  carriers  on  a  semiconductor’s  susceptibility  results  from  the 
characteristic  absorption  of  free  carriers  in  the  infrared.  The  peak  of  the  absorption 
increases  in  energy  with  carrier  density,  thus  increasing  the  carrier  density  will  decrease 
the  susceptibility,  and  the  index  of  refraction,  below  the  peak  of  the  absorption.  The 
susceptibility’s  dependence  on  carrier  concentration  is  given  as^^’’* 

X  =  — ^  0-9) 

€o 
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where  N  is  the  carrier  concentration,  e  is  the  electron  charge,  e©  is  the  free  space 
permittivity,  o  is  the  angular  frequency  of  the  radiation,  and  me*  is  the  effective  carrier 
mass.  The  effective  mass  of  the  electrons  is  much  smaller  than  that  of  the  holes,  thus 
only  the  change  in  the  susceptibility  caused  by  the  electrons  need  be  taken  into  account. 
The  susceptibility  is  related  to  the  index  of  refraction  by  the  approximate  equation,  see 
appendix  A, 

n^=l  +  %  (110) 


Using  a  change  in  the  carrier  concentration  to  produce  a  change  in  the  index  of 
refraction  yields  the  perturbation  equations 


An  =  -^^AN 


2n^  dN 


An  =  - 


871^  Sp 


(1.11) 


AN  in  an  optical  excitation  case  would  be  the  change  in  carriers  caused  by  the  optical 
beam. 


1 .3.3  CARRIER  CONCENTRATION  CHANGE  IN  A  GAIN  DEVICE 


The  changes  in  the  index  of  refraction  derived  above  depend  upon  a  change  in  the 
carrier  concentration.  This  change,  however,  is  not  constant  in  an  amplifier  and  must  be 
taken  into  account.  To  account  for  this  change,  we  examine  the  change  in  carrier 
concentration  in  a  differential  increment  along  the  optic  axis.  The  power  along  the  axis  is 
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given  by  P(z)  =  Pf„  exp( gz )  where  g  is  the  gain  coefficient,  z  is  the  distance  along  the 

optic  axis  with  z=0  at  the  input,  and  Pi„  is  the  optical  power  at  the  input.  The  change  in 
the  carrier  concentration  in  an  incremental  segment  is  proportional  to  the  increase  in 
optical  power  found  by  differentiating  the  power  equation. 

dP 

^  =  gP  (1-12) 

dz 

The  change  in  carrier  concentration  is  measured  within  the  volume  Adz,  where  A  is 
the  effective  cross  sectional  area  of  the  beam.  Power  relates  to  individual  electronic 
transitions  by  the  optical  quanta  ho,  thus  the  rate  of  carrier  depletion  per  unit  volume  is 

^  =  (1.13) 

dt  Atia 

Noting  that  the  power  per  unit  area  is  intensity,  and  the  change  in  carrier  density 
relaxes  with  a  time  constant  t,  the  change  in  the  carrier  concentration  becomes 

AN  =  ^  (1.14) 

fi(o 

In  the  semiconductor  slab  waveguides  used  in  this  research,  the  cladding  material’s 
bandgap  exceeds  the  energy  of  the  gain  region,  and  is  transparent.  Thus  any  of  the  gain, 
or  absorption  will  be  restricted  to  the  region  of  the  gain.  Thus,  a  smaller  effective  gain 
needs  to  be  used  as  opposed  to  the  calculated  gain.  This  effective  gain  is  the  calculated 
gain  times  the  confinement  factor,  F,  which  is  the  ratio  of  the  amount  of  guided  energy  in 
the  active  region  over  the  total  amount  of  guided  energy.  Applying  equation  to  the  basic 
forms  of  the  change  in  the  index  of  refraction,  equations  and  yields  the  standard  form 
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Aw  =  Bgl ,  where  B  is  the  effective  change  in  the  index  of  refraction  in  the  zero-gain 
transparency  region.  Using  the  standard  Kerr  model  of  the  nonlinearity  w  =  Wq  +  nj ,  the 


nonlinear  coefficient,  w^ ,  is  effectively  amplified  by  the  gain  provided  by  the  device.  As 
the  gain  of  a  semiconductor  amplifier  can  reach  values  of  10,000  m  ’,  the  possibility 
exists  to  access  nonlinearities  four  orders  of  magnitude  larger  than  those  typically  seen  in 
passive  media. 

Although  this  method  indicates  an  instantaneous  w^,  an  alternate  method  to 
examine  the  gain  uses  the  total  phase  shift  to  calculate  the  nonlinear  coefficient.  This 
total  phase  shift  is  found  by  the  integration  of  the  incremental  phase  shifts  along  the  path 
of  propagation.  This  becomes 


=  Win 


e^-1 

g 


(1.15) 


This  equation  has  the  same  form  as  the  non-gain  equation  if  the  top  and  bottom  of  the 
right  hand  side  are  multiplied  by  the  distance  z  and  a  new  nonlinear  coefficient  of 

'  -1 

»,  - t/i,  (1.16) 

gZ 


is  used.  For  values  of  the  gz  multiple  less  than  2,  this  function  exhibits  linear  behavior 
with  a  slope  of  approximately  V2.  Thus  with  the  length  of  the  device  fixed,  the  gain 
effects  the  nonlinearity  by  approximately  W2(l+‘/4  gz). 
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1.3.4  GAIN  SATURATION 


When  the  optical  beam  traversing  the  amplifier  becomes  so  large  that  the  rate  of 
removal  of  excited  carriers  exceeds  the  electrical  pumping  rate,  the  gain  will  saturate. 
The  gain  spectrum  arising  from  the  radiative  recombination  of  excited  electrons  of  the 
semiconductor  optical  amplifier  is  homogeneously  broadened.  Thus,  gain  saturation  in 
these  devices  h,  is  given  by  the  following,’’ 

g  =  (1,17) 

1  +  — 

Here  go  is  the  small  signal  gain.  As  the  third  order  optical  nonlinearity  depends  on  the 
magnitude  of  the  gain,  the  nonlinearity  itself  eventually  saturates  with  a  characteristic 
saturation  intensity  4. 

1.3.5  EXCITON  EFFECTS 

Exciton  saturation  creates  large  nonlinearities  as  well.'*’  Excitons  are  a  lower  energy 
state  created  when  electrons  and  holes  on  either  side  of  the  bandgap  are  in  close  spatial 
proximity  and  retain  a  coulombic  attraction  to  each  other.  Excitons  are  not  typically 
measured  at  room  temperature  in  bulk  devices;  however,  high  quality,  quantum  well 
structures  provide  enough  spatial  confinement  to  form  excitons.  The  exciton  can  be  seen 
as  one  or  two  absorption  peaks  just  below  the  main  absorption  band  edge. 
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Promoting  carriers  into  the  conduction  band  easily  saturates  excitons.  The  saturation 
alters  the  absorption  spectrum,  and  creates  an  intensity-based  change  in  the  refractive 
index  through  the  Kramers-Kronig  relationship.  As  quantum  structures  were  not 
available  at  the  time  of  this  research,  exciton  saturation  was  not  investigated. 

1.4  SOLITON 


To  describe  a  spatial  soliton,  we  start  with  the  polarization  of  the  material  the  soliton 
exists  in.  The  expression  for  the  polarization  of  a  linear  material  is  P  =€„  where  x 

is  the  linear  response  of  the  material  to  a  field  E.  In  nonlinear  materials,  however,  the 
polarization  varies  in  time,  depends  on  frequency,  and  contains  multiple  cross  terms 
linking  the  polarizability  between  frequency  components.  To  allow  mathematical 
manipulation  of  the  polarization  term,  it  is  common  to  assume  that  the  polarization  is 
instantaneous  and  to  approximate  the  polarization  response  by  a  Taylor  expansion  of  the 
polarization. 


P(o)=Go 


Y  + x'^^  Y®.‘  ®  1 .  ®  2  )P(  ®  1  )P(  ^2)+^ 


(1.18) 


In  this  notation,  x^^\(»;ou0)2)  is  the  second  order  susceptibility  at  frequency  w  arising 
fi'om  stimulation  at  frequencies  ©i  and  ©2.  Of  particular  interest  are  situations  where 
there  is  only  one  input  stimulus,  and  we  are  interested  only  in  the  response  of  the  material 
at  that  frequency.  This  situation  is  called  the  self-action  nonlinearity  and  equation  above 
reduces  to 
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P(a>)=^o  +  +...)  (1.19) 

this  term  can  be  inserted  into  the  nonlinear  Schroedinger  Wave  equation  (NSWE), 
described  in  appendix  A,  to  derive  the  closed  form  solution  for  a  bright  spatial  soliton. 


assuming  the  term  is  negligible. 


I(x)  =  A^sech^ 


r  ^ 


X 


(1.20) 


Here  Ao  measures  the  amplitude  of  the  soliton,  and  ao  is  the  waist  of  the  soliton  in  the 
dimension  where  diffraction  is  allowed.  As  stated  above,  the  soliton  is  usually  stable  in 
only  one  dimension,  the  confinement  in  the  other  dimension  must  be  provided  by  a  slab 
waveguide.  In  this  dimension,  where  diffraction  is  not  allowed,  w  is  used  to  describe  the 


beam  size.  Inserting  equation  (1.2)  into  the  NSWE,  the  relationship 


Ps^0  = 


2w 


(1.21) 


is  derived,  where  Ps  is  the  linear  power  contained  in  the  soliton  and  no  and  ni  are  the 
linear  and  nonlinear  indices  of  refraction.  These  form  the  expression  for  the  index  of 
refraction  in  a  Kerr  nonlinear  material. 

n(x)  =  na+  nj(  x)  (1 .22) 


Noting  that  all  of  the  terms  to  the  right  of  the  equality  are  constants,  the  ratio  between 
soliton  power  and  waist  remains  constant.  Therefore,  in  low  loss  systems,  the  waist 
remains  constant. 
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The  soliton  can  be  described  logically  in  this  manner.  Propagating  light  diffracts. 
Because  the  index  of  refraction  in  a  Kerr  like  material  varies  with  the  intensity  of  light  at 
a  given  position,  the  index  of  refraction  of  a  propagating  beam  will  be  highest  at  the 
center  of  the  beam  and  the  least  at  the  edges.  This  forms  a  positive  lens  causing  the  beam 
to  focus.  For  any  beam,  the  tighter  the  focus  generally  means  the  larger  the  diffraction. 
For  the  nonlinearity,  the  larger  the  optical  intensity  at  the  peak,  the  tighter  the  focusing 
will  be.  Eventually,  a  stable  condition  can  be  reached  in  one  dimension  and  the  beam’s 
natural  diffraction  is  exactly  balanced  by  the  nonlinear  self-focusing. 

1.4.1  MATHEMATICAL  VS  OTHER  SOLITONS 

The  spatial  soliton  described  in  the  previous  section  is  a  closed  form  solution  to  the 
NSWE  in  a  perfect  Kerr  nonlinear  material.  This  type  of  soliton  has  many  specific 
mathematical  properties  such  as  when  N  solitons  collide,  N  solitons  will  emerge  which 
travel  in  the  same  direction  as  the  input  solitons,  but  endure  a  lateral  shift.  Additionally, 
mathematical  solitons  are  reciprocal  meaning  that  a  series  of  beams  launched  backwards 
along  the  path  used  to  create  a  forward  pattern,  will  recreate  the  launched  pattern  at  the 
input.  Many  other  soliton  like  beams  can  be  created  that  possess  the  self-guided  nature  of 
a  perfect  mathematical  soliton  but  may  not  necessarily  possess  all  of  the  characteristics  of 
a  mathematical  soliton.  Throughout  the  remainder  of  this  dissertation,  the  term  soliton 
will  be  used  to  label  a  self-focused  beam;  however  the  beam  may  not  conform  to  the 
strict  properties  of  the  mathematical  soliton. 
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1.4.2  DEMONSTRATED  SOLITONS 


Kerr  nonlinearity  below  the  two-photon  absorption  bandgap  was  used  to  create 
fundamental  bright  solitons.  Vector  solitons,  and  Manakov  solitons  in  AlGaAs  slab 
waveguides.  This  was  remarkable  as  very  few  materials  have  the  self  and  cross  phase 
modulation  terms  nearly  equal  as  AlGaAs  does.^ 

In  second  order  nonlinear  materials,  the  second  order  nonlinearity  can  cascade  to 
produce  the  fundamental  wavelength  and,  if  the  nonlinearity  is  strong  enough,  can  result 
in  soliton  like  self-guiding.  As  discussed  above,  this  has  been  demonstrated  in  KTP.®* 
Khitrova  and  Gibbs  demonstrated  spatial  solitons  in  InP,  using  the  same  mechanism 
used  in  this  research  to  produce  spatial  solitons  in  GaAs.®^ 

Two-dimensional  solitons  have  been  created  in  materials  with  a  strong  negative  fifth 
order  nonlinearity.  This  nonlinearity  prevents  the  beam  from  focusing  too  tightly  thereby 
avoiding  damage  to  the  optical  medium. 

Dark  solitons  can  be  generated  in  materials  with  negative  third  order  nonlinearities, 
and  photoreffactive  materials  can  also  self-guide  optical  beams.  All  of  these  alternate 
methods  of  creating  soliton  like  self-guided  beams,  with  the  exception  of  the  pure  Kerr 
material,  have  some,  but  not  all,  of  the  properties  of  the  mathematical  soliton.  As  such, 
the  term  soliton  will  be  used  throughout  the  remainder  of  this  dissertation  with  the 
understanding  that  the  solitons  measured  in  this  work  are  not  pure  mathematical  solitons. 
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1 .4.3  SOLITON  BEAM  STEERING  AND  TRAPPING 


The  profile  of  the  index  of  refraction  creates  a  waveguide  as  long  as  the  soliton  exists. 
This  waveguide  can  be  used  to  guide,  steer,  and  trap  other  beams.  The  form  of  the  index 
profile  of  the  soliton  waveguide  for  Kerr  like  materials  will  be  sech^  in  shape.  The  index 
shift  is  not  typically  large,  however,  the  graded  index  nature  of  the  waveguide  permits 
better  guiding  without  a  large  shift  in  the  index  of  refraction.  Kang  used  solitons  to  trap 
beams  in  AlGaAs  waveguides,  and  Lawrence  used  2  dimensional  solitons  to  steer  beams 
through  a  nonlinear  crystal.  A  prime  concern  in  passing  a  soliton  through  a  steering 
prism  is  the  soliton’s  ability  to  remain  intact.  Kang  studied  this  problem  and  concluded 
that  solitons  remained  stable  even  when  Vi  of  the  soliton  passed  through  the  prism  area 
while  the  other  half  did  not.^'* 

1 .5  STUDY  TO  FORM  SOLITONS  IN  SOA  AND  TRAPPING/GUIDING 
OF  SECOND  BEAM 

As  shown  in  the  previous  sections,  a  positive  third  order  optical  nonlinearity  of 
significant  magnitude  exists  in  the  spectral  region  of  the  semiconductor  gain  for  long 
pulses.  This  dissertation  undertook  the  task  of  creating  the  conditions  necessary  to 
demonstrate  spatial  solitary  waves  using  this  nonlinearity,  characterizing  the  nature  of  the 
nonlinearity  and  the  resulting  spatial  solitons,  and  finally  using  these  solitons  to  guide  a 
signal  beam.  These  experiments  form  the  basis  for  continued  work  on  a  proposed  single 
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stage  optic  switch  using  electro-optic  steering  prism  to  steer  the  information,  and  soliton 
waveguides  to  contain  the  information. 

The  work  carried  out  in  this  dissertation  was  performed  on  an  undoped  bulk  GaAs 
gain  region  sandwiched  between  n  and  p  type  AlGaAs  clad  regions.  Although  quantum 
well  devices  could  potentially  provide  larger  third  order  nonlinearities,  the  wafer  foundry 
was  not  able  to  undertake  the  task  of  optimizing  an  optical  quantum  well  structure. 

1.6  OUTLINE  OF  TfflS  DOCUMENT 

This  dissertation  comprises  6  chapters.  Chapter  1  outlines  the  basic  theory  behind  the 
generation  of  third  order  nonlinearities  in  semiconductor  optical  amplifiers.  Chapter  2 
discusses  the  fabrication  details  involved  in  fabricating  the  devices  used  throughout  this 
research.  Chapter  3  describes  the  experiments  preformed  to  characterize  the 
semiconductor  material  and  confirm  the  positions  and  magnitudes  of  nonlinearities. 
Chapter  4  reports  on  the  successes  achieved  in  creating  and  characterizing  a  spatial 
soliton  created  by  end  fire  coupling  an  input  beam  into  the  cleaved  end  of  a  wide-stripe, 
active,  semiconductor  waveguide.  Chapter  5  builds  upon  the  soliton  measurements  by 
guiding  a  second  signal  beam  along  with  the  soliton  beam.  Finally,  chapter  6  reviews  the 
important  findings  discovered  throughout  this  research  and  lays  out  recommended 
additional  experiments  advancing  this  research  towards  the  ultimate  goal  of  creating  a 
usable  soliton  matrix  switch. 
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CHAPTER  2:  DEVICE  STRUCTURE  AND  FABRICATION 


Although  the  optical  and  electrical  characterization  of  the  materials  and  devices  is 
a  prime  concern,  the  ability  to  reproducibly  fabricate  these  devices  was  an  integral 
portion  of  this  work.  Particularly  so  as  a  new  clean  room  facility,  along  with  totally  new 
equipment,  required  installation,  checkout,  calibration,  and  application.  As  such,  this 
chapter  will  outline  the  calculations  and  procedures  used  to  design  the  optoelectronic 
devices  used  throughout  this  work.  In  addition,  critical  process  steps  used  to  fabricate 
these  devices  will  also  be  explained. 

2.1  DEVICE  STRUCTURE 

2.1.1  ACTIVE  SEMICONDUCTOR  OPTICAL  GAIN 

As  shown  in  the  introduction,  for  photon  energies  near  the  bandgap,  the  self-focusing 
nonlinearity  in  a  semiconductor  optical  amplifier  exists  only  in  the  spectral  range  where 
the  semiconductor  possesses  gain.  Although  the  gain  can  be  optically  induced,  this 
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research  employed  electrically  injected  carriers  to  create  the  population  inversion 
necessary  to  provide  gain. 

Several  requirements  need  to  be  fulfilled  for  optical  gain  to  be  realized.  First,  the  k- 
selection  rule  must  be  met,  i.e.  the  momentum  of  the  electrons  and  holes  must  match  for  a 
radiative  recombination  event  to  occur.  Second,  impurities  that  result  in  deep  level  traps 
within  the  forbidden  energy  gap  must  be  minimized  because  they  will  compete  for  the 
excited  electrons.  Third,  the  photons  that  are  emitted,  as  result  of  radiative 
recombination,  must  be  confined  to  the  gain  region  containing  a  high  density  of  inverted 
carriers.  Finally,  this  inversion  is  obtained  by  electrical  injection  of  electrons  into  the 
conduction  band  of  a  p-type  material  that  contains  a  high  population  of  vacancies  in  the 
ground  state  (holes). 

The  first,  second,  and  third  requirements  are  easily  satisfied  through  proper  material 
selection  and  device  design.  The  k-selection  rule  is  met  by  using  direct  bandgap 
semiconductors,  in  our  case  gallium  arsenide  (GaAs).  The  second  requirement  is  met  by 
using  a  p-doped  material  in  the  gain  region  for  greatest  gain  efficiency.  Injection  of 
minority  electrons  into  the  p-type  gain  region  ensures  high  probability  of  radiative 
recombination  since  there  is  already  an  abundance  of  holes  The  third  requirement  of 
confining  the  created  photons  to  the  gain  region,  is  fulfilled  by  waveguiding  the  optical 
wave  within  the  gain  medium.  This  means  that  the  gain  region  must  possess  a  higher 
index  of  refraction  than  the  adjoining  regions.  This  index  is  adjusted  in  the  vertical 
dimension  by  substituting  aluminum  for  gallium  during  the  crystal  growth  to  form 
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AlxGa/.xAs  where  x  is  the  molar  fraction  of  aluminum  in  the  crystal.  The  resulting  index 
of  refraction  is  related  to  the  volume  fraction  (x)  by  the  empirical  formula 

w("x;  =  3.590-0.710x  +  0.091x^  (2.1) 

Samples  used  in  this  study  contained  Al^Ga^.^cAs  cladding  layers  that  had  x=0.20. 
Confinement  in  the  horizontal  dimension  can  either  be  done  by  directly  fabricating  a 
fixed  index  profile  across  the  device  or  by  generating  a  nonlinearly  induced  spatial  index 
profile  with  a  self-focused  wave.  To  satisfy  the  final  requirement  of  having  a  spatially 
localized  carrier  inversion,  a  forward  biased  injection  diode  is  used.  Additionally,  the  use 
of  the  higher  bandgap  AlGaAs  for  the  cladding  layers  has  the  effect  of  confining  the 
electrons  within  a  potential  well  located  across  the  intrinsic  region.  The  resulting 
structure  and  refractive  index  profile  is  shown  in  Figure  2.1.  This  type  of  device  is  called 
a  double  heterostructure  device  due  to  the  n  and  p  doped  Alo^Gao.sAs  along  with  the 
intrinsic  GaAs  region.  The  gain  region  in  the  actual  structure  was  grown  undoped 
because  the  optical  quality  of  undoped  GaAs  is  far  superior  to  that  of  p-doped  GaAs.  On 
the  other  hand  the  efficient  diffusion  of  holes  across  the  p-doped  AlGaAs  layer  ensures 
that  the  gain  medium  effectively  contains  a  large  population  of  holes  for  efficient 
recombination. 
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G2-2329  Structure 


bidex  Profile 


GaAs  Contact  Layer  p  doped 


Alo^GaAs  Clad  Layer  p  doped 
1.0  urn  thick 


GaAs  Active  Region  0.2um 


AL^GaAs  Clad  Layer  n  doped 
1.5  umfiiick 


GaAs  Buffer  n  doped  0.3um 


Substrate  n  doped 
62Sum  thick 


^neff=  3.4943 


Figure  2.1  Cross  section  of  the  semiconductor  optical  amplifier.  Right  graph  shows 

index  profile  of  stack. 


To  describe  the  operation  of  this  device,  we  must  compare  the  band  diagrams  of  the 
static  and  electrically  pumped  diode.  When  n  and  p  regions  of  a  diode  are  connected 
through  a  metallurgical  junction,  the  holes  in  the  p  region  diffuse  to  the  n  side  and 
electrons  in  the  n  region  diffuse  to  the  p  side  leaving  a  depletion  region,  void  of  fi'ee 
carriers,  between  the  two  materials.  This  process,  however,  leaves  unbalanced  donor  ions 
or  net  positive  charges  in  the  n  region  and  unbalanced  acceptor  ions  or  net  negative 
charge  in  the  p  region.  This  resulting  space  charge  field  restricts  further  diffusion  of 
holes  and  electrons  and  equilibrium  results.  The  addition  of  an  intrinsic  region  between  n 
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and  p  layers  extends  the  length  the  depletion  region.  This  condition  is  depicted  in  Figure 
2.2a.  When  an  external  electric  field  is  applied  to  the  device,  the  depletion  region  grows 
for  reverse  bias  or  shrinks  for  forward  bias.  It  is  assumed  that  the  externally  applied 
electric  field  is  dropped  exclusively  across  the  intrinsic  region.  When  a  large  enough 
forward  bias  is  applied,  the  region  becomes  very  small,  and  holes  and  electrons  easily 
difiuse  across  the  junction  and  recombine.  This  condition  is  called  the  forward  biased 
condition,  shown  in  Figure  2.2b.  As  more  forward  bias  is  applied,  the  bands  continue  to 
move  such  that  the  field  is  dropped  across  the  intrinsic  region,  the  device  eventually 
enters  a  degenerate  state.  Figure  2.2c.  In  this  state,  an  excess  of  electrons  exists  locally 
with  an  excess  of  holes,  creating  an  inversion. 


2.1.2  OPTICAL  CONFINEMENT 

To  design  an  efficient  optical  amplifier,  the  beam  propagation  in  the  slab  (vertical) 
dimension  must  be  determined.  Treating  the  amplifier  as  a  symmetric  slab  waveguide  in 
the  TE  polarization  state,  we  set  the  boundary  conditions  across  the  gain  region  and 
infinity  points  then  solve  the  continuity  equations  to  derive  the  transcendental  equation 
relating  the  effective  index  and  size  of  the  gain  region. 

tan(  idi  +  nm)=  ^  (2.2) 

K  -  ps 
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where  h  is  the  actual  thickness  of  the  guiding  layer,  m  is  the  order  of  the  propagated 
mode,  K^=nf\vac^  -  s^=p^-nc^kvac^,  and  p^=p^-ns\yac^.  Within  these  equations,  nf  is  the 


Fermi  Level 


I - ► 

Distance 


V=0  volts 


Inversion 

I  Region 

•« - ► 


Distance 


V=Small  Forward 
Bias 


V=Large  Forward 
Bias 


n doped 

GaAs 

p doped 

Alo^GaAs 

Gain 

AIo^GaAs 

Figure  2.2  Band  Diagram  of  a  Double  Heterostructure  Amplifier.  Top  diagram  (a)  no 
bias.  Middle  diagram  (b)  small  forward  bias  and  bottom  diagram  (c)  large  forward  bias. 
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index  of  refraction  of  the  guiding  layer,  nc  is  the  index  of  the  top  cladding  layer,  ns  is  the 
index  of  the  lower  cladding  layer,  kvao  is  the  vacuum  wave  vector  of  the  optical  wave,  and 
p  is  the  propagation  constant  of  the  wave  inside  the  waveguide,  p  relates  to  the  effective 
index  of  refraction  of  the  waveguide  by  neif=P/kvac.  To  ensure  a  single  mode  device,  a 
thickness  must  be  selected  that  is  below  the  minimum  thickness  required  for  the  m=2 
mode.  Additionally,  to  achieve  the  largest  nonlinear  action  requires  maximizing  the 
optical  intensity  within  the  active  region  by  making  the  guided  region  as  small  as 
possible,  while  also  maximizing  the  overlap  of  the  nonlinear  region  with  the  guided  wave 
by  making  the  guided  region  as  large  as  possible.  Maximizing  these  divergent 
requirements  is  accomplished  by  minimizing  the  effective  thickness  given  by 

A  .  =A  +  i+i  (2,3) 

s  p 

where  h,  s,  and  p  are  as  defined  previously.  For  the  parameters  used  to  design  the 
devices  used  in  this  research,  this  formula  minimizes  at  a  thickness  of  0.2  microns. 

Although  lateral  confinement  by  solitary  wave  propagation  is  studied  in  Chapters 
4  and  5,  device  characterization  outlined  in  Chapter  3  required  a  two-dimensionally 
confined  single  mode  waveguide  design.  The  single  mode  confinement  in  the  vertical 
dimension,  as  described  above,  is  accomplished  by  fabricating  the  appropriate  AlGaAs  - 
GaAs  layers  during  the  MBE  growth.  Confinement  in  the  horizontal  direction  is 
achieved  by  delineating  a  ridge  waveguide.  The  propagation  index  of  this  structure  is 
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found  by  using  the  effective  index  method.  In  this  method,  infinite  slab  waveguides  are 
assumed  for  each  of  the  propagation  regions,  the  effective  index  is  calculated  for  each, 
then  a  slab  mode  in  the  orthogonal  direction  is  analyzed  using  these  calculated  refractive 
indices.  The  index  profile  for  the  ridge  waveguide  is  formed  by  partially  etching  the  top 
AlGaAs  clad  layer  from  the  surface  of  the  wafer.  The  deeper  the  etch,  the  lower  the 
index  of  refraction  for  the  area  below  the  etch.  The  effective  index  in  each  of  these  slabs 
can  be  calculated  by  altering  the  boundary  conditions  for  the  top  layer.  Instead  of 
assuming  the  field  exponentially  decays  to  zero  in  the  top  clad  layer,  we  allow  some  of 
the  field  to  enter  the  regions  above  the  cladding  layer.  An  alternate  method  of  calculating 
the  effective  index  in  multiple  layer  systems  is  to  evaluate  the  transmission  using  ABCD 
matrices  through  a  multiple  layer  stack  at  differing  angles.  Propagation  modes  appear  as 
resonant  spikes  in  the  transmission.  Using  this  method,  the  index  of  refraction  vs  etch 
depth  of  our  one  micron  upper  cladding  is  shown  in  Figure  2.3.  Computer  code  used  for 
this  method  appears  in  appendix  C. 
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Top  Clad  Thickness  (Angstroms) 

Figure  2.3  Dependence  of  effective  index  on  top  clad  layer  thickness.  ns=3.45164, 
nf=3.59  df=0.2  urn,  nc=3 .45164  thickness  varies,  top=air. 

Using  these  design  constraints,  a  wafer  was  grown  using  a  Varian  360  Molecular 
Beam  Epitaxy  (MBE)  chamber  at  the  Air  Force  Wright  Laboratories,  Dayton  OH. 
Structure  as  grown  is  shown  in  Figure  2.1.  The  buffer  layer  shown  in  this  figure  helps 
prepare  the  wafer  surface  for  high  quality  epitaxial  growth  and  reduces  dislocation 
diffusion.  This  layer  has  no  optical  function. 

2.2  FABRICATION 

The  devices  used  in  this  research  were  fabricated  at  CREOL  using  our  class  100- 
cleanroom  facility.  An  overview  of  the  steps  used  for  fabrication  is  shown  in  Figure  2.4. 
Each  of  the  steps  was  developed  for  this  clean  room  and  a  compendium  of  the  fabrication 
follows  below.  The  first  step  (not  needed  for  the  devices  used  in  self-guiding 
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Ni/Ge/Au  p  type  Otimic  Metal 
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Figure  2.4  Fabrication  Steps  for  Device  Soliton  Device 


experiments)  involved  the  formation  of  a  ridge  waveguide  by  wet  chemical  etching. 
Following  this,  a  p-type  top  ohmic  contact  was  applied  using  thermal  vacuum 
evaporation.  In  some  devices,  current  isolation  trenches  were  required,  and  these  were 
formed  by  wet  chemical  etching  on  either  side  of  the  p-type  ohmic  contact.  Next  the 
sample  was  thinned  down  to  200  microns  to  facilitate  cleaving.  After  thorough  cleaning, 
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an  n-type  ohmic  contact  was  applied  to  the  entire  backside  of  the  wafer.  After  this 
metallization,  the  wafer  was  cleaned  and  rapid  thermal  annealed  to  form  the  eutectic  alloy 
ohmic  contacts.  Finally,  the  sample  was  cleaved  to  form  clean  facets  for  end-fire  optical 
coupling  and  mounted  on  a  copper  post  using  conductive  epoxy  to  enable  an  electrical 
contact  to  the  n  side  of  the  sample  and  to  permit  efficient  removal  of  heat  generated  in  the 
active  amplifier.  Each  of  the  tools  used  in  this  fabrication  procedure  is  described 
separately  below. 

2.2.1  PHOTOLITHOGRAPHY 

Regions  for  different  processing  were  defined  using  photoresist  and 
photolithography.  In  this  process,  a  photoactive  polymer  liquid,  the  photoresist,  is 
applied  to  the  sample  and  spun  at  high  speeds  to  forma  thin  uniform  layer.  The  layer  is 
cured,  then  selectively  exposed  to  UV  light  using  contact  masks  and  a  Karl  Suss  MJB-3 
mask  aligner.  The  MJB-3  delivers  25  mW/cm^at  405  nm  (G  line)  and  18  mW/cm^  at  365 
nm  (H  line).  The  layer  is  developed  (selectively  removed)  using  a  liquid  developer.  This 
process  results  in  precise  openings  in  a  resistant  polymer  photoresist  coating.  These 
openings  define  the  sections  of  the  wafer  affected  by  subsequent  processing.  The 
photoresists  used  in  this  work  are  as  follows. 

Shipley  1400-17:  This  is  a  thin  positive  photoresist  best  suited  for  high  precision 
work  to  fabricate  structures  with  feature  sizes  approaching  1  micron.  It  is  best  suited  for 
wet  chemical  etching,  and  performs  poorly  for  metal  liftoff  The  photoresist  thickness  is 
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4500A  using  a  3000-rpm  spin  speed  for  30  seconds.  The  photoresist  is  partially  cured 
(soft-baked)  at  100°C  for  5  minutes  on  a  hot  plate.  The  optimum  exposure  time  with  the 
MJB-3  is  1.6  sec  followed  by  a  30  second  develop  in  one  part  Shipley  351  developer 
mixed  with  five  parts  water.  After  successfiil  developing,  the  resist  must  be  hot  plate 
baked  for  5  minutes  at  100°C  for  more  complete  curing.  The  resist  may  delaminate 
without  this  post  bake,  especially  in  wet  chemical  etches. 

Shipley  1823:  This  is  a  medium  positive  photoresist  best  suited  for  moderate  precision 
work  to  realize  structures  with  dimensions  of  2-3  microns  and  larger.  It  is  the  best 
general-purpose  photoresist,  well  suited  for  ohmic  contact  metal  liftoffs.  However,  it 
may  not  be  suitable  for  liftoff  of  thick  metals  or  metal  plating.  The  resist  thickness  is  1 .3 
microns  using  a  3000-rpm  spin  speed  for  30  seconds.  It  cures  at  119°C  for  5  minutes  on 
a  hot  plate.  Exposure  time  with  the  MJB-3  is  3.5  seconds  followed  by  a  30  second 
develop  in  one  part  Shipley  351  developer  mixed  with  six  parts  water.  After  successful 
developing,  the  photoresist  must  be  hot  plate  baked  (post-baked)  for  5  minutes  at  119°C. 
The  resist  may  delaminate  without  the  post  bake,  especially  in  wet  chemical  etches.  This 
resist  is  made  by  mixing  24  parts  Shipley  1827  resist  with  5  parts  Shipley  type  p  thinner, 
then  gently  mixing  for  5  minutes.  Appropriate  care  must  be  taken  to  ensure  that  mixing 
does  not  induce  bubbles  into  the  polymer  liquid. 

Shipley  1827:  This  is  a  thick  positive  photoresist  that  is  useful  for  thick  metal  liftoffs 
and  plating.  This  resist  is  unsuitable  for  high  precision  work.  A  photoresist  thickness  of 
2.5  microns  is  obtained  using  a  3000-rpm  spin  speed  for  30  seconds.  Soft  baking  is 
performed  at  119°C  for  5  minutes  on  a  hot  plate.  The  optimum  exposure  time  is  8 
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seconds  with  the  MJB-3.  Develop  and  post-bake  are  identical  to  the  Shipley  1823 
photoresist. 

Edge  Bead.  A  problem  with  using  partial  wafers  is  the  formation  of  photoresist  edge 
beads.  Spinning  any  photoresist  on  a  wafer  without  beveled  edges  causes  an  edge  bead  to 
form.  The  bead  is  several  microns  in  thickness  above  the  applied  photoresist  and  extends 
for  10-20  microns  from  all  edges.  The  bead  impedes  the  direct  contact  of  the  center  part 
of  the  wafer  and  the  photomask  causing  fuzzy,  undefined  shadow  exposures.  This  results 
in  poor  edge  definitions  in  structures  approaching  one  to  two  microns.  The  effect  is  not 
as  detrimental  to  structures  exceeding  five  microns.  Although  mechanical  tools  exist  to 
remove  edge  beads,  these  are  not  available  at  CREOL.  The  beads  can  easily  be  removed, 
however,  if  positive  photoresist  is  used.  To  remove  the  bead,  the  center  portion  of  the 
wafer  is  covered  and  the  edges  are  overexposed  for  60  seconds.  Because  positive 
photoresist  is  not  removed,  or  effected,  unless  exposed,  the  wafer  can  be  developed  using 
standard  developer.  After  developing,  the  wafer  is  rinsed  in  deionized  water  for  one 
minute  and  dried.  Hot  plate  baking  is  not  required  between  this  step  and  subsequent 
exposure-develop  cycles. 

2.2.2  WET  CHEMICAL  ETCHING 

After  masking  the  substrate,  the  exposed  areas  can  be  thinned  by  using  wet  chemical 
etching.  To  be  effective  in  forming  the  waveguides  used  in  this  research,  the  wet 
chemical  solution  must  uniformly  etch  both  GaAs  and  AlGraAs  and  must  not  be  selective 
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in  crystal  plane  orientation.  The  solution  used  was  one  part  phosphoric  acid  (H3;P04)  at 
1  molar  concentration,  two  parts  hydrogen  peroxide  (H2O2)  in  a  30%  solution  with  water, 
and  50  or  100  parts  deionized  water.  The  etch  rate  in  p  doped  GaAs  and  AlGfaAs  was 
1850A/sec  for  the  1:2:50  solution,  and  llOOA/sec  for  the  1:2:100  solution.  For  best 
results,  the  solution  was  stirred  after  combining  the  constituents  and  allowed  to  stabilize 
for  five  minutes  prior  to  use.  Etching  solutions  were  not  re-used  after  four  hours. 
Temperature  variations  also  impacted  etch  rates  significantly,  and  the  etch  rates  reported 
here  are  for  solutions  at  23°C. 

2.2.3  PLASMA  PROCESSING 

In  plasma  processing,  materials  can  either  be  removed  or  deposited  through  the  use  of 
plasma  exposure.  Plasma  processing  enables  material  removal  or  deposition  not 
available  by  other  means  and  in  geometries  not  available  by  other  means.  For  example, 
organic  polymers  can  be  removed  selectively  leaving  nearly  a  90°  sidewall.  Or,  silicon 
nitride  can  be  deposited  across  the  entire  surface  of  a  semiconductor  wafer  regardless  of 
the  shapes,  geometries,  or  materials  on  the  surface  of  the  wafer.  Both  processes  were 
accomplished  in  a  PlasmaTherm  790  dual  chamber  plasma  processing  system.  In  this 
research,  the  plasma  processing  was  used  in  three  processing  steps:  oxygen  descum, 
deposition  of  the  dielectric  layer  silicon  nitride,  and  selective  removal  of  silicon  nitride. 
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2.2.3.1  OXYGEN  DESCUM 


Ideally,  after  a  photolithographic  step  where  openings  are  created  in  an  organic 
photomask,  the  openings  are  totally  clear  of  any  organic  residue.  In  practice,  however,  a 
very  thin  film,  usually  less  than  100 A  thick,  of  organic  residue,  or  scum,  frequently 
remains.  The  film  is  not  visible  under  a  microscope,  but  the  film  causes  metals  to 
delaminate  after  vacuum  metal  deposition,  or  non-uniform  wet  chemical  etches  due  to 
areas  of  organic  scum.  Removal  of  this  material  is  termed  a  descum  step. 

Various  methods  exist  to  remove  this  film  including  a  dip  in  very  dilute  ammonium 
hydroxide  or  exposure  of  the  material  to  a  weak  oxygen  plasma.  The  ammonium 
hydroxide  method  potentially  can  damage  the  substrate  and  is  very  isotropic,  meaning 
that  the  openings  in  the  photoresist  will  become  slightly  distorted  as  the  ammonium 
hydroxide  removes  material  from  the  side  walls  of  the  openings  as  well  as  from  the  open 
areas  themselves. 

The  oxygen  plasma  descum  method  overcomes  these  shortcomings.  First,  the  weak 
oxygen  plasma  is  reactive  only  to  organic  materials,  thus  inorganic  semiconductors, 
dielectrics,  and  metals  are  not  affected.  Second,  plasma  chemistry,  pressure,  and  power 
combinations  exist  that  are  very  anisotropic  meaning  that  material  is  removed  from  flat 
surfaces,  but  not  from  sidewalls.  Even  though  material  is  removed  from  all  flat  surfaces, 
including  the  areas  that  should  be  covered,  the  removal  represents  only  a  small  fraction  of 
the  total  resist  thickness.  In  a  typical  descum  step,  100 A  of  organic  photoresist  is 
removed  from  a  layer  that  started  at  between  SOOOA  and  ISOOOA  thick. 
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Organic  removal  by  oxygen  plasma  requires  that  oxygen  radicals,  created  in  a  plasma, 
contact  the  surface  of  the  organic  material  and  form  volatile  compounds  such  as  carbon 
dioxide  and  water.  In  the  PlasmaTherm  790  system,  13.8  MHz  RF  excitation  between 
two  parallel  plates  in  the  plasma  chamber  excites  the  plasma.  A  DC  bias  applied  between 
the  plates  accelerates  the  plasma  radicals  towards  the  lower  plate  where  the  samples  are. 
As  the  acceleration  is  directionally  vertical,  etch  profiles  tend  to  be  anisotropic  in  the 
vertical  direction.  The  addition  of  a  light  carrier  gas,  such  as  helium,  to  the  plasma 
controls  the  oxygen  radical  concentration  and  aids  in  the  mechanical  removal  of  materials 
thereby  increasing  the  anisotropic  nature  of  the  etch.  Chamber  pressure,  gas 
concentration,  and  power  primarily  control  etch  rates,  while  gas  flow  rate  controls 
uniformity.  Figure  2.5  shows  the  characterization  of  three  oxygen  etch  setups  on  organic 
photoresist.  In  this  flgure,  the  descum  process,  used  as  a  mild  removal  of  trace  organics 
prior  to  metalizations  or  etches,  is  the  slowest  etch.  The  anisotropic  etch  process  is 
designed  as  an  anisotropic  removal  of  photoresist.  After  etching  through  1.5  microns  of 
photoresist  with  the  anisotropic  process,  no  discemable  increase  in  opening  size  was 
measured.  Finally,  the  ash  process  is  designed  as  a  cleaning  or  blanket  removal  of  all 
organic  materials  on  the  surface  of  a  substrate. 
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Figure  2.5  Etch  Rate  of  Three  Polymer  Etch  Processes:  Descum,  Anisotropic®*,  and 

Ash  Processes. 


2.2.3  2  SILICON  NITRIDE  DEPOSITION 

The  devices  used  in  this  research  require  ohmic  contact  only  over  the  active  region, 
typically  the  size  of  the  waveguide;  four  to  50  microns.  This  ohmic  contact’s  size  makes 
electrical  contacting  nearly  impossible.  Placing  a  large  metal  contact  pad  and  metal 
current  pathway  over  the  ohmic  contact  permits  current  flow  and  contacting,  however, 
any  metal  contact  directly  to  the  p-side  of  the  wafer  causes  current  leakage.  To  overcome 
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this,  a  dielectric  layer  is  inserted  between  the  contacting  pads  and  substrate.  This  layer  is 
created  by  blanket  coating  the  entire  surface  with  the  dielectric,  then  selectively  removing 
the  dielectric  above  the  ohmic  contacts.  Following  this,  the  metal  contact  pads  can  be 
deposited  without  contacting  the  substrate. 

Two  processes  were  designed  to  create  the  dielectric  layer.  The  first  was  to  spin  coat 
a  polymer  glass  suspended  in  a  carrier  solvent,  spin  on  glass,  to  the  wafer.  The  carrier 
solvent  was  evaporated  by  hot  plate  baking  the  wafer  at  100°C  for  5  minutes,  leaving  an 
amorphous  gel  on  the  surface  of  the  wafer.  The  wafer  was  then  heated  to  350°C  for  30 
minutes  to  cure  the  gel  into  a  hard,  inorganic  silicon  dioxide  cap  layer.  Although  spin  on 
glasses  from  two  venders  were  tested  using  a  variety  of  application  techniques,  the  silicon 
dioxide  layer  created  using  this  technique  delaminated  from  the  wafer  during  the  final 
thermal  processing  stages  of  ohmic  contact  formation. 

The  second  method  proved  successful.  In  this  technique,  the  substrate  is  patterned 
with  waveguides.  Prior  to  patterning  the  top  ohmic  contact,  a  uniform  layer  of  silicon 
nitride  was  applied  to  the  wafer  using  a  PlasmaTherm  790  chemical  vapor  deposition 
system.  In  the  system,  silane,  Sillt,  and  ammonia,  NH3,  gasses  combine  in  a  plasma. 
Silicon  radicals  combine  with  nitrogen  radicals  forming  silicon  nitride,  Si3N4. 


3SiH4+4NH3  -^SijN^+llH, 


(2.4) 


Although  the  silicon  nitride  deposits  on  the  first  surface  it  contacts,  the  substrates  are 
heated  to  250°C  to  help  create  uniform  coatings  and  relieve  deposition  stress. 

Adjusting  the  ratio  of  silane  to  ammonia  can  control  the  index  of  refraction  in  the 
resulting  film  from  1.85  to  2.2.  Figure  2.6  shows  the  relationship  between  index  and  the 
silane  to  ammonia  ratio.  Control  of  the  refractive  index  is  useful  in  antireflective 
coatings,  but  does  not  provide  practical  control  for  dielectric  insulators.  Silicon  nitride 
films  grown  in  this  research  have  a  ratio  of  approximately  1.9. 


Figure  2.6  Dependence  of  Refractive  Index  on  Growth  Temperature  and 
Silane/ Ammonia  Ratio^ 

The  deposited  silicon  nitride  films  coat  uniformly  to  most  geometries,  and  adhere  to 
all  materials  used  in  this  fabrication  including  photoresist,  semiconductors,  and  metals. 
Figure  2.7  shows  an  electron  beam  microscope  photo  of  a  typical  conformal  silicon 
nitride  coating  shown  as  a  dark  band,  along  the  side  of  a  waveguide  ridge. 
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Figure  2.7  Electron  micrograph  of  a  uniform  silicon  nitride  layer  coating  the  sides  of 
a  laser  diode  ridge.  Silicon  nitride  layer  (arrows)  appears  as  a  dark  black  layer  between 

the  substrate  and  top  contact  metal.®^ 


Deposition  rates  are  dependent  on  RF  plasma  power.  Uniformity  depends  on 
deposition  and  flow  rates.  Figure  2.8  shows  the  deposition  rates  of  two  silicon  nitride 
deposition  processes.  The  rate  increases  from  74  A/min  for  the  20-watt  plasma  to  157 
A/min  for  the  100-watt  plasma. 

Residual  stress  in  the  film  can  be  changed  from  tensile  to  compressive  by  adding 
varying  amounts  of  helium  to  the  plasma.  Tensile  stress  in  films  deposited  without 
helium  at  13.56  MHz  is  in  the  range  of  2.5x10*^  dynes/cm^.  By  increasing  the  ratio  of 
helium  to  nitrogen  in  the  gas  mixture  used  to  deposit  silicon  nitride,  the  film  stress  can  be 
controlled  from  tensile,  through  zero,  to  compressive.  For  a  100-watt  plasma,  a  ratio  of 
30%  nitrogen  and  70%  helium  will  produce  a  no  stress  film.^* 
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Facilities  to  measure  the  index  of  refraction  and  stress  of  these  films  was  not  available 
at  the  time  this  work  was  preformed,  thus,  the  system  was  not  characterized  for  these 
parameters. 


Figure  2.8  Deposition  Rate  Dependence  on  Plasma  Power 


22.3.3  SILICON  NITRIDE  ETCHING 

Once  a  uniform  silicon  nitride  layer  was  deposited  over  the  substrate,  a  photomask  as 
patterned  over  the  nitride  and  the  nitride  removed.  This  exposed  the  area  requiring  ohmic 
contacts  to  the  surface  of  the  sample,  thus  the  photomask  used  to  pattern  the  silicon 
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nitride  was  also  used  to  pattern  deposited  ohmic  metal  onto  the  substrate.  This  is  termed 
a  self-aligned  process. 

Wet  chemical  and  plasma  etching  processes  were  studied  to  measure  etch  rates.  In 
the  wet  chemical  procedure,  buffered  oxide  etch  consisting  of  hydrofluoric  acid  and 
sodium  fluoride  was  diluted  with  deionized  water  1:5.  Samples  were  immersed  in  the 
etching  solution,  then  rinsed  with  deionized  water  for  two  minutes.  This  method 
produced  etch  rates  of  95lA/min  for  the  20-watt  plasma  film,  and  1012A/min  for  the 
100-watt  plasma  film. 

Plasma  etching  was  preformed  using  a  mix  of  92%  Freonn  and  8%  oxygen  cleaning 
gas  in  a  20-watt  plasma.  The  etch  rate  for  the  20-watt  plasma  film  was  41 1  A/min.  The 
etch  rate  for  the  100- watt  film  was  461  A/min. 

2.2.4  METALLIZATION 

Electrical  contact  is  made  to  the  samples  by  selective  deposition  of  a  series  of  metals 
that,  when  annealed,  combine  with  the  surface  of  the  semiconductor  and  form  a  layer 
without  a  potential  barrier.  The  contacts  exhibit  only  small  resistance  and  are  called 
ohmic  contacts.  The  contact  used  for  p-type  GaAs  is  50A  gold:  200A  zinc  followed  by 
2000A  of  gold.  The  contact  used  for  n-type  GaAs  is  50A  nickel:  200 A  germanium,  and 
2000A  gold.  The  samples  were  annealed  only  once  at  420°C  for  120  seconds  in  a  rapid 
thermal  annealer.  Non-ohmic  contact  pads  were  deposited  after  the  anneals,  and  were 
comprised  of  50A  titanium:  200A  platinum:  2000A  gold.  In  these  contacts,  the  titanium 
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formed  an  adhesion/wetting  layer,  the  platinum  formed  a  gold  diffusion  barrier,  and  the 
gold  was  the  high  conduction  contact  layer. 

To  form  good  ohmic  contacts,  surface  preparation  was  critical.  All  photoresist  and 
oxides  must  be  removed  from  the  area  of  deposition.  To  remove  photoresist,  an  oxygen 
plasma  descum  should  be  used.  This  comprises  the  exposure  of  the  sample  to  a  low 
power  oxygen  plasma  to  remove  50-lOOA  of  surface  polymer  uniformly  across  the  wafer. 
Following  the  descum,  the  surface  oxides  are  removed  by  dipping  the  wafer  in  a  solution 
of  one  part  buffered  oxide  etch  with  50  parts  water  for  30  seconds.  This  gentle  oxide  etch 
removes  surface  oxides  formed  by  exposure  of  the  wafer  to  atmosphere  and  during  the 
descum  procedure.  After  the  acid  etch  dip,  the  wafer  is  rinsed  for  2  minutes  in  a  flowing 
deionized  water  bath,  then  thoroughly  blow  dried  with  nitrogen.  The  sample  is  then 
immediately  loaded  into  the  metallization  chamber  without  further  heating.  Heating  will 
form  additional  oxides. 

Following  metallization,  the  sample  was  soaked  in  acetone  for  an  hour  to  dissolve  the 
photoresist  on  the  surface  of  the  wafer.  This  caused  the  metal  covering  the  photoresist  to 
peal  off  the  surface  of  the  wafer  leaving  only  the  metal  in  the  exposed  areas.  The  edges 
of  the  remaining  metal  pads  possessed  “wings”  formed  when  the  metal  broke  from  the 
metal  over  the  photoresist.  Moving  the  sample  to  a  clean  beaker  of  acetone  and  using 
ultrasound  for  1  to  5  minutes  removed  most  of  the  “wings”.  A  Cole-Parmer  ultrasonic 
cleaner  was  used  because  it  provided  a  lower  ultrasound  intensity  than  other  units  such  as 
the  Fisher  ultrasonic  cleaner  that  has  been  found  to  delaminate  the  metal  pads  and  on 
occasions,  even  shatter  samples. 
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2.2.5  THINNING  AND  CLEAVING 


In  this  research,  optical  beams  were  launched  into  the  devices  using  the  end-fire 
coupling  technique.  End-fire  coupling  requires  a  highly  polished  and  flat  interface. 
Fortunately,  the  use  of  crystalline  semiconductors  provides  an  easy  method  of  forming 
these  near-perfect  optically  flat  facets,  namely  cleaving.  In  the  cleaving  process,  a  small 
guide  nick  or  short  scribe  mark,  about  Vi  millimeter,  is  made  on  the  fi’ont  of  the  wafer 
perpendicular  to  the  edge  of  the  substrate.  The  direction  of  the  nick  must  be  nearly 
parallel  to  a  crystal  plane.  Application  of  point  pressure  immediately  behind  the  nick,  on 
the  backside,  causes  the  wafer  to  split  along  the  nearest  primary  crystal  plane,  starting  at 
the  nick  point.  Cleaving  the  wafer  prior  to  fabrication  permits  the  alignment  of  processed 
structures  to  the  cleavage  planes. 

During  fabrication,  wafers  require  mechanical  strength  to  avoid  breakage.  As  such, 
wafers  are  typically  625  microns  thick.  This  is  advantageous  for  keeping  the  wafer 
whole  during  fabrication,  however  it  makes  cleaving  difficult.  In  order  to  get  better 
cleaves,  wafers  are  thinned  to  200  microns  after  the  majority  of  processing  is  completed. 
Once  thinned  to  this  thickness,  wafers  are  very  fragile  and  must  be  handled  delicately. 
The  first  step  in  the  thinning  process  is  to  mount  the  wafer  face  down  to  a  lapping  block 
using  an  acetone  soluble  wax  called  crystal  bond.  Measuring  the  starting  thickness  of  the 
wafer  and  lapping  flat  combination  provides  the  reference  point  to  determine  the  extent 
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of  the  thinning.  Then  the  lapping  flat  is  mounted  to  a  stabilized  lapping  flxture,  which 
holds  the  flat  parallel  to  the  lapping  surface.  This  fixture  rides  on  a  South  Bay 
Technology,  Model  920  lapping  table  covered  with  a  slurry  of  10  micron  alumina  power 
and  water.  Removing  400  microns  of  material  requires  about  20  minutes,  and  the 
thickness  is  checked  regularly.  Once  the  thickness  nears  225  microns,  the  lapping  flat  is 
removed  fi-om  the  fixture,  and  the  substrate  is  cleaned.  The  10  micron  lapping  slurry 
leaves  a  rough,  flat,  gray  surface  comprised  of  scratches.  These  scratches  interfere  with 
the  cleaving  and  are  removed  by  manual  polishing  using  finer  grit  material.  The  wafer- 
fixture  combination  is  moved  in  a  figure  8  motion  along  a  polishing  pad  covered  with  a 
slurry  of  0.3  micron  alumina  in  a  mixture  of  1  part  Clorox  bleach  with  10  parts  water. 
The  surface  of  the  wafer  becomes  shinny  and  free  of  scratches  within  10  minutes.  The 
wafer  is  thoroughly  rinsed  with  water,  then  soaked  in  acetone  until  the  wafer  separates 
from  the  lapping  flat. 

At  this  point  the  exposed  n-type  backside  of  the  wafer  is  mounted  on  a  carrier  glass 
slide  and  metallized  with  an  n-type  ohmic  contact.  After  metallization,  the  wafer  is 
detached  from  the  glass  slide,  and  thoroughly  cleaned.  To  establish  the  substrate  ohmic 
metal  contact,  the  wafer  is  rapid  thermal  annealed  for  120s  at  420®C,  which  completes  the 
fabrication  process.  Using  the  cleaving  process  described  above,  the  sample  is  cut  to  the 
required  length. 
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2.2.6  MOUNTING 


The  amount  of  electrical  current  passing  through  these  devices  causes  a  significant 
thermal  load.  Removing  the  heat,  as  well  as  making  electrical  contact  to  the  backside  of 
the  samples,  requires  secure  mounting  onto  an  electrically  and  thermally  conductive  post. 
These  posts  were  fashioned  from  1mm  thick  copper  sheets  and  milled  to  the  required 
dimensions.  In  cases  where  a  wide  sample  was  used,  the  mount  was  formed  by  joining 
two  sheets  together  using  thermally  cured  conductive  epoxy.  In  both  situations,  once  the 
platform  was  fabricated,  a  thin  coat  of  Epotech  H20E  epoxy  was  applied  to  the  mounting 
area  of  the  stud,  and  the  sample  was  carefully  placed  into  this  epoxy  ensuring  that  the 
epoxy  did  not  wick  up  onto  the  facets  of  the  sample.  Heating  the  stud/sample  to 
for  30  minutes  permanently  cured  the  epoxy  and  fixed  the  sample  to  the  stud. 

Electrical  contact  to  the  pad  on  the  p-side  of  the  sample  was  done  using  0.3  mill 
(50pm)  diameter  gold  wire.  These  wires  were  attached  simultaneously  to  the  device  and 
an  insulated  standoff,  which  was  mounted  to  the  copper  stud,  using  the  same  conductive 
epoxy  used  to  attach  the  sample  to  the  stud.  A  pair  of  tweezers  mounted  on  a  3-axis 
positioner  held  the  wire  in  place  while  the  epoxy  was  heat-cured.  To  provide  uniform 
current  injection  into  the  device,  multiple  wires  were  used  for  the  p-side  contact.  The 
device-stud  combination  could  then  be  mounted  into  the  test  fixture  for  evaluation. 
Figure  2.9  depicts  a  cross  section  of  a  wide  device  mounted  on  a  T-stud  and  attached  to 
the  water-cooled  testing  fixture. 
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Figure  2.9  Photograph  of  a  mounted  sample  containing  16  devices,  n  side  contact  is 
made  through  copper  standoff,  p  side  contacts  are  made  with  thin  wires  and  insulated 
bonding  pads  shown  to  either  side  of  the  sample. 


CHAPTERS:  CHARACTERIZATION 


Determining  a  wafer’s  characteristics  is  important  not  only  to  avoid  wasting  time 
fabricating  devices  from  material  that  is  unsuitable  for  testing,  but  also  to  provide  the 
necessary  parameters  to  analyze  experimental  results.  Additionally,  the  growth  of 
semiconductor  materials  remains  as  much  art  as  science,  and  each  set  of  wafers  grown  at 
different  times  will  invariably  have  different  characteristics.  Optical  and  electrical 
characterization  techniques  are  designed  to  determine  these  variables. 

3 . 1  PHOTOLUMINESCENCE  CHARACTERIZATION 

To  achieve  the  best  gain,  our  devices  required  a  well-grown  GaAs  intrinsic 
region,  free  of  major  defects  and  compositional  alterations,  between  the  n  and  p  doped 
AlGaAs  cladding  regions.  To  quantify  the  quality  and  growth  composition,  as  well  as  the 
spectral  dependence  of  the  gain  peak  with  carrier  injection  of  this  intrinsic  layer,  the 
photoluminescence  and  electro-luminescence  techniques  were  used. 
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3.1.1  PHOTOLUMINESCENCE 


In  the  photoluminescence  measurement,  excited  carriers  create  electrons  in  the 
conduction  band  and  holes  in  the  valence  band  of  the  semiconductor,  using  an  incident 
pump  beam  with  photon  energy  greater  than  the  bandgap  of  the  semiconductor.  In  this 
case,  a  3  mW  HeNe  laser  emitting  at  632.8-nm  wavelength  was  used  as  the  excitation 
source.  The  electrons  in  the  material  absorb  the  incident  pump,  and  are  raised  in  energy 
into  the  conduction  band.  As  (thermal)  phonon  transitions  within  the  conduction  band 
occur  1000  times  faster  than  the  radiative  decay  of  these  electrons  due  to  band-to-band 
transitions,  electrons  excited  above  the  band  edge  rapidly  settle  to  the  base  of  the 
conduction  band  and  then  recombine  with  holes  thereby  creating  photons.  Vibrational 
broadening  causes  a  Gaussian  energy  distribution  around  the  bandgap  energy  with  the 
peak  remaining  at  the  bandgap.  The  number  of  carriers  excited  into  the  conduction  band 
stays  relatively  small,  thus  band  renormalization  experienced  with  high  current  injection 
is  negligible.  This  allows  an  unperturbed  measurement  of  the  bandgap  energy.  The 
unperturbed  bandgap  of  GaAs  at  room  temperature  is  1.425  eV  corresponding  to  a 
wavelength  of  870  nm;  quantum  effects  are  not  expected  with  a  gain  region  width  of  0.2 
microns. 

The  photoluminescence  photons  emitted  from  the  intrinsic  gain  region  pass  through 
the  AlGaAs  clad  layers,  which  are  transparent  to  the  incident  photons  because  the 
AlGaAs  bandgap  is  larger  than  the  photoluminescence  photon  energy.  The  photons  that 
are  emitted  from  the  sample  in  a  Lambertian  distribution  are  collected  and  focused  into 
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the  entrance  slit  of  a  spectrometer.  A  schematic  of  this  experiment  is  illustrated  in  Figure 
3.1. 


Figure  3.1  Photoluminescence  Measurement 


Although  the  emitted  photons  from  the  active  region  will  pass  out  of  the  material 
without  being  absorbed  by  the  Alo  jGao.gAs  upper  cladding  region,  the  inbound  632  nm 
pump  laser  will  be  absorbed  and  may  not  penetrate  to  the  gmn  region.  Additionally, 
because  the  Alo.2Gao.8As  upper  cladding  region  does  not  absorb  the  photons  that  are 
emitted  from  the  GaAs  active  region,  it  will  efficiently  absorb  the  excitation  photons 
provided  by  the  632nm  HeNe  laser  and  prevent  them  from  reaching  the  active  region.  To 
enhance  the  photoluminescence  measurement,  the  top  GaAs  contact  layer  and  50%  of  the 
AlGaAs  cladding  layer  were  removed  using  wet  chemical  etching. 
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Using  a  three-milliwatt  HeNe  laser  to  excite  electrons  in  the  gain  region,  the 
photoluminescence  spectrum  shown  in  Figure  3.2  was  obtained.  The  peak  of  the 
spectrum  lies  at  868  nm  as  would  be  expected'  theory  predicts  870  nm.  To  determine 
how  the  spectrum  behaves  at  higher  temperatures,  the  photoluminescence  spectrum  was 
re-measured  with  sample  heated  to  80°C. 


Figure  3.2  Photoluminescence  Spectra  of  G2-2329  at  Room  Temperature  (22°C),  left 

trace,  and  at  80°C,  right  trace. 


54 


As  shown  in  Figure  3.2,  the  peak  of  the  photoluminescence  spectrum  of  the  heated 
sample  lies  at  883nm.  Semiconductor  bandgaps  follow  the  Varshni  empirical  fit  equation, 
equation  3.1. 
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(3.1) 


where  the  Varshni  fit  parameters  for  GaAs  have  been  entered  into  this  equation  3.1. 
At  80°C,  this  equation  predicts  that  the  photoluminescence  signal  should  peak  at  886nm. 
This  is  in  good  agreement  with  our  measured  value  of  883 nm.  These  measurements 
confirm  that  the  material  within  the  gain  region  is  good  quality  intrinsic  GaAs. 

3.1.2  ELECTRO-LUMINESCENCE 

A  second  analysis  technique  is  electro-luminescence.  Electro-luminescence 
characterizes  the  spectral  nature  of  the  semiconductor  amplifier’s  gain.  In  this  technique, 
the  output  of  a  laser  diode  fabricated  from  the  material  under  characterization  is 
spectrally  analyzed  at  varying  injected  carrier  levels  as  shoAvn  in  Figure  3.3.  The  earners 
are  injected  electrically  into  the  gain  region  of  the  diode,  and  recombine  with  the  holes 
injected  from  the  opposite  side  thereby  creating  photons. 
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Figure  3.3  Electro-Luminescence  Measurement  Configuration 


Under  forward  biased  conditions,  the  quasi  Fermi  levels  shift  towards  their  respective 
bandedges,  and  eventually  become  degenerate  by  crossing  the  bandedge.  Using  double 
heterostructure  devices  enhances  this  increase  in  population  by  physically  restricting  the 
diffusion  of  carriers  in  the  vertical  dimension  through  the  used  of  higher  energy  bandgap 
AlGa  As  layers  on  either  side  of  the  active  region. 

The  degeneracy  resulting  fi’om  the  respective  quasi-Fermi  levels  approaching  the 
bands,  causes  an  effect  called  bandtailing.  Bandtailing  bends  the  parabolic  band  edges 
into  the  bandgap.  This  results  in  a  narrowing  of  the  bandgap  and  a  concentration  of 
carriers  below  the  expected  bandedge.  Studies  done  show  that  bandtailing  shifts  the 
bandedge  lower  by  as  much  as  lOOmeV.®’  In  GaAs,  this  shift  in  bandgap  corresponds  to 
a  wavelength  change  from  870  nm  to  935  nm.  As  more  carriers  are  injected  into  the  gain 
region,  the  quasi-Fermi  level  continues  to  move  further  into  the  conduction  band.  This 
movement  causes  the  peak  of  the  carrier  concentration  to  reside  at  higher  energies  in  the 
conduction  band  and  lower  energies  in  the  valence  band.  Because  the  spontaneous 
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emission  intensity  depends  on  the  gap  between  these  peak  energies,  as  the  current  is 
increased,  the  peak  of  the  spontaneous  emission  spectra  will  also  increase  in  energy. 

More  carriers  injected  into  the  bands  provide  a  larger  optical  gain.  Eventually  this 
gain  will  overcome  any  losses,  at  which  time  the  device  will  move  into  the  amplification 
regime.  If  a  resonant  cavity  exists,  the  device  will  begin  to  lase  and  the  output  spectrum 
significantly  narrows.  At  this  point,  carriers  are  rapidly  removed  from  the  upper  state  by 
stimulated  emission,  and  the  quasi-Fermi  levels  no  longer  increase  in  energy. 

Electro-luminescence  measurements  reveal  the  behavior  of  the  shifts  in  peak 
wavelength  of  the  emission  versus  current  injection.  The  maximum  gain  available  exists 
at  the  peak  of  the  electro-luminescence  spectrum.  The  behavior  and  location  of  this  peak 
must  be  determined  to  find  the  maximum  nonlinearity. 

Electro  luminescence  experiments  were  conducted  on  a  single  mode,  4-p.m  wide 
waveguide,  with  an  offset  angle  of  0.9°  to  the  primary  cleavage  plane  to  fiiistrate  lasing. 
The  0.9°  inhibits  lasing  except  at  very  high  currents  thereby  permitting  the  study  of  the 
behavior  of  the  device  well  within  the  degeneracy  region.  The  device  was  cleaved  to  an 
overall  length  of  5  mm. 

As  shown  in  Figure  3.4(a),  the  low  injection  case,  130mA,  produces  a  broad  spectrum 
with  a  peak  of  920  nm  consistent  with  the  expected  reduction  of  the  bandgap. 
Application  of  additional  current,  plots  (b)  through  (d),  fills  available  band  states  and 
shifts  the  wavelength  to  higher  energies.  Finally,  the  gain  available  in  the  device 
overcomes  the  angular  slant  of  the  output  facets,  and  the  device  begins  to  lase,  Figure  3.4 
(e).  Figure  3.4(f)  shows  the  output  optical  power  versus  current  for  this  device.  By 
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graphical  extrapolation,  the  lasing  threshold  current  is  deduced  to  be  650mA.  Also 
shown  in  Figure  3.4(f)  is  a  plot  of  the  position  of  the  peak  of  the  electro-luminescence 
spectra  with  injection  current.  In  the  region  below  lasing  threshold,  the  peak  shifts 
linearly  with  the  increase  in  current  then  shifts  very  slightly  in  the  opposite  direction  at 
current  levels  above  the  lasing  threshold.  Remember  that  the  quasi-Fermi  levels  no 
longer  shift  with  increasing  current  injection.  Thus  the  increase  in  the  amplifier’s  peak 
wavelength  after  the  onset  of  lasing  is  most  likely  be  a  result  of  device  heating. 
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Figure  3.4  Electro-Luminescence  Measurements  on  G2-2329  5.5mm  long  4um  wide 
waveguide.  Plots  (a)  -  (d)  show  spontaneous  emission  behavior,  plot  (e)  shows  the  device 
after  the  onset  of  lasing.  Plot  (f)  shows  the  output  power  and  emission  peak  versus 

injection  current. 
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3.2  MEASUREMENT  OF  GAIN  SPECTRUM 


Critical  to  the  successful  operation  of  the  soliton  device  is  the  ability  of  an  input 
optical  beam  to  extract  electrons  from  the  conduction  band  thereby  altering  the  carrier 
population.  As  discussed  in  section  3.1.2,  the  peak  of  the  gain  will  follow  the  electro¬ 
luminescence  peak,  however,  for  a  fixed  current  injection  level,  the  magnitude  of  the  gain 
will  vary  by 

3 

g(v = / 1/-,«  j  , ;]  (3-2) 


where  m*  is  the  reduced  mass  of  the  carriers  in  the  gain  region,  fcr  is  the  oscillator 
strength  of  the  transition  from  conduction  to  valence  band  (approximately  23  eV  for 
GaAs),  Wr  is  the  linear  index  of  refraction  which  is  3.590  for  GaAs,  nio  is  the  electron  rest 
mass,  ^  is  the  energy  where  the  gain  is  calculated,  is  the  unperturbed  bandgap  energy, 
and  finally  /„(^2)  and  are  the  values  of  the  quasi  Fermi  functions  at  weighted 
energies  ^2  for  the  conduction  band  and  ^1  for  the  valence  band.  Equation  3.2  was 
developed  in  Chapter  1.  Note  that  and  the  energies  are  set  such  that 

momentum  is  conserved  by  the  transition.  Using  this  equation  is  cumbersome,  as  the 
quasi  Fermi  functions  require  knowledge  of  the  carrier  concentration  in  the  active  region, 
which  in  turn,  is  directly  proportional  to  the  injection  current.  Assuming  a  quantum 
efficiency  of  unity,  the  carrier  concentration  is  estimated  by  assuming  that  the  injected 
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carriers  recombine  only  in  the  active  region  with  a  characteristic  recombination  time  To, 
typically  0.6  nanoseconds  for  GaAs.  Note  that  the  lifetime  used  here  is  the  combination 
of  the  spontaneous  radiative  and  nonradiative  recombination  lifetimes.  Thus  the  average 

cfix 

carrier,  electron  and  hole,  concentrations  wie  N  =  ,  where  /  represents  the  injected 

Ad 


current  in  amps,  A  is  the  cross  sectional  area  of  the  topside  contact,  and  d  is  the  thickness 
of  the  active  region.  Next  the  quasi-Fermi  levels,  ^  and  are  found  by  fixing  the 
valence  band  energy  to  zero  and  solving 


=-kTH 


JL 


(3.3) 


where  Nc  and  Nv  are  the  effective  density  of  states  values  for  the  conduction  and  valence 
bands  respectively.  In  this  case  Nc=4. 4x10*^ /cm^  and  Nv=8.2xl0**/cm*  were  used.  Now 
the  Fermi  distribution  function  can  be  evaluated  using 


fjv  =  - 
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The  last  step,  the  weighted  energies  are  calculated  by  maintaining  the  conservation  of 
momentum  in  the  transition. 
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(3.5) 


m, 

* 

The  gain  spectrum  was  measured  by  launching  a  beam  through  a  single  mode 
waveguide  amplifier  formed  using  a  4  pm  wide  waveguide  electrically  contacted  at  the 
top  for  earner  injection.  The  beam  was  end  fire  coupled  into  the  waveguide  using  a  40x- 
microscope  objective,  and  the  output  of  the  waveguide  was  focused  onto  a  detector  using 
a  20X-microscope  objective.  To  isolate  the  gain  fi'om  the  spontaneous  emission,  the 
input  was  chopped  and  measured  using  a  detector  connected  to  a  lock-in  amplifier.  The 
experimental  set-up  for  this  experiment  is  shown  in  Figure  3.5. 


Variable  Current 


Figure  3.5  Experimental  Setup  for  Gain  Measurements 


The  absolute  gain  of  this  device  could  not  be  measured  as  there  was  no  way  of 
measuring  coupling  and  cavity  losses,  however,  the  shape  and  spectral  profile  of  the  gain 
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are  valuable  in  fitting  to  the  model.  Results  of  the  measurement  and  a  theoretical  fit  are 
shown  in  Figure  3.6.  The  modeled  fit  assumed  unity  quantum  efficiency,  oscillator 
strength  of  46  eV,  and  a  carrier  lifetime  of  0.6  ns.  As  described  in  chapter  1,  the  valence 
band  quasi-Fermi  level  is  a  strong  function  of  the  hole  injection  fi’om  the  AlGaAs 
cladding  material.  As  such,  the  valence  quasi-Fermi  level  is  modeled  as  a  constant 
energy  of  10.5  meV  above  the  valence  band  edge.  This  value  was  determined  by  best 
curve  fit.  Fitting  parameters  used  here  agree  well  with  established  constants  for  high 
quality  GaAs. 


Figure  3.6  Experimentally  measured  gain  response  spectrum.  Solid  line  represents 

theoretical  curve. 
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3.3  MEASUREMENT  OF  OPTICAL  NONLINEARITY 


The  third  order  optical  nonlinearity  was  measured  by  using  an  interferometer  to 
measure  a  phase  shift  induced  onto  an  optical  beam  in  one  leg  of  the  interferometer.  In 
this  experiment,  a  Mach-Zehnder  interferometer  was  used,  see  Figure  3.7. 

If  the  input  beam  is  passed  through  an  element  with  an  intensity  dependent  index  of 
refraction,  «  =  «„  +  n^I ,  the  beam  will  encounter  an  intensity  dependent  change  in 


Figure  3.7  Mach-Zehnder  Interferometer 


optical  phase,  ^  =  LknJ ,  where  L  is  the  length  of  the  device  and  k  is  the  propagation 
constant.  An  interference  pattern  is  established  by  offsetting  the  two  arms  of  the 
interferometer  slightly  to  create  multiple  peaks  of  spatial  frequency  where  0  is 


64 


the  angle  between  the  two  beams  recombining  on  the  image  plane.  A  one-dimensional 
CCD  array  captured  the  interference  pattern  and  dumped  the  output  to  a  digital 
oscilloscope  with  averaging  capability.  The  averaging  function  was  necessary  to  average 
out  jitter  in  the  interference  pattern  that  is  caused  by  air  movements  and  mechanical 
vibrations  in  the  room.  The  shift  in  the  pattern  was  measured  between  two  input  optical 
powers  of  12.5mW  and  50mW  at  a  variety  of  wavelengths.  The  input  current  to  the 
device  was  pulsed  using  a  1.5  microsecond  on  to  48.5  microseconds  off  ratio  to  maintain 
consistency  with  the  measurements  of  spatial  solitons.  The  optical  beam  was  also 
temporally  shaped  using  an  acousto-optic  modulator  to  track  the  electrical  modulation. 
The  optical  pulses  were  derived  from  a  cw  titanium  sapphire  laser  pumped  by  a  mode- 
locked,  fi’equency  doubled  Nd:YAG  laser.  Each  modulated  optical  pulse  was  one 
microsecond  in  duration  and  timed  to  pass  through  the  amplifier  Vi  microsecond  after 
electrical  current  started  to  inject  into  the  amplifier.  The  delay  allowed  the  amplifier  to 
stabilize  prior  to  the  optical  pulse  arriving.  Figure  3.8  shows  the  measured  output  from 
the  interferometer.  The  maximum  observable  shift  in  the  interference  pattern  was  3.57t. 
Using  the  phase  shift  equation,  A<|>  =  kLnJ^,  this  corresponds  to  a  nonlinear  coefficient 
of  5.7  X  lO"'”  cm^AV  and  occurred  at  the  peak  of  the  gain  spectrum.  This  corresponds 
with  the  normalized  plasma  and  Kramers-Kronig  theoretical  fits  shown  also  in  Figure  3.8. 
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Figure  3.8  Mach-Zehnder  phase  shift  measurements. 

3.4  MEASUREMENT  OF  INDEX  SHIFT  WITH  CARRIER  INJECTION 

Using  the  same  setup  used  in  the  nonlinearity  measurement,  the  phase  shift 
induced  by  carrier  injection  was  also  measured.  In  this  experiment,  a  constant  optical 
power  of  5mW  was  maintained  at  904  nm,  while  changing  the  electrical  injection  into  the 
waveguide.  The  results  are  shown  in  Figure  3.9.  Using  the  linear  relationship  An  =  w^J 
where  ric  is  a  constant  relating  the  change  in  current  concentration  to  the  change  in  the 
index  of  refraction  An,  =  -7.521x10”’ c/n^  /  A  .  This  is  used  to  predict  the  effect  of 
the  electro-optic  steering  mechanism  used  to  steer  solitons  in  devices  considered  later  in 
this  research. 
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Current  Density  (amp/cm^) 


Figure  3.9  Device  phase  change  with  varying  current.  4  um-wide  waveguide,  1mm 
long  measured  at  904nm,  Current  injection  10-40ma. 


3.5  CHARACTERIZATION  OF  THE  OHMIC  CONTACTS 

Although  not  critical  to  the  preliminary  devices  studied  here,  the  quality  of  the 
electrical  ohmic  contact  made  to  the  p-side  of  the  wafer  is  critical  in  the  optimization  of 
the  device  efficiency  and  in  reducing  resistive  heating  due  to  contact  resistivity.  To 
measure  the  resistivity  to  the  contacts,  a  test  pattern,  called  a  Transmission  Line  Pattern 
(TLM),  was  fabricated  onto  the  p-side  of  the  wafer.  The  n-side  contact  was  not  analyzed 
as  it  was  a  large  area  metallization,  in  direct  contact  with  the  heat  sink  and  separated  from 
the  devices  by  several  hundred  microns.  The  p  side  contacts,  on  the  other  hand,  were 
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within  1  micron  of  the  junction  and  were  narrow  contacts  passing  a  large  amount  of 
current. 

Resistance  in  the  semiconductor  can  be  modeled  as  a  series  of  resistances.  With  an 
ohmic  contact,  the  resistances  are  the  resistance  the  contact  forms  with  the  surface,  and 
the  bulk  resistance  of  the  material  between  the  contact  pads.  The  resistance  between  two 
closely  spaced  pads  where  the  distance  between  the  pads  is  much  less  than  the  distance 
laterally  along  the  interface,  is  R  =  where  Rcontact  is  the  resistance 

between  the  contact  and  the  material,  d  is  the  distance  between  the  two  contacts,  and  RbuVc 
is  the  resistivity  of  the  material.  The  contact  resistance  can  be  found  by  measuring  the 
resistance  between  two  pads  at  different  spacings  and  fitting  the  results  to  the  equation. 

The  TLM  pattern  consists  of  a  row  of  50x300  micron  ohmic  pads  separated  by 
increasing  distances  of  2.5, 10, 20,  and  40  microns.  A  diagram  of  this  pattern  is  shown  in 
Figure  3.10.  Note  that  the  distance  between  the  pads  is  much  less  than  the  distance  fi-om 
side  to  side.  Resistivity  measurements  are  made  between  adjacent  pads. 

The  specific  contact  resistivity  can  be  calculated  by  assuming  that  the  current  passing 
between  the  two  pads  crowds  through  an  area  at  the  front  portion  of  the  each  pad.^  ’ 
This  greatly  reduces  the  effective  area  the  current  passes  through,  and  minimizes  the 
increase  in  resistance  created  by  the  current  passing  underneath  the  contact  through  the 
semiconductor.  The  distance  into  the  pad  is  referred  to  as  the  transfer  lengthy  denoted 
here  by  Lt.  The  width  of  the  area  is  the  width  of  the  pad,  w.  The  x  axis  intercept  of  the 
distance  versus  resistance  plot  denotes  twice  the  value  of  this  transfer  length.  As  the 
transfer  length  is  small  compared  to  the  length  of  the  contact,  the  short  surface  resistivity 
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under  the  pad  is  ignored.  Therefore,  the  specific  contact  resistivity  in  (Q  cm^)  of  the  pad 
becomes 

^c=Rccn,actL,^  (3.6) 

Additional  formulas  and  measurement  techniques  exist  to  pinpoint  exact  specific 
contact  resistances^^’^^  however,  the  measurement  and  analysis  preformed  here  are 
sufficient  to  judge  the  relative  health  of  the  p  type  contact. 

A  typical  plot  of  these  results  is  shown  in  Figure  3.11.  The  pad  resistivity  is  V2  of  the 
y  intercept  of  the  plot  seen  here  as  0.4  Q.  The  transfer  length  is  V2  of  the  x-axis  intercept 
or  5  microns.  Therefore,  the  pad  resistivity  is  6.0x1 .  Using  this  equation  3.6, 


a  4-micron  pad,  1-mm  long  will  have  a  total  resistance  of  0. 1 5  Q. 


Pads  Separated  Transfer  Length  Lj 


by  Varying 


Distances 


Figure  3.10  Diagram  of  a  TLM  pattern.  Pad  separation  is  2.5,  5, 10, 20,  and  40 

microns. 
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Resistance  (ohms) 


Figure  3.11  TLM  p  type  pad  resistivity  measurements 
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CHAPTER  4:  SOLITON  MEASUREMENT 


The  theory  and  measurements  outlined  in  Chapters  1  and  3,  indicate  that  a 
nonlinearity  strong  enough  to  produce  viable  solitons  exists  in  GraAs  semiconductor 
optical  amplifiers.  This  chapter  describes  the  device  that  was  fabricated  along  with  the 
measurements  that  demonstrate  and  characterize  these  solitary  waves. 

4.1  DEVICE 

The  active  semiconductor  gain  medium  used  in  this  work  was  a  double 
heterostructure  Alo.2Gao.8As/GaAs  p-i-n  structure  grown  on  a  (100)  oriented  n-doped 
GaAs  substrate,  with  an  intermediate  3000A  thick  n-doped  GaAs  buffer  layer.  The  lower 
n-doped  AlGaAs  layer  was  1.5  pm  thick,  the  intrinsic  GaAs  active  region  was  0.2  pm 
thick,  and  the  p-doped  upper  cladding  layer  was  1.0  pm  thick.  A  lOOOA  p-doped  contact 
layer  covered  the  upper  clad  layer.  The  devices,  shown  in  Figure  4.1,  were  formed  by 
depositing  a  60  pm  wide  Au/Zn/Au  contact  over  the  wafer  and  then  etching  5  pm  wide, 
2000  A  deep  current  isolation  trenches  to  either  side  of  the  gold  contact.  The  sample  was 
then  thinned  down  to  200  pm  and  a  Ni/Ge/Au  n-type  ohmic  contact  was  applied  to  the  , 
backside.  After  annealing  at  420“C  for  120  sec  to  form  the  n  and  p  ohmic  contacts,  the 
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device  was  cleaved  to  650  jim  and  mounted  on  a  copper  stud.  To  frustrate  lasing  in  the 
active  device,  the  contact  stripes  were  angled  at  4°  from  the  primary  cleavage  plain  and 
then  aligned  to  the  corresponding  15°  launch  angle.  Anti-reflection  coatings  were  not 
applied.  The  size  of  the  semiconductor  amplifier’s  active  area,  the  lack  of  facet 
antireflection/protective  coatings,  and  the  amount  of  current  required  to  drive  the  device 


Cross-Section  View 

60iiTn  widp,  Au/Zn/An  contact 
GaAs  capping  layer 
p  AlGaAs  cladding  layer 
Undoped  GaAs  active  layer 
n  AlGaAs  claddinp  laver 

n+  GaAs  substrate 
Au/Ge  contact 

Figure  4.1  Diagram  of  soliton  device 


72 


into  inversion  mandated  the  use  of  a  pulsed  current  driver.  The  current  pulses  were 
1.5  ps  in  duration,  had  a  repetition  frequency  of  20  kHz  (50  ps  cycle)  and  had  a  peak 
amplitude  of  750  mA,  unless  otherwise  stated  in  the  experiment.  The  device  was  found 
to  lase  at  an  instantaneous  current  of  1 .2  Amps. 

4.2  EXPERIMENTAL  SET-UP 

The  experimental  setup  that  was  used  to  measure  spatial  solitons  is  shown  in  Figure 
4.2.  This  experiment  used  the  same  tunable  Lexel  Model  480-D-lX,  Ti-sapphire  light 
source  pumped  by  a  Coherent  Antares  Model  76-S,  frequency  doubled,  mode  locked, 
Nd:YAG  laser  as  the  Mach-Zehnder  experiment  described  in  Chapter  3.  As  with  the 
previous  experiment,  the  beam  was  temporally  shaped  to  1  ps  in  duration  using  a 
Quantum  Technology  Model  304A  acousto-optic  modulator  to  match  the  1.5  ps  current 
pulses  applied  to  drive  the  amplifier.  This  modulator  provided  300-mW  power  at  the 
peak  of  each  pulse.  Following  a  polarization  filter  and  rotator,  the  beam  diameter  was 
then  spatially  reduced  using  a  2X-reducing  telescope,  so  that  it  would  fit  within  the 
launch  objective  lens.  No  attempt  was  made  to  adjust  the  spatial  profile  of  the  beam  to 
match  the  input  astigmatism  of  the  waveguide.  The  beam  was  end-fire  coupled  using  a 
40X  long  working  distance  fiber  launch  objective,  and  the  output  facet  of  the  amplifier 
was  imaged  onto  a  Watek  Model  09  calibrated  CCD  camera  using  a  20X  microscope 
objective.  The  output  beam  profiles  were  measured  by  frame  grabbing  a  unity-gamma 
calibrated  camera  image  of  the  amplifier  output  with  a  Snappy  frame  grabber  and  the 
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beam  profiles  were  analyzed  using  Image  Tool,  a  software  downloaded  fi’om  the 
University  of  Texas  Health  Science  Center. 


Accousto-Optic 

Modulator 

Half-Wave  Plate 

Polarization 
Beam  Splitter 
Half-Wave  Plate 


Figure  4.2  Soliton  Measurement  Experimental  Setup 


4.3  COUPLING  AND  PROPAGATION  LOSSES 

Coupling  of  a  laser  beam  into  a  semiconductor  slab  waveguide  using  a  single 
objective  lens  is  invariably  an  inefficient  process.  Firstly,  the  uncoated  facet  of  the  high 


refractive  index  semiconductor  leads  to  large  Fresnel’s  reflection  losses.  Secondly, 
focusing  the  beam  from  the  Ti:  sapphire  laser  to  a  spot  size  equivalent  to  the  thickness  of 
the  active  waveguide  core  (0.2nm)  is  impracticable.  Furthermore,  the  asymmetric 
numerical  aperture  of  a  slab  waveguide  is  another  tough  obstacle  to  overcome  and  indeed 
the  optimum  launch  condition  for  a  slab  waveguide  does  not  correspond  to  the  ideal 
condition  for  launching  a  soliton  waveguide.  However,  in  order  to  measure  the  coupling 
efficiency,  the  waveguide  propagation  losses  needs  to  be  known.  There  are 
insurmountable  obstacles  that  prevented  us  from  determining  the  actual  propagation 
losses.  In  particular,  the  propagation  losses  need  to  measured  at  either  some  specific 
injection  current  density  or  at  a  wavelength  where  the  corresponding  intrinsic  band-to- 
band  absorption  losses  are  less  severe.  Although  in  a  single  mode  waveguide,  the 
propagation  losses  could  have  been  measured  by  thermal  scanning  of  the  Fabry-Perot 
transmission  resonances  of  the  waveguide  resonator  cavity,  these  losses  would  not  be 
meaningful  for  the  slab  waveguide.  Consequently,  the  coupling  losses  due  to  mode 
mismatch  were  estimated  by  calculating  the  overlap  integral  of  the  waveguide  mode  with 
the  focused  laser  beam  spot  size.  The  coupling  losses  due  to  Fresnel’s  reflection  at  the 
dielectric  interface  was  also  evaluated  theolretically.  Reflectivity  Losses 

Light  was  coupled  into  the  soliton  devices  by  end-fire  coupling.  This  type  of 
coupling  suffers  from  losses  induced  by  the  reflectivity  created  at  the  interface  of  air  and 
the  semiconductor  fecet.  The  reflectivity  losses  can  be  calculated  by  using  the  linear 
index  of  refraction  of  the  GaAs  and  the  angle  of  coupling  in  Fresnel’s  equations, 
equations  (4.1)  and  (4.2).  Although  the  input  beam  propagates  in  the  TE  mode  inside  the 
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slab  waveguide,  the  polarization  is  TM  to  the  physical  interface,  which  is  positioned  at 
15°.  The  power  transmittance  at  this  interface  is  found  using  the  Fresnel’s  equation  for 
dielectric  interfaces. 


|2 


and  the  transmission  coefficient,  x,  is  found  using  the  equation 

^GciAs  cos(e^J  +  n^,cos(QcaAs) 


(4.2) 


using  nGa4s=3.59,  and  the  incident  and  transmitted  angles  of  15°  and  4°  respectively, 
the  power  transmittance  is  T=69.4%.  Note  that  the  Brewster’s  angle  would  be  74.4°. 

In  addition  to  reflection  losses,  the  mismatch  between  the  mode  created  by  the 
launching  objective  and  the  waveguide  mode  reduces  the  energy  coupled  into  the 
waveguide.  The  loss  is  particularly  bad  in  the  vertical  direction  where  the  size  of  the 
waveguide  is  only  0.2  microns,  while  the  size  of  the  beam  at  the  facet  is  6.0  microns  full 
width  half  maximum.  The  loss  in  the  horizontal  direction  is  less  critical  as  the  waveguide 
is  a  slab  and  is  able  to  accept  any  width  of  input  beam.  The  coupling  efficiency  between 
the  input  mode  and  the  waveguide  mode,  F,  is  expressed  as  the  integrated  overlap  of  the 
guided  mode  and  the  launched  beam  mode 
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Uni,n  are  the  normalized  modal  shape  functions  (x  and  y)  of  the  waveguide  and  E(x,y) 
is  the  total  field  of  the  input  beam.  Assuming  that  the  input  beam  is  Gaussian  with  a 
FWHM  of  6  microns,  the  coupling  efficiency  to  a  slab  waveguide  mode  in  a  0.2-micron 
waveguide  is  9.18%.  The  coupling  between  a  6  micron  FWHM  Gaussian  input  beam  and 
a  6-micron  FWHM  sech^  beam  is  93%.  The  resulting  total  coupling  becomes 
r  =  0.93*0.0918* 0.694  =  0.059 or  5.925%.  Other  effects  such  as  aberrations,  dust,  and 
noise  will  decrease  this  number.  As  such,  5.0%  coupling  efficiency  is  used  throughout 
the  rest  of  these  calculations. 

4.4  BEAM  PROFILES  THROUGHOUT  SYSTEM 

In  order  to  characterize  the  magnitude  of  a  soliton  using  the  soliton  stability  equation 
(4.4),  the  magnitude  and  size  of  the  launched  beam  must  be  determined. 

Pa=-^  (4.4) 

n^n^k 

P  is  the  power  in  the  soliton  beam,  a  is  the  waist  of  the  soliton,  w  is  the  waist  of  the 
confined  mode  in  the  slab  waveguide,  no  and  are  the  linear  and  third  order  nonlinear 
indices  of  refi'action,  and  k  is  the  propagation  constant. 

Beam  powers  were  measured  using  a  Newport  Model  835  calibrated  power  meter, 
and  the  beam  size  was  measured  using  a  Watek  Model  09  CCD  camera  to  image  the 
beam,  and  a  Sensor  Physics  Laser  Test  LS-4W  beam  profiler  to  analyze  the  captured 
images.  Calibration  measurements  confirmed  that  the  Watek  camera’s  gamma  correction 
factor  was  0.45,  and  the  gamma  correction  was  accounted  for  in  the  program.  Beam 


77 


waist  positions  and  size  were  determined  by  measuring  beam  shape  at  several  points  and 
fitting  the  waist  to  the  Gaussian  propagation  equation 


\  X'o 


+1 


(4.5) 


J 


where  wo  is  the  beam  waist,  z  is  the  distance  from  the  waist,  zo  is  the  position  of  the 
waist,  and  X  is  the  wavelength  of  the  radiation.  The  familiar  Rayleigh  range  is  given  in 


equation  (4.5)  by  .  The  beam  waist  equation  (4.5)  describes  the  propagation  of 

X 


a  perfect  Gaussian  beam;  deviations  in  the  fit  to  this  equation  are  expected  as  the  beam 
emitted  from  the  Ti:  sapphire  laser  is  close  to,  but  not  perfectly  Gaussian.  Of  prime 
importance  is  the  waist  of  the  beam  at  the  input  to  the  device.  Figure  4.3  shows 
measurements  taken  to  confirm  a  beam  waist  of  5.8  microns.  The  Gaussian  profile  of  the 
beams  measured  by  calibrated  camera  is  shown  in  Figure  4.4.  These  beam  profiles 
matched  a  1/e^  Gaussian  fit  by  96%  and  94%  for  the  x  and  y  axis  respectfully. 
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4.5  SOLITARY  WAVE  FORMATION 


The  solution  of  the  nonlinear  Schroedinger  wave  equation  (NSWE)  states  that  the 
spatial  transverse  envelope  of  the  soliton  is  sechi^.  The  input  beam,  however,  is  Gaussian. 
Thus,  the  beam  must  transform  from  a  Gaussian  profile  to  the  sech^  profile  along  the 
length  of  the  device.  The  nonlinearity  that  leads  to  the  formation  of  the  soliton  within  the 
device  assists  in  the  transformation,  as  one  of  the  properties  of  solitons  is  their  ability  to 
form  in  the  presence  of  noise.  Energy  not  conforming  to  the  soliton  mode  is  not 
contained  in  the  solitary  wave  and  diffracts  away  from  the  soliton.  Conducting  soliton 
formation  measurements  along  the  device  would  not  be  feasible,  however,  using  a 
simulation  (beam  propagation  code,  appendix  C),  the  solitary  wave  formation  can  be 
modeled.  Figure  4.5  shows  a  graphical  representation  of  the  evolution  of  a  Gaussian 
beam  with  the  same  waist  as  that  launched  in  the  experiments  as  it  propagated  through 
the  sample. 


Lateral  Dimension  (urn) 


Figure  4.5  BPM  Simulation  of  a  Gaussian  beam  launched  into  a  nonlinear  optical  slab 

waveguide. 
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In  the  experiment,  904-nm  wavelength  light  was  launched  into  the  device  with  a 
focused  beam  waist  of  six  microns.  The  measured  input  intensity  was  1.1  MW/cm^  in 
this  experiment  and  the  propagation  distance  was  650  microns.  18mW  of  instantaneous 
output  power  was  measured  coming  out  of  the  device  after  subtracting  the  spontaneous 
emission  the  device  produced  without  any  optical  input.  Figure  4.6  shows  an  actual 
measured  output  of  the  device  along  with  the  theoretical  prediction  made  by  the  BPM. 
Using  the  BPM  code  to  iteratively  fit  the  peak,  the  nonlinear  112  coefficient  was 
n2=8.83xl0'  cm  AV.  From  the  measured  soliton  beam  waist,  of  6  microns,  and  the 
input  optical  intensity  of  l.lMW/cm^  and  taking  into  account  the  coupling  losses 
described  above,  the  nonlinear  coefficient  n2,  was  calculated  to  be  5.6x10’*°  cm^AV. 


Microns 

Figure  4.6  Output  beam  profile  of  a  slab-guided  soliton.  Solid  line  is  theoretical  fit  to 
the  expected  beam  profile  generated  using  a  beam  propagation  code.  The  wide  function 
is  the  beam  output  without  the  nonlinearity. 
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4.5.1  NONLINEARITY  VS.  WAVELENGTH 


The  soliton  waist  was  measured  at  a  variety  of  wavelengths  while  keeping  the  device 
electrical  current  and  input  optical  intensity  constant.  These  optical  and  electrical 
parameters  were  carefully  selected  such  that  the  solitons  that  were  formed  contained  the 
largest  percentage  of  the  launched  energy  at  the  lowest  input  optical  intensity.  This  was 
necessary  in  order  to  produce  a  solid  soliton,  but  not  to  saturate  the  nonlinearity. 
According  to  the  theory  developed  in  Chapter  1,  the  nonlinearity  should  be  directly 
proportional  to  the  gain  measured  in  the  material.  This  hypothesis  is  confirmed  by  the 
plot  of  Figure  4.7,  which  shows  the  normalized  nonlinearity  as  a  function  of  the  input 
wavelength.  A  normalized  shape  of  the  measured  gain  curve  is  plotted  on  the  same  axes 
for  comparison. 


Figure  4.7  Nonlinearity  versus  wavelength.  Nonlinearity  found  by  taking  the  inverse 
of  the  product  of  the  beam  power  and  the  beam  width.  Normalized  gain  curve  displayed 

as  the  solid  line. 
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4.5.2  NONLINEARITY  VS.  CURRENT 


The  same  measurement  as  the  previous  section  was  performed,  but  in  this  experiment 
the  current  was  varied  while  keeping  the  intensity  and  wavelength  fixed.  The  theoretical 
calculation  in  this  situation  as  shown  previously  predicts  that  the  nonlinearity  should 
follow  the  gain  curve.  The  peak  of  the  gain  curve,  however,  shifts  with  changing  current. 
These  effects,  taken  together,  may  result  in  the  nonlinearity  remaining  relatively  constant. 
This  effect  was  seen  in  beam  waist  measurements  taken  against  increased  current,  shown 
in  Figure  4.8.  The  inset  of  Figure  4.8  shows  the  dependence  of  the  output  power  and 
beam  waist  on  the  injection  current.  At  very  high  currents,  the  beam  waist  remained 
constant,  while  the  output  power  continued  to  increase.  At  this  point,  side  lobes  began  to 
form  around  the  main  lobe  possibly  signaling  the  onset  of  higher  order  self-guided 
modes.^"*  The  falloff  of  the  nonlinearity  results  fi'om  energy  transferring  into  the  side 
lobes  of  the  output  beam. 
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Figure  4.8  Nonlinearity  versus  Current  Density.  Inset  shows  the  relationship  between 
beam  waist  and  output  power.  At  high  currents,  the  beam  waist  no  longer  shrinks, 
however,  side  lobes  begin  to  form  around  the  main  lobe. 


4.5.3  NONLINEARITY  VS  INTENSITY 

In  the  next  series  of  measurements,  the  best  soliton  was  formed  using  the  lowest  input 
optical  intensity.  The  input  optical  intensity  was  then  varied  and  the  soliton  waist  was 
measured.  Figure  4.9  depicts  the  measured  soliton  beam  waists  versus  input  intensity.  In 
theory,  the  gain,  and  thus  the  nonlinearity  must  eventually  saturate. 
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optical  Input  Power  (mW) 


Figure  4.9  Nonlinearity  versus  input  optical  power.  Inset  shows  how  beam  waist  and 
power  change  with  input  optical  power. 

This  measurement,  although  producing  a  small  roll-ofF  in  the  nonlinearity,  may  not 
give  a  precise  measure  of  the  saturation  intensity.  As  the  intensity  was  increased,  a 
pedestal  and  side  lobes  began  to  form  around  the  soliton.  This  pedestal  may  be  excess 
energy  not  able  to  be  contained  within  the  soliton,  or  may  be  the  onset  of  a  higher  order 
mode  as  suggested  in  other  soliton  experiments.^^  Although  a  second  order  soliton  in  a 
pure  Kerr  like  material  would  not  form  3  distinct  peaks,  the  nonlinearity  in  this  material 
is  not  entirely  Kerr  like. 
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4.5.4  NONLINEARITY  VS.  WAVELENGTH:  PEAK  INTENSITY 


As  a  final  confirmation  of  the  nature  of  the  nonlinearity,  the  peak  of  the  soliton  was 
measured  versus  wavelength.  The  power  contained  in  the  beam  is  proportional  to  the 
peak  of  the  function,  thus,  the  behavior  of  the  peak  provides  insight  into  the  spectral 
nature  of  the  nonlinearity.  Figure  4. 10  displays  the  measured  peak  intensity  of  the  soliton 
against  the  wavelength.  As  expected,  the  nonlinearity  tracks  the  gain  curve.  The 
measured  and  normalized  gain  curve  is  overlaid  for  reference.  As  expected,  the 
nonlinearity  tracks  the  gain  curve.  As  noted  in  the  experiments  above,  the  final  width  of 
the  output  will  only  reduce  to  the  width  of  the  diffracted  mode,  or  approximately  34 
microns  FWHM.  This  produces  the  non-zero  intensities  to  the  lower  energy  side  of  the 
gain  peak.  To  the  high-energy  side,  no  signal  appears  due  to  absorption  in  the 
semiconductor. 


Figure  4.10  Measurement  of  absolution  signal  peak  with  varying  wavelength. 
Nonlinearity  continues  to  follow  the  rough  gain  outline. 
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CHAPTERS:  GUIDING  OF  OPTICAL  BEAM  BY  SOLITON 


WAVEGUIDE 


The  solitons  studied  in  the  previous  chapter  have  demonstrated  the  ability  to  form 
waveguides  capable  of  guiding  other  beams  within  the  index  of  refraction  wake  of  the 
soliton.  The  ability  to  guide  unrelated  beams  leads  to  the  formation  of  devices  capable  of 
switching  information  from  one  input  to  a  variety  of  outputs  in  a  single  stage  if  the 
soliton  can  be  steered  electro-optically. 

5.1  CROSS  POLARIZED  SIGNAL  BEAM 


If  coherent  beams  are  employed  then  the  addition  of  the  second  beam  coincident  upon 
the  soliton  beam  requires  cross  phase  modulation  to  be  taken  into  account.  The  coupled 
equations  governing  the  interaction  are 


dz 


i 

2k 


i  d^E 


2k  d^x^  2 


^E„+ik 


E„+AJEJ 


0 


(5.1) 
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In  these  coupled  equations,  Ea  represents  the  TE  polarized  field  creating  the  soliton, 
E„  represents  the  weak  TM  signal  beam  guided  by  the  first  beam,  k  represents  the 
propagation  wave  vector  of  the  two  beams.  The  /12  terms  are  the  nonlinear  coefficients 
for  the  two  polarization  states  (self  phase  modulation),  and  the  A  terms  are  the  ratio  of  the 
cross  phase  to  self  phase  modulations.  As  the  cross  polarized  signal  beam  is  well  beyond 
the  range  of  the  gain  creating  the  nonlinearity,  the  n^ic  and  g„  terms  are  expected  to  be 
very  small,  and  the  term  large.  Additionally,  the  signal  beam  is  weak  compared  to  the 
soliton  beam,  thus  the  signal  beam  is  not  expected  to  contribute  significantly  to  the  phase 
modulation  of  the  soliton  beam.  Finally,  because  the  index  profile  created  by  the  soliton 
results  fi'om  the  depletion  of  unpolarized  carriers  in  the  conduction  band,  the  change  in 
the  index  of  refraction  for  the  cross  polarization  state  is  anticipated  to  be  nearly  equal  to 
that  created  by  the  soliton.  The  index  will  be  ratioed  by  a  factor,  B,  due  to  the  dispersive 
nature  of  the  index  shift  with  wavelength  as  discussed  in  section  5.2.  Thus,  as  a  first 
order  hypothesis,  the  coupled  mode  equations  reduce  to 


dE,  i 

dz  2k^ 


+  +^^oW2.cKi 


E„ 


dz 


i  d^E, 
2k ^  d^x^ 


+  ikoBn^„\E„\^  E„ 


(5.2) 
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5.2  REFRACTIVE  INDEX  PROFILE  (K-K  RELATIONSHIP)  WITHIN 


SOLITON  TRACK 

The  relationship  between  the  localized  index  of  refraction  change  and  the  solitary 
wave  was  developed  in  Chapter  1  and  confirmed  by  measurements  in  Chapters  3  and  4. 
This  index  of  refraction,  however,  was  the  self-action  index,  or  the  index  that  the  photons 
creating  the  soliton  see.  Photons  outside  this  spectral  range  encounter  a  different  index  of 
refraction.  As  shown  in  Chapter  1,  the  mechanisms  creating  the  change  in  the  index  fall 
into  two  categories;  contributions  from  a  carrier  plasma  effect,  and  contributions  from  the 
shift  in  the  absorption  edge  due  to  the  gain  depletion.  In  the  first  effect,  the  index 
contribution  from  the  plasma  effect  is  proportional  to  the  wavelength  squared,  as  long  as 
the  wavelength  remains  far  from  the  resonance,  as  is  the  case  at  960nm.  The  spectral 
dependence  of  the  change  in  reflective  index  resulting  from  the  shift  in  the  absorption, 
however,  is  related  by  the  Kramers-Kronig  relationship,  developed  in  appendix  B,  to  the 
spectrum  of  the  absorption  change.  The  spectral  relationship  between  the  band  edge  shift 
and  refractive  index  is  given  by 

Aw(m)  =  -  f -  — -  (5.3) 

L  ft)  -  o' 


where  Aa((D)  is  the  change  in  the  absorption  spectrum.  This  change  is  modeled  by 
changing  the  carrier  concentration  in  the  semiconductor  gain  equation,  and  taking  the 
difference  between  the  gain  curves.  The  Kramers-Krbnig  transform  is  applied  to  the 
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resulting  gain  difference  yielding  a  spectral  dependence  of  the  change  in  the  index  of 
refraction.  Combining  these  two  effects  yields  the  plot  of  the  index  of  refraction  curve 
shown  in  Figure  5.1.  Using  this  spectrum  and  the  fact  that  the  maximum  calculated 
change  in  the  index  across  the  soliton  is  5x1 0"^,  the  maximum  expected  change  in  the 
index  of  refraction  experienced  at  960nm  would  be  2.49x10“*.  The  sech^  shape  of  the 
soliton  defines  the  refractive  index  profile  of  the  waveguide. 


400  600  800  1000  1200  1400  1600  1800 

Wavelength 


Figure  5.1  Change  in  refractive  index  due  to  carrier  induced  nonlinearities. 
Nonlinearity  caused  by  a  depletion  of  carriers  at  904nm  by  an  injected  pump  beam. 


5.3  GRADED  INDEX  WAVEGUIDING  (EXPECTED  BEAM  SHAPE) 


Waveguide  modes  can  be  calculated  analytically  for  only  a  very  few  specific  set  of 
waveguide  geometries.  A  development  for  square  profile  waveguides  is  shown  in 
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appendix  C.  In  this  method,  trial  functions  are  established  for  each  of  the  zones  of  the 
waveguide,  and  boundary  conditions  are  matched  to  provide  a  solution.  This  method  can 
be  extended  to  multiple  square  index  systems  using  a  computer  to  approximately  solve 
the  complete  set  of  simultaneous  equations.  Using  this  method,  the  hyperbolic  secant 
squared  shaped  index  profile  can  be  approximated  by  a  series  of  small  square  index 
waveguides,  shown  in  Figure  5.2,  and  boundary  conditions  set  to  solve  for  the  mode 
profile  of  this  system,  overlay  on  Figure  5.3. 


Figure  5.2  Square  approximation  used  to  solve  mode  profiles  for  a  secl^  shaped 
waveguide.  Solid  line  is  the  calculated  field  profile  for  this  waveguide. 


As  a  second  confirmation  of  the  mode  profile  of  this  waveguide,  the  Beam 
Propagation  Method,  appendix  C,  is  used  to  simulate  the  index  profile  of  the  solitary 
wave  waveguide.  In  this  simulation,  the  hyperbolic  secant  squared  shape  is  assumed 
fixed  with  an  end-fire  coupled  signal  wave  guided  down  the  length  of  the  waveguide. 
This  simulation  and  the  resulting  output  beam  shape  are  shown  in  Figure  5.3. 


91 


L 


Figure  5.3  Calculation  and  Simulation  of  Intensity  Profile  in  a  6  micron  FWHM  sech^ 
graded  index  waveguide.  BPM  assumed  a  6-micron  FWHM  Gaussian  launch  beam. 

The  peak  magnitude  of  the  guiding  region  was  set  at  5x10"^  above  the  base  index  of 
refraction.  In  both  of  these  simulations,  it  is  seen  that  the  waveguide  guides  the  beam, 
and  the  final  confinement  beam  has  a  FWHM  of  6  microns  using  a  FWHM  on  the  soliton 
beam  of  6  microns. 

The  analysis  shown  above  is  accurate  for  a  perfectly  Gaussian,  single  mode  input 
beam.  The  device  used  to  create  the  signal  beam,  however,  did  not  emit  a  single  mode 
output.  The  laser  used  was  an  InGaAs  wide  stripe  laser  diode  mounted  and  bonded  in  our 
laboratory.  The  beam  profile  was  measured  using  a  Watek  09  calibrated  camera 
(r=0.45)  and  analyzed  by  a  Spectra  Physics  LM-4W  software  package.  Figure  5.4 
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Figure  5.4  960nm  input  beam  before  2x  reduction  telescope. 


The  waist  size  of  this  beam  was  set  to  match  the  waist  position  of  the  Ti:  sapphire 
laser  beam  that  was  used  to  generate  the  soliton  wave.  The  success  of  the  match, 
however,  is  suspect  as  small  adjustments  that  were  required  to  align  the  launch  optics, 
lens  material  dispersion,  and  material  dispersion  through  the  15°  launch  angle  make  an 
exact  match  of  this  direct  launch  system  impossible.  Using  this  beam  profile  as  the  input 
to  the  soliton  waveguide,  a  BPM  analysis  of  the  output  of  a  soliton  waveguide  device  is 
shown  in  Figure  5.5.  As  can  be  seen  in  Figure  5.4,  the  input  beam  irregularities  provide  a 
large  array  of  possible  launch  conditions.  The  simulation  results  of  Figure  5.5  were 
iteratively  matched  to  the  actual  output  by  varying  the  input  launch  position  and  angle. 
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Figure  5.5  BPM  simulation  of  the  960nm  diode  beam  passing  through  the 
semiconductor  optical  amplifier  with  and  without  the  guiding  soliton  beam. 


5.4  BROAD  AREA  SOLITON  DEVICE 

The  device  characterized  in  Chapter  4  was  used  to  demonstrate  the  ability  of  a  soliton 
to  guide  a  second  unrelated  beam.  Remember  that  this  device  consisted  of  a  60-micron 
wide  active  region  across  a  650-micron  long  device.  The  stripe  was  angled  at  4°  to  the 
cleaved  facets  to  frustrate  lasing.  Due  to  the  high  refractive  index  (3.5)  of  GaAs,  the  non- 
orthogonal  orientation  of  the  stripe  made  it  necessary  to  maintain  a  15®  launch  angle 
external  to  the  cavity.  This  launch  angle  contributed  to  difficulties  in  confining  the  signal 
beam  within  the  soliton  beam  due  to  the  fact  that  the  index  of  refraction  dispersion  of  the 


94 


semiconductor  made  the  launch  angles  different.  The  use  of  a  single  mode  injection 
waveguide  is  anticipated  to  alleviate  this  problem,  however,  a  thorough  investigation  of 
this  type  of  device  shall  be  the  subject  of  additional  research. 

5.4.1  EXPERIMENTAL  SETUP 

An  experiment  to  demonstrate  the  ability  of  the  solitary  wave  waveguide  to  carry 
signal  information  was  conducted.  The  experiential  setup  that  was  used  to  measure 
spatial  soliton  beam  guiding  is  shown  in  Figure  5.6. 

This  experiment  used  the  same  setup  as  the  soliton  measurement  scheme  used  in 
Chapter  4  with  the  inclusion  of  an  orthogonally  polarized  960nm  signal  that  is  added  at 
the  polarization  filter/beam  splitter.  A  tunable  Ti-sapphire  laser  served  as  the  soliton 
creating  light  source.  The  Ti-sapphire  laser  was  optically  pumped  using  a  frequency 
doubled,  mode  locked  Nd:  YAG  laser.  As  with  the  previous  experiment  used  to  measure 
the  soliton  characteristics  shown  in  Chapter  4,  the  beam  was  temporally  shaped  to  1  ps  in 
duration  using  an  acousto-optic  modulator  to  match  the  1.5  ps  current  pulses  applied  to 
drive  the  amplifier.  This  arrangement  provided  300-mW  power  at  the  peak  of  each  pulse. 
The  beam  diameter  was  then  spatially  reduced  using  a  2X-reducing  telescope,  so  that  it 
would  fit  within  the  clear  aperture  of  the  launch  objective  lens. 
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This  beam  was  end-fire  coupled  into  the  wide  stripe  semiconductor  optical  amplifier 
using  a  40X  long  working  distance  fiber  launch  objective,  and  the  output  facet  of  the 
amplifier  was  imaged  onto  a  calibrated  CCD  camera  using  a  20X  microscope  objective. 
To  demonstrate  self-guiding,  a  second  laser  beam,  62  nm  longer  in  wavelength  than  the 
soliton  writing  beam,  was  added  to  the  first  beam  using  a  polarization  beam  splitter.  This 
second  laser  was  a  partially  collimated,  InGaAs  multimode,  wide-stripe  laser,  operating 
at  960nm.  At  the  output,  the  960nm  guided  beam  was  imaged  independently  fi’om  the 
904nm  soliton  forming  beam  using  a  narrow  bandpass  interference  filter.  The  beam 
output  profiles  were  measured  by  fi'ame  grabbing  a  calibrated  camera  image  of  the 
semiconductor  amplifier  output  and  plotting  the  intensity  distribution  along  the  active 
region  of  the  amplifier. 

5.4.2  BEAM  GUIDING  MEASUREMENTS 

Figure  5.7  displays  beam  profile  plots  of  the  best  confined  output  of  this  device 
revealing  the  region  where  the  soliton  waveguide  contained  the  960nm  signal  beam. 
Nearly  3dB  increase  in  signal  strength  is  obtained  in  the  region  corresponding  to  the  peak 
of  the  soliton  profile.  The  input  signal  beam  was  highly  multimode  contributing  to  the 
lack  of  confinement  in  the  output  of  the  device.  This  type  of  poor  collection  is  predicted 
in  the  BPM  simulation  shown  in  Figure  5.5.  In  Figure  5.7,  the  dotted  line  represents  the 
theoretical  shape  of  the  output  mode  of  the  signal  beam  if  a  perfect  launch  beam  was 
used.  It  is  observed  that  the  shape  of  the  theoretical  curve  matches  the  section  of  the 
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output  corresponding  to  the  area  guided  by  the  soliton  beam.  The  guiding  soliton 
generated  in  this  experiment  had  a  FWHM  waist  of  8  microns. 


Figure  5.7  Soliton  Guiding  showing  theoretical  shape  of  the  guided  beam. 


5.4.3  BEAM  TRAPPING  MEASUREMENTS 

For  optimum  beam  guiding,  the  soliton  beam  and  signal  beam  must  be  spatially 
overlapped  at  the  input  to  the  device.  An  analog  to  this  would  be  the  alignment  of  the 
signal  beam  into  a  single  mode  waveguide.  As  with  a  fixed  waveguide,  the  soliton  will 
guide  the  portions  of  the  input  beam  that  leak  into  the  soliton  waveguide.  This  can  be 
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viewed  as  the  overlap  of  that  portion  of  the  signal  beam’s  wave  vectors  with  the 
propagation  modes  of  the  soliton  waveguide.  This  is  also  referred  to  as  beam  trapping. 

To  measure  the  angular  acceptance  of  the  soliton  trapping,  the  experimental  setup 
shown  in  Figure  5.6  was  used.  To  maintain  coupling  at  various  angles,  the  soliton  and 
signal  beams  should  overlap  exactly  at  the  input  to  the  device.  Offsetting  the  input  angle 
of  the  signal  beam  from  a  point  2  meters  from  the  launch  objective  approximated  a 
collinear  input  to  the  launch  objective.  The  percentage  of  power  guided  by  the  soliton 
beam  was  calculated  using  a  beam  propagation  code  with  the  actual  measured  signal 
beam  at  incremental  angles  as  the  input. 

The  measured  amounts  of  energy  at  various  input  angles  contained  in  the  guided 
signal  beam,  along  with  a  theoretical  simulation,  are  graphed  in  Figure  5.8.  Deviations 
from  the  shape  of  the  curve  are  the  result  of  the  misalignment  between  soliton  and  signal 
when  the  signal  beam  was  deviated,  and  also  from  the  aberrations  caused  by  the  deviation 
of  the  signal  beam. 
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Figure  5.8  Ratio  of  signal  diode  power  trapped  by  the  soliton  waveguide  at  the  output 
of  the  semiconductor  optical  amplifier. 

5.5  SINGLE  MODE  LAUNCH  DEVICES 

A  primary  shortcoming  in  the  evaluation  of  the  solitons  was  that  the  launch  angles 
and  modes  could  not  be  precisely  controlled  and  aimed.  To  overcome  this,  a  soliton 
propagation  device  was  designed  with  a  single  mode  waveguide  as  the  input  to  the  soliton 
guiding  region.  This  single  mode  waveguide  provided  gain  as  well  as  confinement  to  the 
launched  beam.  To  determine  the  ability  of  the  soliton  to  be  deflected,  a  triangular  region 
was  inserted  between  the  end  of  the  single  mode  waveguide  and  wide  area,  soliton 
propagation  region.  As  the  index  of  refi’action  could  be  altered  by  carrier  injection,  the 
triangular  region  could  be  used  to  steer  waves  across  it.  Finally,  the  input  to  previous 
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devices  required  arranging  the  stripe  ofF-normal  to  the  facet  to  prevent  the  device  from 
lasing.  Although  the  design  of  this  device  permits  the  waveguide  to  be  normal  to  the 
cleaved  facet,  a  slight  angle  of  0,5°  (1.75°  external  launch  angle)  was  built  into  the 
waveguide  to  ensure  lasing  would  not  occur  if  the  prism  angle  was  set  such  that  a  lasing 
cavity  could  be  created  across  the  device. 

5.5.1  DEVICE  LAYOUT 

Following  the  successful  beam  guiding  experiments,  a  device  designed  to  launch, 
steer,  and  guide  a  signal  beam  from  a  single  input  to  any  number  of  outputs  was 
fabricated.  The  basic  device,  shown  in  Figure  5.9,  incorporates  a  single  mode 
waveguide,  a  steering  prism,  and  finally  a  wide  area,  soliton  propagation  region.  The 
input  single  mode  waveguide  acts  as  an  amplifier  to  the  soliton  generation  beam  and  as  a 
low  loss  waveguide  to  the  signal  beam,  assuming  that  the  wavelength  of  the  signal  beam 
is  below  the  absorption  band  of  the  semiconductor.  Thus,  the  waveguide  is  a  fixed  ridge 
waveguide  with  an  active  contact  to  allow  current  injection  into  the  ridge.  This  first 
waveguide  also  restricts  the  soliton  generating  beam,  as  well  as  the  signal  beam,  to  a  very 
specific  launch  point  and  direction. 
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Steered  Output 


Figure  5.9  Prototype  Soliton  Steering  Device 


An  electro-optic  steering  prism  immediately  follows  the  single  mode  launch 
waveguide.  Creating  a  ridge  weighted  and  electrically  active  triangular  region  following 
the  single  mode  device  forms  this  steering  region.  Deviation  of  the  beam  occurs  when 
the  beam  encounters  a  spatial  linear  phase  chirp  created  by  a  difference  in  the  index  of 
refraction  between  the  prism  region  and  the  surrounding  region. 

Two  geometries  of  prisms  were  designed  into  the  testing  mask.  In  the  first,  a  right 
angle  prism  is  placed  close  to  the  waveguide.  This  prism  minimizes  the  distance  the 
soliton  pump  beam  travels  without  amplification,  however,  it  is  the  most  inefficient  for 
steering.  The  second  uses  an  equilateral  prism  which  does  just  the  opposite,  by  providing 
most  loss,  however,  with  the  best  steering.  If  this  prism  is  operated  into  the  gain  region. 


the  maximum  distance  the  soliton  launch  beam  travels  without  amplification  is  eight 
microns. 

To  permit  independent  control  of  the  soliton  steering,  the  prism’s  electrical  contacts 
are  isolated  from  those  of  the  other  two  regions.  This  isolation  requires  an  electrically 
inactive  gap  on  either  side  of  the  prism.  The  gap,  however,  forms  an  absorption  region 
for  the  soliton  beam  and  must  be  minimized.  Physical  fabrication  limits  impose  the 
tightest  restrictions  upon  this  separation.  Our  fabrication  tools  provided  the  ability  to 
faithfully  reproduce  devices  possessing  a  feature  size  of  one  micron  or  larger. 
Consequently,  the  separation  between  the  single  mode  waveguide  and  the  closest  edge  of 
the  steering  prism  was  selected  at  2  microns  to  permit  high  yield  during  subsequent 
fabrication  runs.  Finally,  a  wide  area  waveguide  forms  the  gain  region  necessary  for  a 
soliton  to  form.  This  region  is  identical  to  the  regions  tested  in  Chapter  4. 

5.5.2  THEORETICAL  DEFLECTION 

Because  etching  a  ridge  around  the  prism  changes  the  upper  cladding  boundary 
conditions,  there  is  a  constant  fixed  shift  in  the  index  of  refraction.  This  is  the  same  shift 
as  formed  by  the  ridge  of  the  single  mode  waveguide.  Appendix  B  describes  methods  for 
calculating  this  index  shift.  In  addition  of  the  fixed  shift,  a  controllable  shift  can  be 
imparted  onto  the  prism  using  either  the  electro-optic,  or  reverse  biased,  effect  where 

(5-4) 

a 
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where  /•41,  the  electro-optic  coefficient,  for  GaAs  is  1.4x10  and  V  is  the  voltage 

applied  across  the  device^^  Since  the  bulk  of  the  field  drops  across  the  active  region, 
d=0.2  microns,  the  expected  change  in  the  index  of  refi'action  is  An  =  3xl0^F.  This 
method  for  steering  the  beam  cannot  be  used  because  the  prism  acts  as  a  strong  absorber 
when  reverse  biased.  Additionally,  imparting  an  index  change  of  0.01  on  the  material 
requires  33  volts  across  the  pn  junction.  This  is  well  beyond  the  reverse  bias  breakdown 
voltage  and  would  destroy  the  device. 


On  the  other  hand,  forward  biasing  the  region  can  reduce  or  eliminate  absorption  and 
induce  a  change  in  the  index  of  refraction  through  both  the  Burstein-Moss  effect  and  the 
plasma  effect,  discussed  earlier.  Phase  shift  measurements  on  a  single  mode  device, 
described  in  chapter  3,  indicate  that  the  shift  in  the  index  will  approximately  be 
A«  =  7.52 X 10“’  J  where  J  is  the  current  density  through  the  device.  Using  this  number 
and  the  right  angle  prism  scheme,  the  angular  for  a  right  prism  with  an  apex  angle  a  will 
be 


AO  =  a  sin 


(n+An) 


n 


sm(90-a) 


-a 


(5.5) 


The  deviation  angle  created  by  the  equilateral  prism  is  given  by^^ 


AO  =  a  Sim 


sm(60°). 


(n^+An) 


K 


-sin^  (30°)  -  sin(  30°  )cos(  60°) 


■30° 


(5.6) 
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assuming  the  apex  angle  of  the  prism  is  60°  and  the  input  beam  is  parallel  to  the  base  of 
the  equilateral  prism.  Decreasing  the  apex  angle  of  the  right  prism  to  very  small  number 
can  yield  large  deviations  for  small  changes.  The  negative  aspect  of  using  this  prism  is 
that  the  distance  the  wave  must  travel  prior  to  entering  the  gain  region  increases 
drastically  for  decreasing  angles.  As  such,  only  the  45-degree  right  and  60-degree 
equilateral  prisms  were  analyzed  during  this  research. 

5.5.3  MEASURED  DEFLECTION 

A  soliton  creation  beam  was  launched  into  this  device  and  the  output  imaged  in  the 
same  manner  as  the  measurements  taken  in  Chapter  4.  In  this  experiment,  however,  the 
single  mode  area,  triangular  prism  area,  and  soliton  area  all  required  separate  contacts 
and  drivers.  The  input  waveguide  and  soliton  creation  regions  were  pulsed,  but  high  to 
destructive  continuous  currents  drove  the  triangular  steering  region.  The  position  of  the 
soliton  measured  against  the  injection  current  in  the  triangular  region  is  shown  in.  Using 
the  electro-optic  data  presented  in  Chapter  3,  the  deviation  angle  versus  change  in 
refractive  index  along  the  theoretical  deviation,  is  shown  in  Figure  5.1 1.  The  departure  of 
the  measured  deviation  from  theory  results  from  leakage  of  current  through  the  sides  of 
the  device  and  in  the  area  leading  up  to  the  prism.  Other  deviations  result  from  the 
current  spreading  away  from  a  strictly  triangular  region  just  below  the  contact  area 
resulting  in  a  “fuzzy”  prism.  Finally,  the  current  density  entering  the  device  is  very  high. 
As  such,  thermal  or  carrier  saturation  effects,  which  increase  the  index  of  refraction,  may 
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limit  the  ultimate  shift  in  the  index.  The  angular  deviation  theory  can  be  fit  to  the  data  by 
including  an  efficiency  term  that  ratios  the  injection  current.  Fitting  to  the  measured  data 
yields  an  efficiency  of  0.33,  which  is  the  solid  line  shown  in  Figure  5.11. 


0  20  40  60  80  100 

Injection  Current  (ma) 


Figure  5.10  Measured  Soliton  Deflection  at  various  injection  currents 


Figure  5.11  Beam  deviation  with  injected  current.  Theoretical  deviations  shown  for  a 

right  angle  and  equilateral  prism. 
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5.6  CRITERIA  FOR  FUTURE  DESIGNS 


The  beam  guided  using  this  experiment  had  a  FWHM  beam  waist  of  8  microns  with  a 
very  multimode  guided  beam.  As  other  experiments  have  shown  nearly  perfect  guiding, 
it  can  be  assumed  that  a  6-micron  soliton  could  guide  a  6-micron  wave.  If  this  distance  is 
used  to  separate  output  waveguides,  this  device  could  theoretically  switch  a  single  input 
beam  to  4  output  waveguides  while  maintaining  a  cross  talk  ratio  less  than  2  to  1. 

The  amount  of  beam  deviation  created  by  the  steering  prism  depends  upon  the  change 
in  the  index  of  refraction  in  the  prism.  Increasing  or  decreasing  the  number  of  carriers  in 
the  prism  area  through  electrical  injection  accomplishes  this  task.  It  should  be  noted, 
however,  that  the  carrier-induced  shift  is  dispersive  in  nature  and  will  steer  the  soliton 
and  signal  beams  differently.  The  amount  of  deviation  the  prism  can  effect  no  more  than 
0.6  degrees,  and  the  variance  in  the  index  of  refraction  between  the  signal  and  soliton 
beam  has  been  calculated  to  be  about  60%.  Using  this  60%,  the  variance  in  angle 
between  soliton  and  signal  will  be  0.2  degrees  and  the  amount  of  energy  confined  in  the 
soliton  beam  will  drop  from  80%  for  a  perfect  launch  to  70%.  Although  the  energy 
confined  will  decrease,  the  increase  in  cross  talk  is  small,  and  the  device  continues  to 
show  promise  as  an  effective  method  of  steering  information  between  points. 

The  angle  of  the  output  soliton  to  the  cleaved  edge  of  the  device  is  the  angle  of  launch 
at  the  input  side.  If  the  information  guided  through  the  device  is  to  be  coupled  to  output 
waveguides,  the  angle  of  the  incoming  information  must  be  matched  to  the  angle  of  the 
output  waveguides.  Using  a  IxN  device  permits  the  output  waveguides  to  be  angled  to 
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match  the  angle  of  the  approaching  soliton.  However  if  an  NxN  device  is  designed,  the 
soliton  could  approach  from  a  variety  of  angles.  Placing  prisms  on  the  input  as  well  as 
the  output  vrill  provide  angular  compensation  for  an  array  of  input  angles.  Note  that  the 
continuous  matching  of  input  and  output  prisms  may  prove  a  daunting  task  for  a 
commercial  device. 

The  maximum  reconfiguration  time  of  the  device  is  the  time  required  to  return  the 
material  to  its  equilibrium  condition.  This  is  the  time  to  pump  carriers  into  the 
conduction  band  to  replace  those  depleted  by  the  soliton  beam.  Measurements  done  on 
semiconductor  optical  amplifiers  indicate  that  the  re-pumping  time  constant  is  on  the 
order  of  1  ns.  It  should  be  expected  that  this  would  be  a  worst  case  time  constant  due  to 
the  fact  that  a  large  pool  of  carriers  spatially  exist  in  the  conduction  band  to  either  side  of 
the  generated  soliton  that  may  return  the  soliton  region  to  near  equilibrium  in  much 
shorter  times.  Thus,  this  would  permit  the  device  to  be  reconfigured  at  a  minimum  rate 
oflGHz. 

Experiments  measuring  these  effects  are  the  subject  of  a  future  research  project. 
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CHAPTER  6:  CONCLUSION 


6.1  OVERVIEW 

During  the  introduction  of  this  dissertation,  the  mathematical  derivation  behind 
solitary  wave  formation  was  reviewed.  This  provided  the  minimum  requirements  that  a 
material  should  possess  to  support  the  possibility  of  solitary  wave  formation.  In  any 
system,  the  major  factor  to  overcome  was  the  amount  of  optical  intensity  needed  to 
produce  a  large  enough  shift  in  the  spatially  localized  index  of  refi'action  to  overcome 
diffraction.  Typical  materials  used  require  optical  intensities  in  excess  of  1  GW/cm^  to 
form  solitary  waves.  Although  new  polymer  materials  possessing  high  nonlinearity 
coefficients  have  been  fabricated  and  tested,  they  tend  to  be  unstable  over  long  periods  of 
time.  In  addition,  materials  compatible  with  current  electronic  circuitry  and  processing 
technology  could  provide  an  advantage  in  future  monolithic  integration  of  electrical 
driving  circuits  with,  optical  sources,  photonic  processing,  and  transmission  circuitry  as 
well  as  optoelectronic  receivers. 

Theoretical  studies  indicated  that  semiconductor  optical  amplifiers  may  provide 
enough  intensity  dependent  refractive  index  changes  to  produce  spatially  self-confined 
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waves.  In  particular,  GaAs  devices  possess  these  optical  qualities  as  well  as  the  potential 
for  monolithic  integration  of  this  material  with  electronic  circuitry  in  a  relatively  mature 
fabrication  technological  base.  Although  time  resolved  studies  of  the  response  of  the 
electronic  nonlinearities  in  the  femptosecond  regime  indicated  that  a  large  negative 
nonlinearity  existed  during  the  first  50  femptosecond  of  the  pulse,  this  nonlinearity 
quickly  vanished  and  a  large,  slower  speed,  positive  self-focusing  nonlinearity  resulting 
both  fix)m  the  direct  influence  of  a  decrease  of  excited  carriers  in  the  device  and  the 
causality  effect  due  to  the  shift  in  the  absorption  edge.  The  penalty  paid,  however,  is  that 
the  speed  of  the  nonlinearity  required  that  the  optical  pulse  creating  this  nonlinearity  must 
last  several  times  longer  in  duration  than  the  period  of  the  initial  negative  going  intensity 
dependent  index  shift,  or  several  tens  of  picoseconds.  Furthermore,  the  index  of 
refraction  shift  created  by  the  incident  solitary  wave  existed  within  the  device  until  the 
depleted  carriers  could  be  replaced  by  electrical  pumping,  commonly  acknowledged  as 
requiring  1  nanosecond.  A  quick  analysis  shows  that  a  device  employing  this 
mechanism  to  transport  information  from  an  input  to  a  corresponding  output  requires 
approximately  1.0  nanoseconds  for  reconfiguration,  allowing  such  a  device  to  be 
reconfigured  at  a  rate  of  1/2  ^gahertz.  As  this  would  be  a  single  stage  switch,  this 
technology  could  provide  a  lOOx  improvement  in  switch  configuration  over  traditional 
fixed  networks  fabricated  using  2x2  couplers  in  a  tiered  network. 

In  our  venture  to  demonstrate  optical  self  focusing  and  beam  guiding  within  GaAs 
semiconductor  optical  amplifiers,  a  characteristic  double  heterostructure  waveguide  was 
designed,  grown,  and  fabricated  using  the  techniques  outlines  in  chapter  2. 
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Characterization  of  these  devices  showed  that  a  strong  optical  nonlinearity  exists  at  the 
peak  of  the  gain  spectrum.  This  not  only  agreed  with  theoretical  calculations  of  the 
device,  but  quantitative  measurements  of  the  nonlinearity  indicated  that  the  effect  was 
strong  enough  to  produce  self  focusing  with  relatively  small  input  optical  intensities. 
Following  this,  a  device  was  fabricated  to  demonstrate  and  analyze  spatial  self  focusing 
and  the  subsequent  successful  creation  of  spatially  self  focused  waves.  After  analysis  of 
the  characteristics  of  these  waves,  a  second  beam  spectrally  separated  from  the  first  one 
by  60  nanometers  was  added  to  the  beam  and  was  found  to  be  spatially  contained  by  the 
index  profile  created  by  the  solitary  wave. 

As  a  demonstration  of  a  device  using  this  effect,  an  electro-optic  steering  prism  was 
added  into  the  device,  and  used  to  steer  the  output  of  the  device  between  an  array  of 
outputs.  These  demonstrations  showed  the  ability  of  the  device  to  direct  an  output  over  a 
30  microns  output  region  using  as  little  as  5  volts  on  the  steering  prism,  thereby 
completing  the  initial  phase  of  the  overall  project  and  this  dissertation. 

6.2  PRIMARY  CONTRIBUTION  AND  MEASUREMENTS 

Specific  measurements  of  relevance  were  made  using  several  characterization  and 
soliton  measurement  techniques.  The  primary  measurements  made  using  single  mode 
devices  produced  the  following  important  conclusions 

A)  The  nonlinearity  follows  the  gain  spectrum  and  the  peak  of  the  nonlinearity  and  gain 

correspond. 


Ill 


B)  The  magnitude  of  the  nonlinearity  when  the  device  was  biased  near  lasing  was 
measured  to  be  n2=5.7xl0'*®  cm^AV  using  a  single  mode  waveguide  and  a  discrete 
Mach-Zehnder  interferometer. 

C)  Spatial  solitary  waves  were  demonstrated  using  an  input  intensity  of  l.lMW/cm*. 
Measuring  the  spatial  profiles  of  these  waves,  the  nonlinearity  was  estimated  to  be 
2.8x10  cm^AV  which  was  in  good  agreement  with  the  value  measured  with  the 
interferometer. 

D)  The  change  in  the  index  of  refraction  was  theoretically  shown  to  be  5.6x10"^  by 
comparing  the  output  mode  profile  of  a  second  wave  guided  by  the  primary  solitary 
wave.However  using  the  nonlinear  refractive  index  found  above  would  give  a  net 
change  in  the  refractive  index  of  1.4x10’^. 

6.3  FUTURE  WORK 

In  the  quest  to  answer  the  basic  questions  concerning  spatial  soliton  creation  in  a 
semiconductor  optical  amplifier,  many  more  questions  have  arisen.  In  an  attempt  to 
promote  and  guide  future  research,  the  following  is  a  list  of  issues  I  feel  should  be 
explored  in  depth  to  answer  the  questions  of  whether  this  technology  is  an  interesting 
scientific  research  topic,  or  a  potentially  viable  commercial  product  technology. 
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6.3.1  SOLITON  LIFETIME 


The  lifetime  of  the  induced  index  of  refraction  created  by  the  soliton  requires 
analysis.  This  information  determines  the  amount  of  time  the  created  waveguide  can 
exist  and  still  be  used,  as  well  as  determining  the  time  the  waveguide  region  requires  to 
completely  recover  and  return  to  the  non-soliton  state.  This  experiment  requires  a  pump 
probe  type  of  measurement  where  the  pump  and  probe  are  synchronized.  The  easiest 
method  initially  considered  is  to  pass  a  pump  and  orthogonal  probe  through  a  single 
mode  active  waveguide  placed  in  the  arm  of  an  interferometer.  This  method,  however, 
may  not  work  due  to  the  large  amount  of  spontaneous  emission  in  the  orthogonal 
polarization  state,  thereby  swamping  the  probe  signal.  A  method  that  may  work  better, 
but  one  that  requires  more  laser  power,  is  to  pass  a  chopped  probe  as  a  signal  beam  just 
behind  the  large  soliton  creating  pump  pulse.  The  polarization  filtered  peak  of  the  probe 
output  would  be  measured  and  analyzed  as  the  probe  is  delayed  further  from  the  pump. 
The  drawback  here  being  that  the  response  of  the  device  would  be  convoluted  with  the 
moderately  long  probe  pulse.  As  the  general  shape  of  the  response  function  can  be 
assumed  to  be  exponential  and  the  temporal  shape  of  the  probe  in  known,  the  two 
functions  may  be  reasonably  deconvolved. 
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6.3.2  WAVEGUK)E  LAUNCH  AND  COLLECTION  EFFICIENCY 


Although  we  have  demonstrated  that  solitons  created  in  a  semiconductor  optical 
amplifier  can  be  used  to  guide  other  information,  the  question  of  launch,  guiding,  and 
collection  efficiency  remains  to  be  analyzed.  In  addition,  the  distance  an  optical  wave 
can  travel  without  gain  and  be  resurrected  into  a  soliton  should  be  analyzed.  In  this 
experiment,  the  goal  is  to  analyze  the  difference  between  waveguiding,  launch,  non-gain 
propagation  loss,  and  soliton  carrier  loss.  In  changing  the  variables  of  input  waveguide 
length,  non-gain  spacing,  and  soliton  region  spacing  each  of  these  effects  may  be 
decoupled.  It  should  be  stressed  that  sample  to  sample  deviations  force  as  many  variants 
of  the  structure  to  be  fabricated  onto  the  same  sample  to  ensure  uniformity  between 
devices. 

6.3.3  SOLITON  INTERACTION 

In  a  more  advanced  NxN  device,  multiple  solitons  may  be  crossing  the  gain  region, 
and  each  other,  simultaneously.  As  the  address  of  the  output  from  the  input  may  be  a 
fixed  voltage  on  an  electro-optic  prism,  it  is  critical  to  determine  the  effects  of  crossing 
solitons  on  the  target  stability  of  the  output.  In  this  experiment,  shown  in  figure  6.3,  two 
solitons  are  launched  into  the  same  area  and  set  to  cross  at  varying  angles.  The  soliton 
interaction,  angular  deviation  caused  by  index  profiles,  or  possibly  trapping  (although 
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this  seems  unlikely  for  a  non-instantaneous  nonlinearity),  should  be  analyzed  on  the 
output  of  the  device. 

6.3.4  CROSS  TALK  (IXN)  DEVICE 

Cross  talk  between  crossing  solitons  and  in  adjacent  outputs  will  drive  the  packing 
density  of  the  device.  As  such,  experiments  must  be  performed  to  determine  the  cross 
talk  isolations.  This  experiment,  shown  in  figure  6.4,  uses  the  same  a  complete  2x2 
device  with  varied  input  and  output  spacing.  Two  spacings  may  be  sufficient  to 
determine  cross  talk  vs.  distance,  however,  at  least  three  varying  spacings  will  help 
confirm  trends  in  data. 

6.3.5  COMPENSATING  RECEPTION  ANGLES 

In  a  device  with  one  input  and  N  outputs,  each  of  the  outputs  can  be  angled  precisely 
toward  the  apex  of  the  launching  electro-optic  prism  ensuring  the  optimum  reception 
efficiency  possible.  In  the  case  with  multiple  inputs  and  outputs,  however,  the  receiving 
waveguides  cannot  be  angled  at  all  inputs  simultaneously.  Thus,  the  placement  of 
reception  prisms  at  the  receiving  waveguides  can  help  compensate  the  variance  in 
reception  angle  deviation  fi-om  the  angle  of  the  reception  waveguide.  Factors  to  be 
determined  are  whether  the  angles  can  be  compensated  and  what  efficiency  can  be 
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achieved,  as  well  as,  what  precision  of  electro-optic  control  voltage  is  necessary  to 
achieve  a  minimum  signal  to  noise  ratio. 

6.3.6  FULL  DEMONSTRATION  DEVICE  (NXN) 

Eventually  most  of  the  critical  design  issues  will  be  answered,  and  a  device  complete 
with  multiple  inputs  and  outputs  needs  to  be  demonstrated.  A  mass  extension  of  this 
device,  shown  in  figure  6.6,  incorporates  an  8x8  geometry.  Each  launching  waveguide 
requires  a  steering  electro-optic  prism,  as  does  each  receiving  waveguide.  Aside  from  the 
ability  to  be  used  bidirectionally,  the  steering  electro-optic  prisms  on  each  side  of  the 
soliton  region  will  compensate  for  variances  in  reception  angles  caused  by  the  differing 
launch  positions. 
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APPENDIX  A:  DERIVATION  OF  THE  NON-LINEAR  WAVE  EQUATION 

FOR  SPATIAL  SOLITONS 


In  developing  the  non-linear  wave  equation,  we  start,  of  course,  with  Maxwell’s 
equations  for  non-magnetic  materials,  which  are 


Vx£  = 


VxH  =  J  + 

VD  =  p 
VB  =  0 


a 


D  =  e^E  +  P  «  =  sE 

b  =  m,h 


An  optical  field  is  commonly  expressed  as  a  fimction  of  a  traveling  envelope  (A  in  this 
equation)  and  a  traveling  phase  wave. 


E(x,y,z,t)  ^^A(x,y,z,t)e^^^~^^  +c.c. 
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The  first  step  in  developing  the  wave  equation  is  to  take  the  curl  of  the  first  Maxwell’s 


equation  (and  assume  a  source  free  medium  ( v  -  D  =  0 )). 

Vx:Vx£  =  VxC%) 
V(V  E)-V^E  =  ~ 


Although  we  set  V'D  to  zero  for  a  source  fi’ee  medium,  the  assumption  that  the  V*E 
term  is  also  zero,  in  some  cases,  is  not  valid.  Assuming  that  s  contains  polarization  effects 


and  varies  spatially 


D  «  fiE 

VD  =  V-(fi£:) 

0  =  Ve-E  +  fiVE 
eWE  =  -WeE 


VE 


This  tells  us  that  if  there  is  a  gradient  in  the  permitivity,  or  index  of  refraction,  the  term 
V’E  is  not  zero.  Many  waveguide  structures  such  as  graded  index  waveguides  or  third 
order  nonlinear  materials  hosting  solitary  waves  have  such  a  gradient.  Fortunately,  the 
change  in  permittivity  is  very  small,  and  in  most  cases  can  be  totally  ignored,  as  we  shall 
do. 
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V=£  =  «>^Vxff^ 

/^^(7  + 

V^E  = 


where 

c  is  the  speed  of  light  in  a  vacuum 
J  is  assumed  approximately  constant  with  time 

If  this  was  free  space  (P=0),  the  final  equation  above  would  be  the  Helmholtz  wave 
equation. 


We  deal,  however,  in  polarizable  material.  In  the  time  domain,  the  polarization  is  the 
convolution  of  the  nonlinear  response  with  the  input  electric  field. 


it)  =  ^0 1 x^\t  -  r)Ejir)dv 

—00 

=  £01/ h)Ejir,)E,ix^)dv,dT^ 

—00—00 

f‘*'(0  =  j  j  i -  ’■2.'  -  t,)E,{r,)E,(v^)E,(z,)ffr,gT,d^, 
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where 

units  are  in  MKS 

jl^  is  the  susceptibility  of  the  material  and 
is  the  x'*  order  non-linearity 


One  way  to  handle  the  complexity  of  these  convolutions  is  to  work  with  the 
polarization  in  the  frequency  plane,  where  the  convolutions  of  time  functions  become 
multiplication  of  the  frequency  functions.  A  second  method  is  to  assume  that  the 
susceptibility  responds  much  faster  than  the  electric  field  and  can  be  approximated  as  an 
instantaneous  function  in  time.  This  reduces  to  polarization  to 


P'”‘(>)  =  e.z*“E,.(/)E,(i)E,(0 

This  equation  can  again  be  simplified  by  assuming  a  specific  set  of  experimental 
conditions.  If  we  are  interested  only  in  one  input  field,  we  set  a  fixed  and  specific  crystal 
orientation,  and  we  analyze  only  the  fundamental  frequency  output  (self-action  term),  we 
can  replace  the  susceptibility  tensors  with  effective  susceptibilities  valid  only  for  our  given 
experimental  setup.  This  can  be  justified  in  semiconductor  waveguide  experiments  where 
the  waveguide  is  oriented  on  a  specific  crystal  direction  and  only  the  light  that  is  confined 
to  the  waveguide  is  of  interest.  Although  the  total  polarization  is  a  complex  function,  the 
reduced  equation  is  a  power  series  driven  by  the  applied  field. 
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P  =  s,(x"'E  +  +  X'^'E’  +  x'‘'E'^■..) 

where  is  the  n***  order  effective  susceptibility  for  our  specific  setup  (MKS  units). 
Inserting  this  into  the  Maxwell’s  equation  above 


In  the  class  of  materials  we  are  concerned  with  is  called  centro-symmetic,  third  order 
nonlinear.  This  means  that  the  material  is  electronically  symmetric  about  the  center  of  the 
molecule,  thus  the  even  order  non-linearities  cannot  exist.  The  odd  order  nonlinearities 
can  exist.  Of  these  nonlinearities,  only  the  third  order  nonlinearity  vrill  be  significant, 
unless  the  field  strength  is  extremely  high.  In  this  case  the  wave  equation  reduces  to 

y 2  p _ ^  ^ K/  1  woy (')  ^ E/  I  ^ E  / 

The  term  contains  the  cube  of  the  electric  field.  Cubing  the  field  defined  at  the 
beginning  of  this  paper 
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Eix,y,z,tf  =^A{x,y,z,t)e^^"*-^'> 

A{x,y,z,(fe^^''“*~'^'*  +  'iA{x,y,z,iy  A{x,y,z,t)  *  + 

}A{x,y,z,t)e^^‘^-^'>A{x,y,z,t)  +  Aix,y,z,ty^e-^^^‘^-^^ 

=  <)’«’■'*"■'”  ■K:.c]  +  ^^{x,y,z,(f  A(x,y,!,t)e'^“-‘‘' 

=  i  [^(x,;<,r,  ,)’««»-■=>  +  c.  c.]  +  \  |£f  E 


+  C.C. 


1 


In  our  case,  we  are  interested  only  in  the  self-action  terms,  or  the  terms  that  result  in 
the  same  frequency  as  the  input  term.  The  |  E  ]  ^  term  is  estimated  to  be  independent  of 
time  and  can  be  removed  from  the  derivative.  This  term  is  also  directly  related  to  the 
positional  intensity  1=  j  E  |  where  )/  =  Vp/e  or  l/(ceono).  Regrouping,  the  equation 
now  becomes 


Now,  the  interesting  term  is  a  unitless  term  that,  in  effect,  ratios  the 

speed  of  the  optical  wave  in  a  material  to  that  of  the  wave  traveling  in  free  space.  This 
might  be  recognized  as  the  square  of  the  optical  index  of  refraction  n.  It  is  common  to 
represent  the  index  of  refraction  as  a  linear  and  nonlinear  component  of  the  form  n(I)  =  no 
+  n2l.  To  find  the  value  of  n2,  we  square  the  common  form  equation  n(I)^  =  no^  +  2non2l  + 
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n2^I^.  Because  n2  is  a  very  small  number,  we  can  ignore  the  final  term  as  insignificant  and 
solve  the  simple  equation 


nl  +  2n^J  =  1  +  /  2 


A3). 


fio=yli  +  Z 


(1) 


«2  = 


(3) 


4w„ 


\ce^l 


w(/)«w„  +V 


Finally,  the  wave  equation  becomes 


n(jf  /fF.  y 


'a'‘ 


=  0 


Now  let’s  start  executing  the  derivatives.  We  assume  that  the  spatial  and  time 
variables  are  separable.  Additionally,  we  choose  to  analyze  the  wave  in  a  frame  moving 
with  the  envelope,  the  time  dependence  of  the  envelope  drops  off  the  and  second 
derivative  with  time  produces  the  original  function  times 


'^^E+n{lfklE^Q 

k  =^/ 

*0  /c 
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The  Laplacian  represents  the  spatial  derivatives.  The  Laplacian  in  rectangular 


coordinates  is 


Referring  to  the  form  of  the  field  stated  at  the  start  of  the  derivation,  the  x  and  y  terms 
are  independent  and  will  be  referred  to  as  the  transverse  coordinates.  The  Laplacian  for 
the  transverse  coordinates  will  be  the  symbol  V^j^.  The  z  term  appears  both  in  the 
amplitude  and  traveling  terms  and  the  Laplacian  of  the  z  term  is  (dropping  the  non  z 
dependent  terms) 

The  first  and  second  term  of  this  equation  contain  the  term  k  which  is  the  wave 
propagation  constant  in  the  z  direction  (k  =  nko=  n27r/A).  Values  for  k  in  the  optical  range 
are  large.  For  a  near  infra-red  beam  (800  nm)  in  glass  (n=l .5),  k  =  1 . 178  •  lO’  m■^  Thus 
any  term  without  k  will  be  quite  insignificant,  unless  the  pulses  are  short. 
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Why  short?  To  examine  the  approximate  magnitude  for  the  derivative  terms,  lets  take 
a  snap-shot  of  the  pulse  (again,  use  a  moving  reference  frame)  and  set  the  traveling 
envelope  function  A(z,t)  in  terms  of  z.  For  a  typical  short  laser  pulse,  the  function  is 
gaussian,  or  very  close  to  it,  which,  at  any  given  time,  has  the  form 


-i—J 

A(z)  oce  ^  ^ 

where 

z  is  the  dependent  variable 

zo  represents  the  displacement  along  the  z  axis,  0  for  a  mo\ing  reference  frame 
w  represent  the  relative  spatial  width  of  the  pulse 

Laser  pulses  vary  greatly  in  temporal  pulse  width,  but  pulses  from  active  mode  locked 
lasers  are  typically  1-100  picoseconds.  Kerr  lens  mode  locked  Titanium  Sapphire  laser 
pulse  widths  can  be  as  short  as  6  femptoseconds,  but  typically  range  from  100-3000 
femptoseconds.  The  relationship  between  temporal  width  and  spatial  width  (w)  is  given 
by 


w 


where 

to  is  the  e‘*  temporal  pulse  width 
c  is  the  speed  of  light 
n  is  the  optical  index  of  refraction 
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For  a  1  picosecond  pulse  in  glass  (n=1.5)  w=200  microns.  Now  examine  the 


derivatives  of  A(z) 


The  (z-Zo)  term  in  the  above  equations  is  localized  about  the  Gussian  function  and  is  of 
the  order  of  w.  Thus  the  first  derivative  has  a  maximum  on  the  order  of  1/w,  or  in  the 
case  of  a  1-picosecond  pulse,  5000  m‘^  The  second  derivative  has  a  maximum  on  the 
order  of  1/w^  or  2.5  •  lO’ m'^.  Now  looking  at  the  fiill  equation  ( i.e.  multiplying  by  and 
k),  the  magnitude  the  third  term  is  still  less  than  three  orders  of  magnitude  less  than  the 
second  term.  Thus,  the  second  derivative  can  be  ignored.  This  approximation  is  also 
called  the  slowly  varying  envelope  approximation,  meaning  that  the  envelope  of  the 
function  varies  much  slower  than  the  frequency  of  the  electromagnetic  wave.  Substituting 
the  terms  for  A(z)  and  its  first  derivative  back  into  the  wave  equation 
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'^^E  +  kln{lfE  =  0 
V\A-ljk^y^-eA  +  kln{If  A  =  0 
^{A-  2jk  -  nXA  +  kln{lY  A  =  0 

VIA- 2jk ^y^  +  kl{n{,lf  - nl)A  =  0 


We  can  further  simplify  this  equation  by  noting  that  n2l  is  very  small  in  comparison  to 
no  so  that  the  index  square  terms  can  be  approximated  by 


n(/y  - nl  =  (w(/)  +  n^Ynil)  - n^) 

=  (w,  +  n^I  +  «oX«o  +  -  Wo) 

=  (2«oXW2/) 

=  (2Wo)w2 

=  n^n^\Aflrj 

Substituting  this  back  into  the  wave  equation 


VIA- Ijk %  +  kl{r^If  - =  0 
VlA-2)k^y^  +  kl^\Af  A  =  I3 

V\A*kl’^\Af  A  =  2jk^/^ 
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If  taken  in  a  1 -dimensional  waveguide,  this  equation  has  an  exact,  steady  state. 


solution,  called  a  spatial  soliton,  of  the  form 


^(x,z)  =  4)  sec/» 


I{x,z)  = 


2rj  2rj 


where 
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APPENDIX  B:  KRAMERS-KRONIG  RELATIONSHIP 


The  polarizability  of  a  material  is  commonly  written  as  a  Taylor  expansion  about  the 
applied  electric  field 

P  =€o  +X^'^E^  +...) 

The  susceptibility  tensors  %  are  frequency  dependent  have  both  real  and  imaginary 
parts.  In  simple  situations,  such  as  linear  materials  Mrith  only  terms  or  in  Kerr  like 
materials,  the  real  part  of  the  susceptibility  relates  to  the  index  of  refraction,  and  the 
imaginary  part  relates  to  the  absorption  or  gain.  The  Kramers-Kronig  relationship  defines 
a  specific  relationship  between  the  real  and  imaginary  parts  of  this  or  any  real  function. 

As  with  any  function  defined  in  the  frequency  domain,  the  Fourier  transform  links  this 
function  to  a  corresponding  time  domain  function.  It  is  this  function  that  forms  the  basis 
for  the  Kramers-Kronig  Relationship.  The  causal  principle  states  that  the  response  to  a 
stimulus  cannot  occur  prior  to  the  application  of  the  stimulus.  In  our  case,  a  polarization 
response  to  an  electric  field  cannot  occur  prior  to  the  application  of  the  field.  Assume  that 
the  stimulus  is  applied  at  time  zero  and  the  response  to  the  stimulus  is  an  arbitrary  function 


129 


f(t).  By  the  causal  principle,  f(t)  must  be  zero  for  times  less  than  zero.  As  such,  f(t) 
multiplied  by  a  step  function  at  zero  is  stiliy(t)  or  / (t)  =  fif)  *  step(0) . 

The  frequency  domain  response  is  found  be  executing  the  Fourier  transform. 

F((i))  =  3{/'(r)}®  3{step(0)} 

—00 

where  the  integral  is  the  execution  of  the  convolution  function.  Now,  allowing  F(a)) 
to  have  both  real  and  imaginary  parts,  F(m)  =  F'(<u)  +  iF  "(co) ,  equation  3  becomes 


F'(g))  +  F"(®) 


F'((o')  +  iF"(o)) 

CD -O) 


dcD' 


Real  and  imaginary  parts  of  an  equality  must  be  independently  equal,  thus 

— oO 

F''(®)=  hcD 

^a-CD 


Thus,  the  real  and  imaginary  parts  of  the  frequency  response  are  related  through  the 
Kramers-Kronig  integral.  In  practice,  this  integral  is  relatively  easy  to  graph  by  noting  that 
the  integral  is  the  convolution  of  the  real  or  imaginary  part  of  the  response  with  the 
function  within  the  integral. 
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As  an  example,  apply  the  Kramers-Kronig  relationship  to  a  Gaussian  response 
function,  shown  in  figure  2.  As  the  integral  function  approaches  from  the  low  frequency 
side,  the  convolution  wall  follow  the  same  sign  as  the  function.  As  the  convolution  nears 
the  center  of  the  Gaussian,  the  convolution  turns  towards  the  axis,  and  hits  the  axis  when 
the  two  functions  have  the  same  center  point.  As  the  integral  function  proceeds  further, 
the  result  is  negative,  then  turns  towards  the  axis  as  the  two  functions  separate. 

An  important  corollary  to  this  equation  is  that  changes  in  the  absorption  or  gain  can  be 
directly  related  to  changes  in  the  refractive  index  directly  wdthout  measuring  the  entire 
absorption  spectra. 

^^(0-0)^ 

—  00 


If  the  change  in  the  absorption,  or  gain,  is  caused  by  an  incident  optical  intensity,  then 
this  equation  calculates  the  change  in  optical  refractive  index,  nj . 

Although  powerful  in  its  simplicity,  the  greatest  drawback  to  this  equation  is  that  to 
produce  an  exact  form  of  the  change  in  the  index,  the  change  in  the  absorption  must  be 
known  throughout  the  entire  spectrum.  Often  this  is  either  impossible  or  impractical. 
Fortunately,  the  change  to  the  index  at  a  given  frequency  drops  dramatically  as  the 
absorption  measurement  range  fiirther  fi'om  the  wavelength  of  interest.  Therefore,  this 
method  can  be  used  to  determine  a  qualitative  measure  of  the  change  in  the  index  of 
refi'action  given  a  limited  spectrum  measurement  of  the  change  in  the  absorption. 
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APPENDIX  C:  WAVEGUIDE  MODE  CALCULATIONS 


Electric  fields  traveling  through  a  waveguide  must  conform  to  continuous  boundary 
conditions.  The  first  condition  mandates  that  the  fields  decay  to  zero  at  large  distances 
away  fi’om  the  core  of  the  waveguide.  Secondly,  the  field  in  the  core  may  have  oscillatory 
behavior.  Finally,  that  the  normal  components  of  the  electric  flux  density  must  be  equal 
(Z)„,  =D„2),  and  that  the  tangential  electric  field  and  its  derivative  must  be  equal 


i^n  = 


dx 


\ 

dx 


These  assume  zero  charge  at  the  boundaries. 


With  the 


coordinate  system  set  at  x=0  for  one  boundary,  and  x=h  at  the  other,  the  solution  of  this 
system  is 


f(x)  =  e’^ 

(  s  \ 

f(x)-  cos(-)Kx)-\-—sin(Yx) 

\  K  ; 


f(x)  = 


cos(  ydi)  +— sm(  ydi ) 

K 


g  ■p( x-h)^H<»t-Pz) 


lower  cladding  region 
guiding  region 

upper  cladding  region 


The  relationship  between  the  decay  and  oscillation  constants  and  the  propagation 
constant,  p,  is  found  by  substitution  of  these  equations  into  the  SWE  discussed  in 
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appendix  A.  Application  of  the  boundary  conditions  to  these  equations  yields  the 
dispersion  relationship 

^  ,  n .  K(p-\-s) 

tcm(Kh  +  m—)  =  -^ - 

2  K  -ps 

m  in  this  equation  is  the  order  of  the  waveguide  mode. 

This  solution  applies  only  to  configurations  where  the  cladding  layers  are  infinite.  If 
the  cladding  is  not  infinite,  such  as  in  the  case  of  the  weighted  waveguide,  an  additional  set 
of  boundary  conditions  must  be  applied.  The  math  quickly  become  tedious  and 
computerized  solution  methods  exist  to  determine  the  propagation  mode  constants. 

A  lateral  transmission  method  can  be  used  to  calculate  propagation  constants  for  a 
structure  of  any  number  of  layers.  The  method  uses  the  multiple  beam  interference 
method  to  calculate  resonances  in  a  multiple  layer  stack.  Light  at  an  interface  of  high 
index  to  low  index  will  be  total  internally  reflected  for  all  angles  greater  than  the  total 
internal  reflection  angle.  Certain  angles,  however,  will  cause  the  multiply  reflected 
continuous  waves  within  the  waveguide  to  overlap  constructively.  These  are  the  modes  of 
the  waveguide.  At  these  specific  angles,  there  exists  a  singularity  in  the  angular 
transmission  spectrum.  Locating  the  angle  of  these  singularities  yields  the  propagation 
constant  through  the  expression  P  =  A:„  sin(Q),  where  A:„is  the  propagation  constant  in 

the  bulk  material.  The  angles  are  found  by  dividing  the  structure  into  a  series  of  regions 
and  applying  a  matrix  method  to  determine  the  series  solutions  for  the  structure.  The 
matrix  consists  of  an  interface  step  defined  by  the  matrix  Hn  (interface  between  layers  1 
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and  2)  and  a  propagation  matrix  step  defined  by  the  matrix  L2  (propagation  across  layer 

2). 


11 


Pv 


where  p^^  andTj,  are  the  reflection  and  transmission  coefficients;  kj  is  the  material 

propagation  constant. 

The  total  transfer  matrix  is  then 


^12^2^23^3  n-\,n 


Sn  S,, 

s^, 


where  the  angular  dependence  of  the  equation  is  found  using  Snell’s  Law.  The  final 
transmission  and  reflection  of  this  stack  is  then 


To  implement  this  method,  an  angular  ray  at  the  first  level  of  the  waveguiding  region  is 
used  to  find  the  series  of  angles  throughout  the  rest  of  the  system.  These  angles  are  used 
to  calculation  the  transmission  and  reflection  coefficients  for  the  H  and  L  matrix  elements, 
and  the  S  matrix  is  calculated.  The  initial  launch  angle  is  varied  to  approach  the 
singularities  of  the  waveguide  modes. 

A  Matlab  code  implementing  this  method  appears  at  the  conclusion  of  this  appendix. 


134 


function  duinb=inain; 

% 

%  index  =  index  of  refraction  of  each  layer 
%  thick  =  thickness  of  each  layer  in  meters 
%  guide  =  the  layer  the  angles  will  be  set  in.  All 

%  other  layer  angles  depend  on  the  guide  angle 

%  lambdaO  =  wavelengh  in  free  space 

%  counts  =  number  of  wavelengths  that  will  be  analyzed 

%  Note:  the  angle  range  is  set  by  the  pTh 

%  statement  inside  the  counts  loop. 

% 

%  Program  calls 

%  Angles  ->  calculate  all  angles  (complex)  within  the 
%  stack 

%  Stack  ->  calculate  transmission  and  reflection  of  the 
%  stack 

um=le-6; 

nm=le-9; 

ang=le-10; 

deg=2*pi/360; 

%lambda_start=500*nm; 

% 1 ambda_end= 1 2 0 0 * nm ; 

% 1 ambda_inc= 1 *nm ; 

%lambda0=1000*nm; 

%for  mn=l : 30, 

%layers=mn*2 ; 

%layers=4 ; 

%n_top=l ; 

%n_low=l .5; 

%n_high=2 .0; 

%n_sub=3 .5; 

% index ( 1 ) =n_top ; 

%for  1=2 ; 2 : layers+1 
%  index ( 1 ) =n_high ; 

%  index (1+1) =n_low; 

%end; 

%index (layers+2) =n_sub; 

%k=round( (lambda  end-lambda_start) /lambda_inc) 
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%index=[3.59  3.08024  3.45164  3.590  3.45164  3.590  3.45164 
3.590  3.45164  3.08024  3.590  1]  ; 

%thick=[l  1.39*uiti  0.075*tim  80*ang  80*ang  80*ang  80*ang 
80*ang  0.075*um  0.93*um  1000*ang  1]; 

%index=[3.59  3.08024  3.453164  3.590  3.453164  3.08024]; 
%thick=[l  1.39*um  0.7385*um  0.008*um  0.7385*um  1]; 

index=[3.1  3.58  3.1]; 
thick=[l  0.2*um  1]; 

lainbda0=900*nm; 

%index=[1.40  1.50  1.45]; 

%thick=[l  5*uia  1]; 

%  lainbda0=  1 0  0  0  *nin; 

guide=2; 
counts=2000; 
j=sqrt (-1); 

%theta=0; 

for  i=l: counts, 

pTh (i) = (90- (30* i/counts) ) *deg; 

[rho, tau,R, T]=stack(lambdaO, index, thick, pTh (i) , guide) ; 
pR (i) =abs (R) ; 
pT (i) =abs (T) ; 
end; 

subplot (1,2,1),  plot (pTh/deg, pR, ' r ' ) ; 

%axis([60  90  -0.1  2000]); 

subplot (1,2,2),  plot (pTh/deg, pT, ' g ' ) ; 

%axis([60  90  -0.1  2000]); 

itiodeangles=  [  ] ; 

i=l; 

while  i<counts, 
if  pR (i) >5000, 

low=pTh (i-1) /deg; 
while  pR(i)>5000, 
i=i+l; 
end; 

high=pTh (i) /deg; 

modeangles= [modeangles, (low+high) /2] ; 
end; 
i=i+l ; 
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end; 


kf=2*1.5*pi/lainbda0; 
kx=cos (modeangles*deg) .*kf; 

ef f_index=index (guide) *sin (modeangles (1) *deg) ; 
kx 

modeangles 
ef f_index 
end; 
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function 

[Srho, Stau, R, T] =s tack (lambda, index, thick, theta, seed) / 

thetas=angles (index, theta, seed) ; 

betas=cos (thetas) .*index.*thick* (2*pi/lambda) / 

i=sqrt (-1) ; 
k=length (index) -1; 

a=tan (thetas (k) ) /tan (thetas (k+1) ) ; 
tau=2/ (1+a) ; 
rho= (1-a) / (1+a) ; 

S=[l  rho  ;  rho  1] ./tau; 

for  j=k:-l:2, 

a=tan (thetas ( j-1) ) /tan (thetas ( j ) ) ; 
tau=2/ (1+a) ; 
rho=(l-a)/(l+a) ; 

H=[l  rho  ;  rho  1] ./tau; 

L= [exp (-i*betas ( j ) )  0  ;  0  exp (i*betas ( j ) ) ] ; 

S=H*L*S; 

end; 

Srho=S(l,2)/S(2,2)  ; 

Stau=l/S(2,2) ; 

R=abs  (Srho)  ''2; 

T=  (abs  (Stau)  '"2)  *  (index  (k+1)  *cos  (thetas  (k+1)  )  )  /  (index  (1)  *c 
os (thetas (1) ) ) ; 
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function  angs=angles (index, theta,  seed) ; 


seed  is  the  layer  number  that  is  calling  all  other  angles 

l=length (index) / 
angs (seed) =theta; 
i=sqrt (-1) ; 

for  j=seed+l:l,  %layers  after  seed 

aratio=index ( j-1) *sin (angs ( j-1) ) /index ( j ) ; 
iangs=imag (asin (aratio) ) ; 
rangs=real (asin (aratio) ) ; 

if  abs (aratio) >1 
angs ( j ) =r angs -i angs *i; 
else 

angs ( j ) =rangs+iangs*i; 
end; 

end; 

for  j=seed-l : -1 : 1,  %layers  before  seed 

aratio=index ( j+1) *sin (angs ( j+1 ) ) /index ( j ) ; 
iangs=imag (asin (aratio) ) ; 
rangs=real (asin (aratio) ) ; 

if  abs (aratio) >1 
angs ( j ) =rangs-iangs*i; 
else 

angs ( j ) =rangs+iangs*i; 
end; 

end; 


function  n=kzeros (kappa) ; 
global  nf  ns  nc  h  lambda  k  deg; 

t=tan (kappa*h) ; 

betas= ( ( (k^2) * (nf^2) ) - (kappa^2) ) ; 
GS=sqrt  (betas-  (  (k^2)  *  (ns''2)  )  )  ; 
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GC=sqrt  (betas-  (  (k''2)  *  (nc''2)  )  )  ; 

tr= (GS+GC) / (kappa* (1- (GS*GC/ (kappa^2) ) ) ) / 

n=t-tr; 

end; 
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APPENDIX  D:  DERIVATION  OF  THE  BEAM  PROPAGATION  METHOD 


A  very  powerful  method  of  analyzing  a  beam  propagating  through  a  complicated 
material  is  the  beam  propagation  method  (BPM).  The  method  is  particularly  powerful  in 
circumstances  where  a  closed  form  solution  is  not  available.  Additionally,  the  method 
provides  insight  into  the  evolution  of  steady  state  functions  if  those  functions  are  not  input 
directly  to  the  structure.  To  show  the  power  the  beam  propagation  method,  we  start  with 
the  wave  equation  and  expand  the  polarization  about  the  fundamental  frequency,  using 
effective  susceptibilities,  for  centro-symmetric  (no  even  order  nonlinearities),  third  and 
fifth  order  nonlinear,  non-magnetic  materials. 

P  =  +  x^'^^EfE  +  x^'^eH) 


Combining  these  equations,  yields 
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_.^  (u;ir<"+/’>W  +  y’W)^£/ 

V£=  7^2 


where  we  substituted  n^  for  the  susceptibility  term.  The  term  n,  which  is  a  function  of  the 
electric  field  intensity,  should  not  be  confused  with  the  linear  index  of  refi’action,  no,  which 
is  not  a  function  of  the  field  intensity.  We  assume  a  trial  solution  to  this  equation  as  a 
standard  propagating  harmonic  wave  in  the  z  direction  with  a  spatial  envelope  function. 


E(x,y,z,t)  =  UA(x,y,z)e“'“^‘‘>  + 


where 

A  is  the  spatial  envelope  of  the  wave 
/y  is  the  frequency 
k  is  the  spatial  wave  vector 

Inserting  this  into  the  wave  equation  and  assuming  the  z  dependency  of  the  envelope 
function  is  separable  from  the  x  and  y  dependencies  yields 
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S/^^A-k'^A-2jk^y^+^^A  =  Q 

-  2jk  +  (n^k^^  -  k^)A  =  0 

V,='^-2#%  +  («V  -«oV)^  =  0 

-  2jk  %  +  {re  -  nl)k,^A  =  0 

We  have  dropped  the  second  derivative  of  A  vith  respect  to  z  because  this  term  is 
insignificant  The  term  is  the  transverse  Laplacian,  or  the  Laplacian  in  the  x  and  y 
coordinates  only.  As  this  equation  is  written  there  is  no  closed  form  solution  for  A. 
Additionally,  the  problem  of  solving  for  A  becomes  even  more  complicated  if  there  is  a 
variation  in  the  index  of  refi'action  of  the  material  across  the  wave,  as  in  a  waveguide 
structures.  Using  the  beam  propagation  method,  however,  we  can  approximate  the  form 
of  A  at  a  small  increment  in  the  z  direction  beyond  a  known  A.  To  see  this,  rewrite  the 
equation  using  functional  notation,  isolate  the  envelope  function  and  solve  the  derivative 
at  a  small  distance  Az  beyond  a  known  envelope  assuming  the  index  terms  are  constant 
across  the  Az  interval. 
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+  (w^  -nl)K^ A 

"  ijk^’  ■■"«)*«’)^ 

ln(,4(r  +  Az))  -  ln(^(z))  =  +  (”'  "  «!)*•* +  Az  -  z) 

^(z  +  Az)=«^'''’*‘"’''‘’‘‘’'“-4W 

— JVj^*i4(z}4-(n^-»io)to^/4{z))Az 


^(z  +  Az)  =  e 


- 


^(z) 


This  result  shows  that  the  known  field  and  the  material  parameters  can  approximate 
the  field  envelope  at  a  small  incremental  distance  from  a  known  point.  Thus  the  envelope 
function  can  be  determined  by  repeated  application  of  the  formula.  The  only  remaining 
obstacle  is  the  evaluation  of  the  exponential  term,  particularly  the  transverse  Laplacian. 
Although  we  might  have  a  shot  of  evaluating  the  Laplacian  for  a  real  starting  field,  after  a 
few  increments,  the  complexity  of  the  envelope  function  would  overwhelm.  An  alternative 
is  to  specify  the  envelope  at  a  known  point  along  the  ajds  of  transmission  with  a  sampling 
matrix  giving  the  known  field’s  amplitude  and  phase  in  the  transverse  plane.  This  also 
allows  us  to  specify  the  index  of  refi'action  as  a  sampling  matrix.  Using  this  method,  the 
complex  envelope  function  is  easily  tracked  through  the  real  and  imaginary  amplitudes  in 
each  propagated  sampling  matrix.  An  additional  benefit  to  using  a  sampling  matrix  is  the 
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evaluation  of  the  transverse  Laplacian,  which  is  the  second  partial  derivative  in  rectangular 
coordinates.  Using  Fourier  theory  in  one  dimension 


Vl/ W  = 


:r'{*,‘3{/(x)}) 


where 

^  is  the  Fourier  transform  operator 

is  the  inverse  Fourier  transform  operator 
kx  is  the  spatial  frequency  along  the  x  direction 

The  Fourier  transform  can  be  evaluated  using  the  discrete  Fourier  transform  (DFT), 
usually  implemented  in  most  codes  as  the  fast  Fourier  transform  (FFT).  The  DFT  is  a 
summation  approximation  of  the  integral  form  of  the  transform 


m) = 3{/w} = 


,  2;r(5-l) 

=  - 7; - 


j2!t(,s-\yir-\yN 


where 

f(x)  is  a  spatial  function 

Ffitx)  is  the  spatial  frequency  function  of  f(x) 

kx  is  the  spatial  frequency 
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fr  is  the  sampled  spatial  function  of  f(x) 

Fs  is  the  DFT  spatial  frequency  function  of  f, 

A^is  the  total  number  of  points  in  the  sample 
H  is  the  spatial  extent  of  the  sample 
5  is  the  sample  point  index  in  the  Fourier  plane 
r  is  the  sample  point  index  in  the  real  plane 

An  assumption  in  the  discrete  Fourier  transform  is  that  the  input  waveform  repeats 
itself  in  space,  with  a  period  of  the  spatial  extent  of  the  sampling.  Figure  D.l  demonstrates 
this  problem.  Each  line  in  the  figure  represents  a  BPM  propagation  step.  Note  that  each 
step  requires  the  execution  of  a  DFT  and  an  inverse  DFT,  both  assuming  that  the  function 
is  periodic  beyond  the  sampling  extent.  As  portions  of  the  wave  approach  the  edge  of  the 
sampling  space,  through  repeated  application  of  the  DFT,  a  matching  wave  approaches  the 
opposite  side  of  the  sampling  extent.  Waveform  components  that  travel  across  the 
boundary  of  the  spatial  extent  of  the  sampling  will  reappear  on  the  other  side  of  the  spatial 
sampling,  which  does  not  happen  in  the  actual  device,  so  we  need  to  make  sure  it  doesn’t 
happen  in  the  model. 
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Figure  D.  1  Periodic  Waveform  Assumption  of  the  Discrete  Fourier  Transform 

To  eliminate  this,  a  windowing  filter  is  used  to  attenuate  signals  approaching  the 
sampling  edges  (note;  there  is  no  hard  edge  here,  just  the  limits  of  the  selected  sampling 
region).  This,  however,  must  be  done  with  care  because  the  waveform  must  remain  as 
accurate  as  possible  in  both  spatial  and  fi-equency  planes.  An  abrupt,  severe,  rectangular 
attenuation  preserves  the  spatial  characteristics  across  the  wave,  but  acts  like  a  mirror  at 
the  edges,  reflecting  the  wave  back  into  the  sample  space.  Additionally,  choosing  an 
attenuator  that  is  too  sharp  and  narrow  will  cause  a  Gibbs  Phenomenon,  or  ripples  in  the 
frequency  domain.  Figure  D.2  (b).  Choosing  a  smooth  slowly  varying  window,  such  as  a 
Hamming  function,  is  desirable  in  the  fi-equency  plane,  but  vwU  significantly  distort  the 
wave  in  the  spatial  domain.  Figure  D.2  (a). 

Taken  mathematically,  the  spatial  filter  is  a  multiplication  of  two  functions  in  the 
spatial  plane,  filter  t(x)  and  function  f(x).  Transforming  to  the  spatial  frequency  plane,  this 
becomes  the  convolution  of  the  two  separate  transforms. 
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g(x)  =  t(x)/(x) 
G(K)=m)*F(kJ 


Ideally,  in  the  space  domain,  we  would  like  a  rectangle  function,  however,  this  results 
in  a  sine  function  (which  has  ripples)  in  the  frequency  domain,  giving  rise  to  the  Gibbs 
Phenomena.  In  the  frequency  domain,  we  would  ideally  like  the  filter  transform,  T(kx),  to 
be  a  delta  fimetion,  however,  this  would  represent  no  fihering  at  all.  One  approximation 
for  the  delta  function  is  a  narrow  Gaussian.  This  has  been  shown  by  Harvey  (Am.  Journal 
of  Physics,  Aug  1981)  be  the  best  at  reduction  of  ripples,  at  the  expense  of  smoothing  any 
sharp  profiles  in  the  frequency  domain.  Unfortunately,  a  Gaussian  function  approximating 
a  delta  function  in  the  frequency  plane  may  not  provide  enough  attenuation  at  the  edges  of 
the  spatial  sampling  region.  Combining  these  to  functions  results  in  a  “best  fit”  solution. 
In  the  frequency  domain,  multiply  these  two  functions  together.  If  the  widths  are  properly 
selected,  most  of  the  ripples  in  the  frequency  domain  will  be  substantially  reduced.  The 
resuhs  in  the  spatial  domain  is  the  convolution  of  the  two  functions,  or  softening  the  edges 
of  the  rectangle  function  without  creating  any  ripple  in  the  spatial  domain.  Figure  D.2  (c). 
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Hanntng  Sp«cfciB(  Pterw 


Henning  Spettai  Plane 


Figure  D.2  Windowing  Functions,  Spatial  and  Frequency  Plane  Spectrums.  (a) 
Hanning  Filter,  (b)  Rectangle  Filter,  (c)  Convolution  Filter. 

The  distance  between  sampled  points  and  the  total  size  of  the  sampling  space  directly 
effects  the  accuracy  of  the  BPM  not  only  through  the  accuracy  of  the  waveform  itself,  but 
also  in  the  sampling  of  index  waveguides.  The  Nyquist  sampling  theorem  states  that  the 
largest  frequency  that  can  be  sampled  is  V2  of  the  sampling  rate.  Spatial  frequencies 
exceeding  this  will  alias  onto  the  existing  frequency  space  causing  distortion. 
Additionally,  the  resolution  of  the  frequency  space  or  the  distance  between  discrete 
frequency  points  will  be  determined  by  the  inverse  of  the  size  of  the  sampling  space. 
Using  these  two  criteria,  we  would  like  to  sample  as  many  points  as  possible!  Opposing 
the  desire  to  sample  le6  points  across  a  huge  sample  space,  however,  is  the  time  a 
computer  takes  to  calculate  the  DFT  for  the  sample  space  multiplied  by  the  number  of 
steps  taken  for  the  propagating  beam.  Experience  has  shown  that  even  1000  point 
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samplings  on  a  Pentium  90  processor  using  MATLAB,  propagated  over  1000  steps  will 
take  about  10  minutes  to  complete  the  calculation. 

The  size  of  Az  in  the  calculations  also  effects  the  accuracy  and  speed  of  the  calculation 
because  the  method  assumes  that  the  index  profile  is  constant  over  the  distance  Az. 
VanRoey,  (“Beam-propagation  method;  analysis  and  assessment,”  J.  Opt.  Soc.  Am.,  July 
1981),  derives  the  following  restrictions  for  the  use  of  the  BPM. 


|AnU  sin^(a) « 1 
p  »  —yll 

«o 

p»—^i27rXM^^^z 

^0 

Az.j^tan(a)|Aw|„„  «1 
tan(a)  ^ 

n - Az  « 1 


where 

A  is  the  free  space  wavelength 

a  h  the  maximum  angle  corresponding  to  a  significant  part  of  the  angular 
spectrum,  usually  the  Numeric  .^erture  angle  of  the  waveguide 
(N.A.=sin(«)) 

p  is  the  highest  spatial-frequency  component  of  the  refractive-index  profile  An\x) 


VanRoey  illustrated  the  effects  of  varjdng  step  size  by  propagating  a  2.5  pm  wide  (1/e 
width)  Gaussian  beam  (A=l  pm)  through  10  pm  wide  parabolic-index  profile  waveguide 
with  a  peak  refractive-index  variation  of  1%.  He  then  plotted  the  quadratic  error  of  the 
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BPM  propagated  field  compared  to  the  closed  form  solution.  The  results  are  shown  in 
Figure  D.3. 


Figure  D.3  Percent  Error  versus  Propagation  Distance  for  Varying  Step  Sizes  (J.  Opt 
Soc.  Am.A^ol.  71,  No.  7/July  1981  pg.  806) 
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APPENDIX  E:  PROCESS  FOLLOWER  -  SOLITON  DEVICE 


In  this  process,  the  thinning  of  the  sample  to  200  microns  was  done  first.  This  allowed 
the  high  temperature  Ni/Ge/Au  ohmic  contact  to  be  applied  and  annealed  without 
damaging  the  lower  temperature  Au/Zu/Au  topside  ohmic  contact.  Additionally,  the  SiN 
used  to  insulate  the  topside  remains  intact  in  tests  up  to  350“C,  however  it  occasionally 
delaminates  at  450®C,  the  processing  temperature  of  Ni/Ge/Au  contacts.  Thus,  the  SiN  is 
not  exposed  to  the  420®C  Ni/Ge/Au  anneal  temperature.  The  major  disadvantage  of  this 
method  is  that  the  200-micron  sample  is  VERY  fi'agile  and  must  be  handled  with  extreme 
care. 

Thinning  after  the  topside  processing  can  be  done,  however,  the  silicon  nitride  should 
be  checked  carefully  after  the  420“C,  120  sec  rapid  thermal  anneal. 


Initial  Cleave 

1 .  Prior  to  cleaving,  use  a  scribe  to  lightly  maik  the  backside  of  the  wafer  where  the  sample  will  be  cut 
from.  After  the  cleave,  the  top  and  bottom  of  the  sample  will  be  indistinguishable  without  this  mark. 

2.  Cleave  a  section  of  the  wafer  by  making  a  1-mm  scribe  mark  at  the  edge  of  the  wafer,  turning  the 
wafer  over,  and  applying  pressure  to  the  backside  of  the  wafer  immediately  above  the  scribe  mark 
with  a  pointed  object.  The  edge  of  the  subsequent  cleaved  sample  must  be  smooth  along  most  of  one 
side.  TTiis  will  be  used  to  align  the  mask.  The  sample  should  te  no  smaller  than  3/8”  square  - 
remember  that  the  edges  will  not  be  used. 


Thinning 
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3.  Measure  the  starting  thickness  of  your  sample  using  a  micrometer  caliper  niler.  G2-2329is625 
microns  thick.  This  indicates  how  much  material  needs  to  be  removed. 

4.  Mount  the  sample  onto  a  thin  glass  coverlet  using  low  temperature,  acetone  soluble  wax  such  as 
crystal  bond. 

5.  Moxmt  the  sample/coverlet  combination  onto  the  top  surface  of  a  lapping  chuck  using  crystal  bond 

6.  Measure  the  chuck/sample/coverlet  combination.  Zero  out  a  micrometer  at  this  thickness.  Use  this 
zeroed  micrometer  as  a  guide  to  measure  material  removed. 

7.  Mount  the  lapping  chuck  into  a  lapping  jig.  If  the  jig  has  a  thickness  setting  available,  set  this  to 
remove  the  desired  amoimt  of  material. 

8.  Mix  a  slurry  of  2  tablespoons  of  10  micron  alumna  lapping  grit  and  water  on  the  surfece  of  a  lapping 
plate.  The  slimy  should  be  the  consistency  of  mud.  Spread  the  slurry  into  a  2”  circle  on  the  plate. 

9.  Place  the  lapping  jig  into  the  slurry.  If  using  the  lapping  machine,  set  or  a  slow  table  and  chuck 
rotation.  Alxmt  1.5  on  both  dials. 

10.  If  lapping  by  hand,  use  figure  8  motions,  rotating  the  chuck  by  120“  every  100  strokes,  to  remove  the 
desir^  amount  of  material. 

11.  Once  the  sample  is  thinned,  clean  with  a  toothbrush  to  remove  arty  10  micron  grit. 

12.  Using  a  polishing  pad,  mix  a  slurry  of  1  or  3  micron  polishing  grit  with  a  mixture  of  1  part  Clorox 
bleach  to  3  parts  water.  This  is  an  aggressive  chemical  polish  that  will  stain  the  back  of  the  sample  if 
not  removed  quickly  after  polishing. 

13.  Using  the  figure  8  motion,  polish  the  sample  until  aO  scratches  are  removed  from  the  surface.  To 
clean  the  slurry,  have  the  water  running,  and  quickly  move  the  sample  from  a  polishing  motion  on 
the  pad  to  the  water.  Ifthe  slurry  does  not  come  oflFinunediately,  a  stain  will  result.  If  this  happens, 
add  a  small  amount  of  water  (few  drops)  to  the  slurry,  and  polish  again  for  100  strokes.  Move  from  a 
polishing  motion  to  the  water  and  quickly  remove  the  slurry.  A  completed  sample  will  not  have  any 
visible  scratches  and  be  mirror  shirmy.  A  mottled  orange  peal  effect  nuty  be  visible.  This  is  not 
detrimental  to  the  device  and  may  be  ignored. 

14.  Remove  the  laj^g  chuck  firom  the  lai5)ing  jig  and  place  in  a  beaker  of  acetone.  After  several  hours, 
the  coverlet  and  sample  will  separate  from  the  lapping  chuck. 

15.  MOVE  TO  THE  CLEANROOM  FOR  FURTHER  PROCESSING.  Set  hot  plates  for  1 19“C  and 
200“C. 

16.  In  the  clean  room,  remove  the  sample  from  the  acetone  bath,  and  spray  with  a  stream  of  acetone  to 
remove  arty  residual  wax  from  the  solution.  Do  not  let  the  sample  dry. 

17.  Place  the  gample  into  a  clean  beaker  of  acetone,  and  ultrasound  the  beaker  for  5  minutes  using  the 
Cole-Parmer  ultrasound  in  the  class  10,000  area.  DO  NOT  use  the  small  Fisher  ultrasound  in  the 
class  100  area.  This  ultrasound  will  break  up  thirmed  wafers;  it’ll  break  up  un-thiimed  wafers  too! 

18.  Spray  with  acetone  for  several  seconds  as  you  remove  the  sample  from  the  acetone  bath,  then  place 
the  sample  into  another  clean  beaker  filled  with  acetone.  DO  NOT  allow  sample  to  dry  during  the 
solvent  cleaning  process. 

19.  Boil  this  beaker  of  acetone  at  1 19°C  for  5  minutes. 

20.  Remove  from  acetone  bath  and  spray  with  acetone  for  several  seconds.  Spray  with  methanol  to 
remove  the  acetone  for  several  seconds. 
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21.  Place  in  a  clean  beaker  fill  with  methanol.  Boil  sample  in  methanol  at  IW^C  for  5  minutes. 

22.  Remove  from  methanol;  do  not  allow  sample  to  dry.  Spray  with  methanol  for  several  seconds  as  you 
remove  the  sample  from  the  methanol  bath,  then  spray  with  isopropyl  alcohol  (IPA)  for  several 
seconds  to  remove  the  methanol. 

23.  Boil  sample  in  IPA  at  200®C  for  5  minutes. 

24.  Remove  from  IPA  do  not  allow  sample  to  dry.  Spray  with  IPA  for  several  seconds  as  you  remove  the 
sample  from  the  bath. 

25.  While  still  covered  with  IPA  place  on  a  lint  free  rag  and  blow  the  IPA  from  the  sample  with  nitrogen 
gun.  Do  not  let  the  IPA  evaporate  without  the  nitrogen  blow. 


N-tvne  Metallization 

This  step  employs  a  ine-metal  dip  to  remove  any  sur&ce  oxicfes  built  up  upon  the  wafer  during  the 
thinning  process.  The  sample  should  be  loaded  into  the  metallizer  as  soon  as  possible  after  the  pre-metal 
dip. 

Prep 

Mix  a  pre-metal  dip  solution  of  1  part  buffered  oxide  etch  to  50  parts  18  megohm  deionized  water 
from  the  class  10,000  area. 

26.  Use  a  drop  of  photoresist  to  attach  the  sample,  backside  up,  to  a  glass  slide. 

27.  Hot  plate  bake  the  sample  at  1 19®C  for  5  minutes  to  cure  the  resist. 

28.  Dip  the  sample  into  the  pre-metal  dip  solution  for  30  seconds. 

29.  DIW  rinse  for  2  minutes. 

30.  Blow  dry  thoroughly. 

31.  Blow  dry  some  more. 

32.  Load  into  the  metallizer  and  press  ‘cycle’. 

33.  Pump  the  metallizer  to  IxlO"®  torr  and  use  standard  posted  techniques  to  deposit  50A  Ni/  200A  Ge/ 
2000AAu. 

34.  Cool  the  metallizer  30  minutes  prior  to  removing  the  sample. 

35.  Remove  the  sample,  and  soak  in  a  beaker  of  acetone  until  the  sample  separates  ftom  the  glass  slide. 

36.  While  spraying  the  sample  with  acetone,  remove  the  sample  from  the  soaking  acetone,  and  place  in  a 
clean  beaker  of  acetone.  Boil  sample  in  acetone  at  1 19®C  for  5  minutes. 

37.  Remove  fi-om  acetone,  do  not  allow  sample  to  dry.  Spray  with  acetone  for  several  seconds  as  you 
remove  the  sample  from  the  acetone  bath,  then  spray  with  methanol  to  remove  the  acetone  for  several 
seconds.  DO  NOT  allow  sample  to  become  dry  Airing  the  solvent  cleaning  process. 

38.  Boil  sample  in  methanol  at  150°C  for  5  minutes. 

39.  Remove  from  methanol,  do  not  allow  sample  to  dry.  Spray  with  methanol  for  several  seconds  as  you 
remove  the  sample  from  the  methanol  bath,  then  spr^  with  isopropyl  alcohol  (IPA)  for  several 
seconds  to  remove  the  methanol. 
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40.  Boil  sample  in  IPA  at  200®C  for  5  minutes. 

41 .  Remove  from  IPA,  do  not  allow  sample  to  dry.  Spray  with  IPA  for  several  seconds  as  you  remove  the 
sample  from  the  bath. 

42.  While  still  covered  with  IPA,  place  on  a  lint  free  rag  and  blow  the  IPA  from  the  sample  with  nitrogen 
gun.  Do  not  let  the  IPA  evaporate  without  the  nitrogen  blow. 

N-contact  Rapid  Thermal  Anneal 

43.  Place  a  clean  GaAs  mechanical  wafer  onto  the  silicon  substrate  of  the  Rapid  Thermal  Annealer 
(RTA).  Place  the  sample,  topside  down,  onto  this  mechanical  substrate. 

44.  Close  the  RTA  chamber  and  flow  nitrogen  gas  at  the  maximum  rate  for  5  minutes. 

45.  Reduce  gas  flow  to  3  seem  wait  1  minute  for  the  chamber  to  stabilize. 

46.  RTA  the  sample  at  420'’C  for  120  seconds. 

47.  Allow  the  RTA  to  cool  to  below  200"C  before  opening. 


This  completes  backside  processing 


If  you  wish  to  process  the  front  side  first,  start  here. 


Clean  and  Kiotoresist 

Cleaning  must  immediately  proceed  photoresist,  exposures,  and  develops.  If  you  cannot  do  all  three 
in  a  single  day,  consider  waiting  until  you  can. 

Prep 

Ensiue  exposure  tool  nitrogen  tank  in  chase  contains  at  least  750  lbs.  of  pressure  and  is  on. 

Ensure  house  nitrogen  tank  in  chase  is  on. 

Ensure  house  vacuum  in  chase  is  on. 

Start  with  one  hot  plate  at  1 19“C  and  the  second  at  200“C. 

Turn  on  exposure  tool  gas,  and  start  exposure  tool  lamps. 

Set  the  exposure  time  of  aligner  to  1.6  seconds. 


48.  Boil  sample  in  acetone  at  1 19“C  for  5  minutes. 

49.  Remove  from  acetone,  do  not  allow  sample  to  dry.  Spray  with  acetone  for  several  seconds  as  you 
remove  the  sample  from  the  acetone  bath,  then  spray  with  methanol  to  remove  the  acetone  for  several 
seconds.  DO  NOT  allow  sample  to  become  dry  during  the  solvent  cleaning  process. 

50.  Boil  sample  in  methanol  at  150®C  for  5  minutes. 

5 1 .  Remove  from  methanol,  do  not  allow  sample  to  dry.  Spray  with  methanol  for  several  seconds  as  you 
remove  the  sample  from  the  methanol  bath,  then  spray  with  isopropyl  alcohol  GPA)  for  several 
seconds  to  remove  the  methanol. 
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52.  Boil  sample  in  IPA  at  200®C  for  5  minutes. 

53.  Remove  fiom  IPA,  do  not  allow  sample  to  diy.  Spray  with  IPA  for  several  seconds  as  you  remove  the 
sample  from  the  bath. 

54.  While  still  covered  with  IPA,  place  on  a  lint  free  rag  and  blow  the  IPA  fiom  the  sample  with  nitrogen 
gun.  Do  not  let  the  IPA  evaporate  without  the  nitrogen  blow. 

55.  Place  on  a  clean  glass  microscope  slide,  and  hot  plate  bake  at  200“C  for  5  minutes. 

56.  Remove  fiom  hot  plate  and  allow  sample  to  cool  for  1  minute. 

57.  Center  sample  on  spinner  chuck,  engage  vacuum.  Check  vacuum  hold. 

58.  Test  spin  sample  to  ensure  acceleration,  spin  speed,  and  spin  time.  Sample  should  spin  iq)  to  speed  in 
1  second,  spin  at  3000  rpm,  and  spin  for  30  seconds. 

59.  With  sample  stoiq)ed,  dispense  enough  1400-17  photoresist  to  cover  sample,  typically  1  to  3  drops, 
but  not  bleed  over  edges.  Immediately  upon  completing  the  dispense,  start  spinner. 

60.  Remove  fiom  spinner  without  disturbing  top  of  resist.  This  can  be  tricky  so  be  careful.  Ifthe  resist  is 
not  completely  covered  and  free  from  dust  eddies  and  holes,  briefly  rinse  sample  with  acetone  and  go 
to  step  2. 

61.  Place  sample  face  up  on  a  clean  microscope  slide  and  hot  plate  bake  for  5  minutes  at  119*’C. 

62.  Remove  from  hot  plate  and  cool  for  5  minutes. 


Exposure 

Prep 

Ensure  lamp  control  power  supply  has  a  number  close  to  274  watts  on  the  display. 

If  the  display  says  cold,  wait  until  a  number  is  display. 

Turn  the  MJB-3  exposure  tool  on  by  pressing  the  power  button 

Set  the  MJB-3  to  HP  (High  Precision)  mode  by  pressing  the  HP/ST  button  until  the  HP  side  lights 
Ensure  the  vacuum  mask  light  is  off 

Clean  the  sample  holder  chuck  with  small  amount  of  acetone  on  a  wipe  if  necessary. 

Mix  a  solution  of  1  part  351  developer  to  6  parts  water.  Typically  25ml  / 150  ml. 

Find  the  mask  you  will  use  and  visually  inspect  to  ensure  mask  is  clean.  If  not,  clean  it. 


63.  With  nitrogen  gun,  blow  surface  of  exposure  tool  to  remove  dust.  Do  not  attach  mask  yet. 

64.  If  small  amoimts  of  photoresist  are  on  the  back  of  the  sample,  place  sample  upside  down  on  the 
sample  chuck  of  the  exposure  tool,  and  expose  for  99  seconds.  If  large  amoimts  of  resist  are  on 
sample,  remove  with  acetone  being  careful  not  to  get  acetone  on  the  front  of  the  sample.  If  this 
happens,  however,  rinse  briefly  with  acetone  and  go  to  step  2. 

65.  This  step  is  used  to  remove  edge  beading.  Turn  sample  over  and  place  a  clean  mechanical  sample 
over  the  sample  such  that  1  mm  of  two  edges  are  exposed.  Expose  for  60  seconds.  Repeat  procedure 
to  exposure  the  two  remaining  edges. 
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66.  Develop  sample  in  developer  for  60  seconds  or  until  all  back  and  edge  photoresist  is  removed.  Use 
caution  not  to  disturb  photoresist  on  the  front  of  the  sample. 

67.  Rinse  in  Deionized  Tap  Water  piW)  for  2  minutes.  Thoroughly  blow  dry  with  nitrogen  gun.  A  hot 
plate  bake  is  not  necessary  at  this  point. 

68.  Attach  mask  to  the  mask  chuck  and  press  vacuum  mask.  This  first  exposure  requires  the  bright  field 
WAVEGUIDE  quadrant  of  the  mask.  The  dull  chromium  side  of  the  mask  must  be  facing  away  from  the 
holder.  Insert  holder  into  the  aligner  and  tighten  thumbscrews. 

69.  Place  sample  face  up  on  MJB-3  sample  chuck.  Sample  must  cover  the  flush  open  holes  on  the  chuck. 
If  not,  place  small  squares  of  tape  over  open  holes.  Note  that  the  large  depressed  hole  on  the  left  of 
the  chuck  must  not  be  covered  or  blocked.  This  is  the  nitrogen  hole  for  separation  mode. 

70.  While  watching  sample  fix>m  the  side,  or  through  the  microscope,  raise  the  sample  to  the  mask.  The 
sample  should  touch  the  mask  when  the  substrate  raiser  lever  is  at  the  90%  position.  If  the  sample 
doesn’t  contact  at  the  90%  position,  raise  or  lower  the  substrate  chuck  with  the  chuck  dial  gauge.  Be 
careful;  you  can  destroy  samples  and  damage  your  mask  if  the  chuck  is  set  too  high;  start  low  to  be 
safe. 

If  the  sample  jumps  to  the  mask,  and  cannot  be  put  into  separation,  one  or  more  of  the  holes  in 
the  substrate  chuck  are  not  covered  properly.  Lower  chuck  and  check  for  good  coverage. 

If  the  sample  rotates  in  separation  mode,  the  chuck  is  set  just  a  bit  too  hi^  and  one  comer  is 
catching.  Lower  chuck. 

If  in  contact  mode,  Fresnel  contact  rings  to  not  form,  the  chuck  is  too  low.  Raise  chuck  holder 
until  rings  form  at  the  comers,  but  go  away  in  separation  mode. 


71.  Align  sample  cleaved  edge  to  a  straight  edge  on  the  mask. 

72.  Press  the  vacuum  chamber  button  and  ensure  vacuum  pressure  drops  to  -0.  8  bar  on  the  gauge.  If  not, 
ensure  vacuum  line  to  the  UV  exposure  tower  is  off.  Next  check  to  ensure  vacuum  holes  on  the 
sample  chuck  are  covered,  and  finally,  check  to  ensure  vacuum  ring  on  the  sample  chuck  is  good  and 
touches  the  mask  around  its  entire  circumference  to  the  mask  during  contact. 

73.  If  pressure  drops  to  a  good  level,  visually  check  the  alignment  once  again.  Note  that  vacuum 
chamber  presses  the  sample  closer  to  the  mask,  and  sometimes  moves  the  sample. 

74.  If  alignment  is  ok,  ensure  exposure  time  is  set  for  1.6  seconds,  and  press  the  expose  button. 

75.  Remove  sample  from  e?qx>sure  tool,  and  place  in  developer  for  30  to  45  seconds.  You  can  see  the 
resist  develop  away.  If  watching  the  sample  in  the  developer,  hold  sample  in  developer  10  seconds 
past  the  point  where  exposed  photoresist  appears  to  be  removed. 

76.  Rinse  in  DIW  for  2  minutes,  then  blow  dry. 

77.  Examine  with  optical  microscope  with  the  green  filter  in  place.  If  resist  is  not  completely  developed, 
develop  again  for  15-second  increments  followed  1^  a  2  minutes  DIW  rinse  and  blow  dry.  Sometimes 
this  can  compensate  for  a  bad  exposure,  however,  be  picky.  If  the  exposure  isn’t  perfect,  rinse  briefly 
with  acetone  and  go  back  to  step  2. 

78.  Place  sample  on  a  clean  microscope  slide  and  hot  plate  bake  for  5  minutes  at  1 19®C.  This  is  critical 
to  ensure  the  resist  does  not  delaminate  during  subsequent  processing. 


Etch  Waveguides 
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Note;  The  Silicon  Nitride  process  immediately  following  the  etch  process  requires  a  non-oxidized 
surface  for  good  adhesion.  As  exposed  GaAs  and  AlGaAs  will  oxide  slowly  at  room  temperature,  you 
should  be  prepared  to  load  the  sample  into  the  deposition  chamber  and  pump  the  chamber  to  vacuum 
immediately  following  the  etch  process. 

Prep 

Mix  a  solution  of  1:2:50  of  H3PO4:  H202:  DIW,  typically  5ml,  10ml,  250ml.  For  this,  vse  the  18 
Meg  water  from  the  Millipore  system  in  the  main  cleanroom.  Clean  glassware  thorou^y  using  18 
Meg  DIW  prior  to  mixing.  Stir  solution  until  viscosity  turbulence  is  no  longer  visible.  Solution  must 
sit  for  at  least  5  minutes  prior  to  use. 

The  etch  rate  of  this  solution  is  {q^Moximately  1850A/min  at  23®C.  Check  temperature  of 
solution.  Ifnot23®C,  adjust  times  accordingly.  You  can  repeat  etches  to  etch  deeper,  the  reverse  is 
not  true.  Err  to  the  short  side! 

Ensure  SiN  Deposition  chamber  on  the  PlasmaTherm  790  at  or  below  50®C. 

Set  the  PlasmaTherm  controller  to  the  right  side  deposition  chamber. 

Load  the  SINIOO  process  and  prepare  to  run  this  process. 

Set  hotplates  at  1 19®C  and  200®C 


79.  Remove  residual  organic  surface  contamination  1^  placing  the  sample  into  the  left  chamber  of  the 
PlasmaTherm  790  and  executing  the  descum  procedure.  Run  the  procedure  for  1  minute. 

80.  Remove  sample  from  PlasmaTherm  790.  Measure  the  resist  thickness  at  several  positions  on  the 
alpha  step  profilometer.  Select  positions  that  you  can  return  to  in  subsequent  measurements.  Record 
resist  thickness  for  later  measurements.  Center  thickness  for  1400-17  resist  at  3000rpm  spin  speed 
should  be  5000 A  ±  500 A  Thickness  will  vary  towards  edges. 

8 1 .  Prior  to  etch,  stir  etch  solution  to  ensure  viscosity  turbulence  is  not  visible.  Immerse  sample  free  up 
or  sideways  into  acid  etch  solution  ft)r  appropriate  time.  Never  allow  the  face  of  the  sample  to  reset 
against  the  walls  or  bottom  of  etch  beaker.  Assume  1850A/min.  Agitate  gently  each  minute  during 
etch. 

82.  Remove  sample  from  etch  solution  and  plunge  into  beaker  of  clean  DIW.  Rinse  for  2  minutes  in 
DIW,  then  blow  dry. 

83.  Measure  etch  depth  +  resist  thickness  at  the  same  positions  you  measured  the  initial  resist  thickness 
and  subtract  resist  thickness  to  determine  actual  etch  depth.  Assume  that  the  resist  will  swell  slightly 
2OO-500A  from  pre-etch  measurement.  Calculate  your  etch  rate.  Use  this  rate  for  subsequent  etches. 

84.  If  the  depth  is  short,  go  to  step  35  using  your  calculated  rate  from  step  37  to  set  the  next  time.  If  the 
depth  is  correct,  remove  resist  with  a  short  acetone,  methanol,  IPA  rinse,  and  then  blow  dry. 

Measure  and  record  final  thickness. 

85.  Examine  under  a  microscope.  It  is  common  for  the  etched  area  to  look  rough  and  nonuniform 
especially  when  viewed  with  a  Newmarski-Differential  Interference  Contrast  (DIC)  system.  The 
areas  not  etched  should  have  sharp  and  straight  edges  and  be  smooth.  Check  the  sharpness  ^d 
continuity  ofthe  etch  on  the  waveguides.  Ensure  separation  between  waveguide  breaks.  Rejection 
criteria  at  this  point  is  up  to  the  fabricator,  however,  this  is  still  early  enough  in  the  process  to  restart 
without  significant  loss  of  time  or  resources. 
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Silicon  Nitride  Deposition 

86.  Boil  sample  in  acetone  at  1 19°C  for  5  minutes. 

87.  Remove  from  acetone,  do  not  allow  sample  to  dry.  Spray  with  acetone  for  several  seconds  as  you 
remove  the  sample  from  the  acetone  bath,  then  spray  with  methanol  to  remove  the  acetone  for  several 
seconds.  DO  NOT  allow  sample  to  become  drv  during  the  solvent  cleaning  process. 

88.  Boil  sample  in  methanol  at  150®C  for  5  minutes. 

89.  Remove  from  methanol,  do  not  allow  sample  to  dry.  Spray  with  methanol  for  several  seconds  as  you 
remove  the  sample  from  the  methanol  bath,  then  spray  with  isopropyl  alcohol  (IPA)  to  remove  the 
methanol  for  several  seconds. 

90.  Boil  sample  in  IPA  at  200®C  for  5  minutes. 

91.  Remove  from  IPA,  do  not  allow  sample  to  dry.  Spray  with  BPA  for  several  seconds  as  you  remove  the 
sample  from  the  bath. 

92.  While  still  covered  with  IPA,  place  on  a  lint  free  rag  and  blow  the  IPA  from  the  sample  with  nitrogen 
gun.  Do  not  let  the  IPA  evaporate  without  the  nitrogen  blow. 

93 .  Do  not  hot  plate  bake.  This  will  cause  oxides  to  form  on  the  surface  of  the  sample. 

94.  Immediately  load  into  the  right  chamber  of  the  PlasmaTherm  790.  Close  chamber  and  manually 
pump  with  tuibopump. 

95.  When  the  chamber  reaches  5x10"’  torr,  manually  flow  50sccm  Na  and  50  seem  He  at  a  chamber 
pressiue  of  100  millitorr  for  2  minutes.  This  is  a  flush  step  to  help  remove  oxygen. 

96.  Stop  manual  mode  and  pump  with  the  turbo  pump  imtil  the  chamber  pressure  reduces  to  1x10"*  torr. 

97.  Set  the  standby  temperature  of  the  ^stem  to  250®C  and  wait  until  the  platen  temperature  is  within 
5"Cof250“C. 

98.  Start  the  SINIOO  process.  Run  this  process  to  6  minutes.  This  will  deposit  approximate  lOOOA  of 
silicon  nitride. 

99.  When  the  process  completes,  set  the  q^stem  standby  temperature  to  50®C  and  allow  the  sanqrle  to  cool 
under  vacuum  until  the  platen  temperature  is  below  100®C.  Note  that  the  system  temperature  may  be 
below  50°  at  the  end  of  Ae  SINIOO  process. 

100. Remove  sample  from  PlasmaTherm  790  and  reseal  and  tuibopump  system. 

101.  Examine  sample  under  optical  microscope.  The  SiN  film  should  follow  the  contom  of  the  substrate 
and  should  be  smooth  without  any  major  holes  or  cracks.  Small-discolored  areas  are  typical  and  have 
not  proven  to  be  detrimental  to  the  processing.  If  any  delaminating  occurs,  the  film  is  defective. 

102.1f  the  film  needs  to  be  removed,  place  the  sample  into  the  deposition  chaniber  and  run  the  etch 
program  for  8  minutes.  Note,  this  process  will  be  updated  with  a  new  process  to  be  used  in  the  left 
etch  chamber  during  the  month  of  June.  Removing  the  film  will  not  damage  the  sample.  Once  the 
film  is  removed,  go  to  step  40. 


P  type  Metallization  -  Photoresist 

Note;  Following  SiN  deposition,  sample  should  be  very  clean  and  does  not  require  additional  cleaning 
at  this  time.  If  just  removed  fiom  the  PlasmaTherm  790,  photoresist  can  be  immediately  applied.  If  the 
sanqile  is  left  out  overnight,  hot  plate  bake  at  1 19°C  for  5  minutes  to  drive  out  any  water  collected. 
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The  photoresist  used  in  this  step  acts  as  a  photomask  for  both  the  opening  of  the  silicon  nitride  and 
the  deposition  of  the  p-type  contact.  This  type  ofdual  use  is  termed  a  self-aligned  process.  Remember  not 
to  remove  the  photoresist  after  the  silicon  nitride  is  removed. 

Prep 

Hot  plate  at  119®C 

Check  availability  of  1813  photoresist.  If  out,  mix  using  28  ml  of  1827  resist  with  10  ml  type  P 
thinner.  Mix  rolling  the  bottle  on  the  table  and  gently  rotating  end  to  end  for  5  minutes.  DO  NOT 
SHAKE  or  ultrasound.  Shaking  will  get  bubbles  into  the  resist;  ultrasound  will  break  up  the  long 
chain  photoresist  polymers. 

Ensure  spinner  rotates  at  4000  rpm. 


103. Center  sample  on  spinner  chuck,  engage  vacuiun.  Check  vacuum  hold. 

104.  Test  spin  sample  to  ensure  acceleration,  spin  speed,  and  spin  time.  Sample  should  spin  up  to  speed  in 
1  second,  spin  at  4000  rpm,  and  spin  for  30  seconds. 

105. With  sample  stopped,  dispense  enough  1813  photoresist  to  cover  sample,  ^ically  1  to  3  drops,  but 
not  bleed  over  edges.  Immediately  upon  completing  the  dispense,  start  spinner. 

106. Remove  firom  spinner  without  disturbing  top  of  resist.  This  can  be  tricky  so  be  careful.  Ifthe  resist  is 
not  completely  covered  and  free  from  dust  eddies  and  holes,  briefly  rinse  sample  with  acetone,  and 
clean  using  the  boiling  solvent  technique.  Hot  plate  bake  for  5  minutes  at  1 19®C,  allow  to  cool  for  1 
minute  and  go  to  step  57. 

107. Place  sample  face  up  on  a  clean  microscope  slide  and  hot  plate  bake  for  5  minutes  at  1 19°C. 

108. Remove  from  hot  plate  and  cool  for  5  minutes. 


P-Tvpe  Metallization  -  Exposure 
Prep 

Ensure  lamp  control  power  supply  has  a  number  close  to  274  watts  on  the  display. 

If  the  display  says  cold,  wait  imtil  a  number  is  display. 

Turn  the  MJB-3  exposure  tool  on  by  pressing  the  power  button 

Set  the  MJB-3  to  HP  (High  Precision)  mode  by  pressing  the  HP/ST  button  until  the  HP  side  lights 
Ensure  the  vacuum  mask  light  is  off 

Clean  the  sample  holder  chuck  with  small  amount  of  acetone  on  a  wipe  if  necessary. 

Mix  a  solution  of  1  part  351  developer  to  6  parts  water.  Typically  25ml  / 150  ml. 

Find  the  mask  you  will  use  and  visually  inspect  to  ensure  made  is  clean.  If  not,  clean  it. 

109.  With  nitrogen  gun,  blow  surface  of  exposure  tool  to  remove  dust.  Do  not  attach  mask  yet. 

1  lO.If  small  amounts  of  photoresist  are  on  the  back  of  the  sample,  place  sample  upside  down  on  the 
sample  chuck  of  the  exposure  tool,  and  ejqwse  for  99  seconds.  If  large  amounts  of  resist  are  on 
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sample,  remove  with  acetone  being  careful  not  to  get  acetone  on  the  front  of  the  sample.  If  this 
happens,  however,  rinse  briefly  with  acetone  and  go  to  step  2. 

1 1 1  .This  step  is  used  to  remove  edge  beading.  Turn  sample  over  and  place  a  clean  mechanical  sample 
over  the  sample  such  that  1  mm  of  two  edges  is  exposed.  Expose  for  60  seconds.  Repeat  procedure  to 
exposure  the  two  remaining  edges. 

1 12. Develop  sample  in  developer  for  60  seconds  or  until  all  back  and  edge  photoresist  is  removed.  Use 
caution  not  to  disturb  photoresist  on  the  front  of  the  sample. 

1 13.  Rinse  in  Deionized  Tap  Water  (DIW)  for  2  minutes.  Thoroughly  blow-dry  with  nitrogen  gun.  A  hot 
plate  bake  is  not  necessary  at  this  point. 

1 14.  Attach  mask  to  the  mask  chuck  and  press  vacuum  mask.  This  first  exposure  requires  the  dark  field 
WAVEGUIDE  quadrant  of  the  mask.  The  dull  chromium  side  of  the  mask  must  be  facing  aw^  from  the 
holder.  Insert  holder  into  the  aligner  and  tighten  thumbscrews. 

1 15. Place  sample  free  up  on  MJB-3  sample  chuck.  Sample  must  cover  the  flush  open  holes  on  the  chuck. 
If  not,  place  small  squares  of  tape  over  open  holes.  Note  that  the  large  depressed  hole  on  the  left  of 
the  chuck  must  not  be  covered  or  blocked.  This  is  the  nitrogen  hole  for  separation  mode. 

1 16.  While  watching  sample  from  the  side,  or  through  the  microscope,  raise  the  sample  to  the  maslf  The 
sample  should  touch  the  mask  when  the  substrate  raiser  lever  is  at  the  90%  position.  If  the  sample 
doesn’t  contact  at  the  90%  position,  raise  or  lower  the  substrate  chuck  with  the  chuck  dial  gauge.  Be 
careful;  you  can  destroy  samples  and  damage  your  mask  if  the  chuck  is  set  too  high;  start  low  to  be 
safe. 


If  the  sample  jumps  to  the  mask,  and  carmot  be  put  into  separation,  one  or  more  of  the  holes  in 
the  substrate  chuck  are  not  covered  jKoperly.  Lower  chuck  and  check  for  good  coverage. 

If  the  sample  rotates  in  separation  mode,  the  chuck  is  set  just  a  bit  too  high  and  one  comer  is 
catching.  Lower  chuck. 

If  in  contact  mode,  Fresnel  contact  rings  to  not  form,  the  chuck  is  too  low.  Raise  chuck  holder 
until  rings  form  at  the  comers,  but  go  away  in  separation  mode. 


1 17.  Align  sample  cleaved  edge  to  a  straight  edge  on  the  mask. 

1  IS.Press  the  vacuum  chamber  button  and  ensure  vacuum  pressure  drops  to  -0.  8  bar  on  the  gauge.  If  not, 
ensure  vacuum  line  to  the  UV  exposure  tower  is  off.  Next  check  to  ensure  vacuum  holes  on  the 
sample  chuck  are  covered,  and  finally,  check  to  ensure  vacuum  ring  on  the  sample  chuck  is  good  and 
touches  the  mask  around  its  entire  circumference  to  the  mask  during  contact. 

1 19.1f  pressure  drops  to  a  good  level,  visually  check  the  alignment  once  again.  Note  that  vacuum 
chamber  presses  the  sample  closer  to  the  mask,  and  sometimes  moves  the  sample. 

120.1f  alignment  is  ok,  ensure  exposure  time  is  set  for  5.0  seconds,  and  press  the  expose  button. 

121.  Remove  sample  from  ejqwsure  tool,  and  place  in  developer  for  30  to  45  seconds.  You  can  see  the 
resist  develop  away.  If  watching  the  sample  in  the  developer,  hold  sample  in  developer  10  seconds 
past  the  point  where  exposed  photoresist  appears  to  be  removed. 

122. Rinse  in  DIW  for  2  minutes,  then  blow  dry. 

123.  Examine  with  optical  microscope  with  the  green  filter  in  place.  If  resist  is  not  completely  developed, 
develop  again  for  15-second  increments  followed  by  a  2  minutes  DIW  rinse  and  blow  dry.  Sometimes 
this  can  compensate  for  a  bad  exposure,  however,  be  picky.  If  the  exposure  isn’t  perfect,  briefly  rinse 
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sample  with  acetone,  and  clean  using  the  boiling  solvent  technique.  Hot  plate  bake  for  5  minutes  at 
1 19“C,  allow  to  cool  for  1  minute  and  go  to  step  57. 

124.Place  sample  on  a  clean  microscope  slide  and  hot  plate  bake  for  5  minutes  at  1 19®C.  This  is  critical 
to  ensure  ^e  resist  does  not  delaminate  during  subsequent  processing. 


P-Tvpe  Mctalli/ation  -  Nitride  Etch 
E*rep 

Ensure  the  temperature  of  the  PlasmaTherm  790  is  SOT  or  lower 
Load  process  etch  and  prepare  to  run  the  process 

125  .Place  sample  into  the  right  chamber  of  the  PlasmaTherm  790.  The  temperature  of  the  platen  must  be 
lower  than  lOOT.  Calibration  runs  of  this  process  were  done  at  SOT. 

126. Run  the  process  etch  for  3  minutes  to  etch  through  and  slightly  undercut  2000A  of  SiN. 

127. D0  not  remove  the  photoresist  for  the  subsequent  metallization. 

ns.Examine  the  sample  using  an  optical  microscope  following  this  step.  The  tops  of  the  waveguides 
must  be  ftee  fiom  nitride.  Remaining  nitride  looks  like  a  smudge  the  follows  the  edges  of  the 
openings.  If  required,  set  the  sample  for  a  1  minutes  nitride  etch  using  the  etch  process. 

P-Tvpe  Metallization  -  Metallization 

The  pre-metal  dip  is  used  to  remove  any  oxides  formed  on  the  substrate  surface  formed  during  the 
nitride  etch  or  by  general  oxidation  over  time.  The  dip  will  slowly  etch  the  SiN,  thus  a  thorough  DIW 
rinse  is  critical  to  ensure  the  SiN  remains  intact. 

Samples  should  be  loaded  into  the  metallizer  as  soon  as  possible  following  the  pre-metal  dip. 

Prep 

Mix  a  pre-metal  dip  solution  consisting  of  1  part  buffered  oxide  etch  and  50  parts  18  Meg  DIW 
Load  Au  and  Zn  into  the  boat  roundtable  at  the  base  of  the  metallizer  enclosure. 

Mark  where  the  metals  are  on  the  rotation  wheel  outside  the  enclosure. 

129. Dip  the  sample  into  the  pre-metal  dip  for  15  seconds. 

130. Flush  with  tap  DIW  for  2  minutes. 

13 1.  Thoroughly  blow  dry  with  nitrogen 

132. B10W  dry  with  nitrogen  some  more.  Do  not  hot  plate  bake. 

133.1mmediately  load  the  sample  into  the  metallizer  and  press  “Cycle”  on  the  control  console. 

134. Pump  system  to  1x10"®  torr  and  use  standard  posted  procedures  to  deposit  50A  of  Au,  200A  of  Zn, 
and  2000A  of  Au. 

135. C00I  the  chamber  for  30  minutes  after  metallization,  then  remove  sample. 
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136. Place  sample  into  beaker  of  acetone  for  10  minutes.  Agitate  with  a  stream  of  acetone  to  remove  metal 
from  substrate.  Do  not  let  acetone  dry  on  substrate,  this  is  DESTROY  you  sample,  guaranteed. 

137. Remove  sample  from  acetone  while  spraying  the  sample  with  acetone.  Place  into  a  clean  beaker  of 
acetone  and  boil  for  5  minutes  at  1 19°C. 

138. Remove  from  acetone,  do  not  allow  sample  to  dry.  Spray  with  acetone  for  several  seconds  as  you 
remove  the  sample  from  the  acetone  bath,  then  spray  with  methanol  to  remove  the  acetone  for  several 
seconds.  DO  NOT  allow  sample  to  become  dry  ^ng  the  solvent  cleaning  process. 

139. Boil  sample  in  methanol  at  119®C  for  5  minutes. 

UO.Remove  from  methanol,  do  not  allow  sample  to  dry.  Spray  with  methanol  for  several  seconds  as  you 
remove  the  sample  from  the  methanol  bath,  then  spray  with  isopropyl  alcohol  (IPA)  to  remove  the 
methanol  for  several  seconds. 

141. Boil  sample  in  IPA  at  119®C  for  5  minutes. 

142. Remove  from  IPA,  do  not  allow  sample  to  dry.  Spray  with  IPA  for  several  seconds  as  you  remove  the 
sample  from  the  bath. 

143.  While  still  covered  with  IPA,  place  on  a  lint  free  rag  and  blow  the  IPA  from  the  sample  with  nitrogen 
giin  Ek)  not  let  the  IPA  evaporate  without  the  nitrogen  blow. 

144. Examine  sample  using  an  optical  microscope.  Ensure  metal  is  removed  from  unwanted  areas.  If 
unwanted  metal  remains,  place  the  sample  into  acetone  ultrasound  for  5  minutes  in  the  Cole-Parmer 
ultrasound  in  the  class  10,000  hood.  DO  NOT  USE  the  small  Fisher  ultrasound.  These  will  breakup 
your  sample,  and  delaminate  the  metal.  Go  to  step  91  following  the  ultrasound. 


P-tvne  Contact  -  Anneal 

Note:  If  you  elect  to  process  the  top  side  first,  omit  this  step. 


Prep 

Hot  plates  at  119®C  and  200®C 


145. Place  the  sample  on  a  clean  glass  microscope  slide  hot  plate  bake  to  1 19®C  for  at  least  1  minute. 

146. Move  the  slide  and  sample  to  a  200®C  hot  plate  and  set  the  hot  plate  for  350®C.  The  plate  will  ramp 
to  350“C  within  20  minutes.  Hold  at  this  temperature  for  20  minutes,  and  set  the  hot  plate  for  100”C. 
The  plate  will  cool  to  100®C  within  20  minutes. 

147. Remove  the  sample  and  cool  for  1  minute. 


Pad  Metal  -  Photoresist 
Prep 

Hot  plate  at  119“C 

Check  availability  of  1813  photoresist.  If  out,  mix  using  28  ml  of  1827  resist  with  10  ml  type  P 
thinner.  Mix  by  rolling  the  bottle  on  the  table  and  gently  rotating  end  to  end  for  5  minutes.  DO  NOT 
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SHAKE  or  ultrasound.  Shaking  will  get  bubbles  into  the  resist;  ultrasound  will  break  up  the  long 
chain  photoresist  polymers. 

Ensure  spinner  rotates  at  4000  rpm. 


148. Place  sample  on  a  clean  glass  microscope  slide,  and  hot  plate  bake  at  1 19°C  for  5  minutes. 

149. Remove  from  hot  plate  and  allow  sample  to  cool  for  1  minute. 

150. Center  sample  on  i^nner  chuck,  engage  vacuum.  Check  vacuum  hold. 

15  l.Test  spin  sample  to  ensure  acceleration,  spin  speed,  and  spin  time.  Sample  should  spin  up  to  speed  in 
1  second,  spin  at  4000  rpm,  and  spin  for  30  seconds. 

152.  With  sample  stopped,  dispense  enough  1813  photoresist  to  cover  sample,  typically  1  to  3  drops,  but 
not  bleed  over  edges.  Immediately  upon  completing  the  dispense,  start  spinner. 

153. Remove  from  spinner  without  distxubing  top  of  resist.  This  can  be  triclg^  so  be  carefiil.  If  the  resist  is 
not  completely  covered  and  free  from  dust  eddies  and  holes,  briefly  rinse  sample  with  acetone  and  go 
to  step  2. 

154. Place  sample  face  up  on  a  clean  microscope  slide  and  hot  plate  bake  for  5  minutes  at  1 19‘’C. 

155. Remove  from  hot  plate  and  cool  for  5  minutes. 


Pad  Metal  -  Exposure 
Prep 

Ensure  lamp  control  power  supply  has  a  number  close  to  274  watts  on  the  display. 

If  the  display  says  cold,  wait  imtil  a  number  is  display. 

Turn  the  MJB-3  exposure  tool  on  1^  pressing  the  power  button 

Set  the  MJB-3  to  HP  (Hi^  Precision)  mode  by  pressing  the  HP/ST  button  until  the  HP  side  lights 
Ensure  the  vacuiun  mask  light  is  off 

Clean  the  sample  holder  chuck  with  small  amount  of  acetone  on  a  wipe  if  necessary. 

Mix  a  solution  of  1  part  351  developer  to  6  parts  water,  locally  25ml  / 150  ml. 

Find  the  mask  you  will  use  and  visually  inspect  to  ensure  mask  is  clean.  If  not,  clean  it. 


156.  With  nitrogen  gun,  blow  surface  of  exposure  tool  to  remove  dust.  Do  not  attach  ma^  yet. 

157.1f  small  amounts  of  photoresist  are  on  the  back  of  the  sample,  place  sample  upside  down  on  the 
sample  chuck  of  the  exposure  tool,  and  e^qwse  for  99  seconds.  If  large  amounts  of  resist  are  on 
sample,  remove  with  acetone  being  careful  not  to  get  acetone  on  the  front  of  the  sample.  If  this 
happens,  however,  rinse  briefly  with  acetone  and  go  to  step  2. 

158.  This  stq)  is  used  to  remove  edge  beading.  Turn  sample  over  and  place  a  clean  mechanical  sample 
over  the  sample  such  that  1  mm  of  two  edges  are  exposed.  Expose  for  60  seconds.  Repeat  procedure 
to  e?^sure  die  two  remaining  edges. 

159. Devel<^  sample  in  developer  for  60  seconds  or  until  all  back  and  edge  photoresist  is  removed.  Use 
caution  not  to  disturb  photoresist  on  the  front  of  the  sample. 
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160.  Rinse  in  Deionized  Tap  Water  (DIW)  for  2  minutes.  Thoroughly  blow  dry  with  nitrogen  gun.  A  hot 
plate  bake  is  not  necessary  at  this  point. 

161.  Attach  mask  to  the  mask  chuck  and  press  vacuum  mask.  This  first  exposure  requires  the  dark  field 
PAD  quadrant  of  the  mask.  The  dull  chromium  side  of  the  mask  must  be  facing  aw^  fi'om  the 
holder.  Insert  holder  into  the  aligner  and  tighten  thumbscrews. 

162. Place  sample  face  up  on  MJB-3  sample  chuck.  Sample  must  cover  the  flush  open  holes  on  the  chuck. 
If  not,  place  small  squares  of  tape  over  open  holes.  Note  that  the  large  depressed  hole  on  the  left  of 
the  chuck  must  not  be  covered  or  blocked.  This  is  the  nitrogen  hole  for  separation  mode. 

163.  While  watching  sample  from  the  side,  or  through  the  microscope,  raise  the  sample  to  the  mask.  The 
sample  should  touch  the  mask  when  the  substrate  raiser  lever  is  at  the  90%  position.  If  the  sample 
doesn’t  contact  at  the  90%  position,  raise  or  lower  the  substrate  chuck  with  the  chuck  dial  gauge.  Be 
careful;  you  can  destroy  samples  and  damage  your  mask  if  the  chuck  is  set  too  high;  start  low  to  be 
safe. 

If  the  sample  jumps  to  the  mask,  and  cannot  be  put  into  separation,  one  or  more  of  the  holes  in 
the  substrate  chuck  are  not  covered  properly.  Lower  chuck  and  check  for  good  coverage. 

If  the  sample  rotates  in  separation  mode,  the  chuck  is  set  just  a  bit  too  high  and  one  comer  is 
catching.  Lower  chuck. 

If  in  contact  mode,  Fresnel  contact  rings  to  not  form,  the  chuck  is  too  low.  Raise  chuck  holder 
until  rings  form  at  the  comers,  but  go  away  in  separation  mode. 


164.  Align  sample  cleaved  edge  to  a  straight  edge  on  the  mask. 

165. Press  the  vacuum  chamber  button  and  ensure  vacuum  pressure  drops  to  -0.  8  bar  on  the  gauge.  If  not, 
ensure  vacuum  line  to  the  UV  exposure  tower  is  off.  Next  check  to  ensure  vacuum  holes  on  the 
sample  chuck  are  covered,  and  finally,  check  to  ensure  vacuum  ring  on  the  sample  chuck  is  good  and 
touches  the  mask  around  its  entire  circumference  to  the  mask  during  contact. 

166.1f  pressure  drops  to  a  good  level,  visually  check  the  alignment  once  again.  Note  that  vacuum 
chamber  presses  the  sample  closer  to  the  mask,  and  sometimes  moves  the  sample. 

167.1f  alignment  is  ok,  ensure  exposure  time  is  set  for  5.0  seconds,  and  press  the  e?qx)se  button. 

168. Remove  sample  from  exposure  tool,  and  place  in  developer  for  30  to  45  seconds.  You  can  see  the 
resist  develop  away.  If  watching  the  sample  in  the  developer,  hold  sample  in  developer  10  seconds 
past  the  point  where  exposed  photoresist  appears  to  be  removed. 

169. Rinse  in  DIW  for  2  minutes,  then  blow  dry. 

170. Examine  with  optical  microscope  with  the  green  filter  in  place.  If  resist  is  not  completely  developed, 
develop  again  for  15 -second  increments  followed  by  a  2  minutes  DIW  rinse  and  blow  dry.  Sometimes 
this  can  compensate  for  a  bad  exposure,  however,  be  picky.  If  the  exposure  isn’t  perfect,  rinse  briefly 
with  acetone  and  go  back  to  step  2. 

ni.Place  sample  on  a  clean  microscope  slide  and  hot  plate  bake  for  5  minutes  at  119'’C.  This  is  critical 
to  ensure  the  resist  does  not  delaminate  during  subsequent  processing. 


Pad  Metal  -  Metallization 

The  pie-metal  dip  is  not  used  in  this  step.  All  sections  of  the  semiconductor  should  be  covered  with 
either  nitride  or  metal  at  this  point. 
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Load  Ti  and  Au  into  the  boat  roundtable  at  the  base  of  the  metallizer  enclosure. 
Mark  where  the  metals  are  on  the  rotation  wheel  outside  the  enclosure. 


172. Place  the  sample  into  the  left  chamber  of  the  PlasmaTherm  790  and  run  process  DESCLIM  for  1 
minute. 

173. Remove  the  sample  from  the  PlasmaTherm,  load  the  sample  into  the  metallizer,  and  press  “Cycle”  on 
the  control  console. 

174. Pump  ^stem  to  IxlO"*  torr  and  use  standard  posted  procedures  to  deposit  lOOA  of  Ti,  and  3000A  of 
Au.  This  metal  combination  adheres  well  to  almost  any  material.  Note  that  Ti  reacts  with  the 
tungsten  boats,  thus  only  3  to  4  Ti  runs  can  be  made  with  a  single  boat. 

175. C00I  the  chamber  for  30  minutes  after  metallization,  then  remove  sample. 

176. Place  sample  into  beaker  of  acetone  for  10  minutes.  Agitate  with  a  stream  of  acetone  to  remove  metal 
from  substrate.  Do  not  let  acetone  dry  on  substrate,  this  is  DESTROY  you  sample,  guaranteed. 

177. Remove  sample  from  acetone  bath  while  spraying  the  sample  with  acetone.  Place  into  a  clean  beaker 
of  acetone  and  boil  for  5  minutes  at  119°C. 

178. Remove  from  acetone,  do  not  allow  sample  to  dry.  Spray  with  acetone  for  several  seconds  as  you 
remove  the  sample  from  the  acetone  bath,  then  spray  with  methanol  to  remove  the  acetone  for  several 
seconds.  DO  NOT  allow  sample  to  become  dry  during  the  solvent  cleaning  process. 

179. Boil  sample  in  methanol  at  119®C  for  5  minutes. 

150. Remove  from  methanol,  do  not  allow  sample  to  dry.  Spray  with  methanol  for  several  seconds  as  you 
remove  the  sample  from  the  methanol  bath,  then  spray  with  isopropyl  alcohol  (IPA)  to  remove  the 
methanol  for  several  seconds. 

151. Boil  sample  in  IPA  at  119®C  for  5  minutes. 

182. Remove  from  IPA,  do  not  allow  sample  to  dry.  Spray  with  IPA  for  several  seconds  as  you  remove  the 
sample  from  the  bath. 

183.  While  still  covered  with  IPA,  place  on  a  lint  free  rag  and  blow  the  IPA  from  the  sample  with  nitrogen 
giin  Do  not  let  the  IPA  evaporate  without  the  nitrogen  blow. 

184. Examine  sample  using  an  optical  microscope.  Ensure  metal  is  removed  from  unwanted  areas.  If 
unwanted  metal  remains,  place  the  sample  into  acetone  ultrasound  for  5  minutes  in  the  Cole-Parmer 
ultrasound  in  the  class  10,000  hood.  DO  NOT  USE  the  small  Fisher  ultrasound.  These  will  breakup 
your  sample,  and  delaminate  the  metal.  Go  to  step  127  following  the  ultrasound. 


This  completes  topside  processing.  The  sample  is  ready  for  cleaving  and  mounting. 
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